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Abstract 


This  paper  documents  a  technique  to  compare  cross¬ 
correlation  parameters  of  binary  sequences  used  for  spread- 

f  w 

spectrum  multiple-access  communication  systems,  by  performing  a 
thresholding  process  on  the  correlation  functions.  The 
performance  of  Maximal  length,  Gold  and  Kasaml  code  sequences  is 
measured  and  analyzed  for  code  lengths  ranging  from  63  to  1023. 
Comparisons  of  optimized  codes  versus  unoptlmlzed  codes  for  each 
type  of  code  sequence  are  analyzed  in  terms  of  the  thresholding 
process.  Comparisons  are  made  of  the  performance  of  code 
sequences  and  code  sets  at  lengths  of  63,  127,  255  and  1023. 

The  software  used  to  analyze  the  codes  is  discussed  in  terms 
of  structure  and  performance,  and  is  included  as  appendices  in 
the  thesis. 


The  results  of  this  investigation  indicate  that  the 


thresholding 
performance . 


[ 


process  can  be  used  to  evaluate  binary  sequence 

^  _ _ 

(  r  V 

'  ' 


I.  Introduction 


Problem 

One  of  the  most  Important  problems  In  spread-spectrum 
multiple-access  communications  is  the  selection  of  periodic 
sequences  which  have  good  cross-correlation  properties.  Previous 
efforts  in  (8,25)  limit  their  evaluation  of  the  code  sequences  to 
the  maximum  peak,  or  a  metric  of  the  peak,  of  the  cross- 
correlation  functions.  The  goal  of  this  thesis  effort  is  to 
investigate  the  parameters  used  to  evaluate  cross-correlation 
properties  and  to  develop  a  robust  technique  which  can  be  used  to 
evaluate  cross-correlation  properties.  To  establish  a  baseline, 
calculations  of  parameters  of  known  periodic  sequences  are 
compared.  Maximal  length.  Gold  and  Kasami  Code  Sequences  are 
used  to  do  this. 


Scope 

The  thesis  effort  is  restricted  to  the  phase-coded  spread- 
spectrum  multiple-access  systems  as  modeled  in  Chapter  II.  The 
known  code  sets,  used  for  comparison,  are  limited  to  lengths  less 
than  or  equal  to  1023. 


Sequence  of  Presentation 

The  order  of  presentation  in  this  thesis  parallels  the 
approach  to  the  problem.  Chapter  II  is  a  review  of  the 
literature  and  presents  the  background  material  necessary  to 
understand  the  code-division  multiple-access  model,  the 
correlation  functions,  the  code  analysis  algorithm  and  the 
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construction  of  Maximal  length.  Gold  and  Rasaml  code  sequences. 
Chapter  III  gives  an  overview  of  the  software  developed  In  this 
thesis  effort.  Chapter  IV  presents  the  analysis  of  the 
performance  of  the  code  analysis  algorithm  for  varying  lengths  of 
Maximal  length.  Gold  and  Kasaml  codes  and  comparisons  of  these 
code  sets  are  made.  Chapter  V  presents  the  actual  performance  of 
the  software  developed  In  the  thesis.  Chapter  VI  presents  the 
conclusions  and  recommendations  for  further  study.  Appendix  A 
Includes  all  the  perfomance  plots  of  the  code  performance 
algorithms.  Appendices  B  -  E  Include  all  the  documented  software 
used  to  generate  the  data  needed  to  obtain  the  performance  plots. 
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Spread  Spectrum  Communications  Is  the  term  used  to  describe 
a  class  of  transmission  techniques,  in  which  the  actual  signal 
bandwidth  Is  much  greater  than  that  which  is  actually  required  to 
pass  the  information  accross  the  channel  (10:855),  and  the 
transmitted  bandwidth  is  not  a  function  of  the  information 
bandwidth  (1:1).  The  modulated  signal  bandwidth  should  be  10  - 
1000  times  greater  than  the  Information  bandwidth  (2:1).  Spread 
Spectrum  systems  are  divided  into  four  categories:  1)  Direct 
Sequence  Modulation;  2)  Frequency  hopping;  3)  Chirp  Modulation; 
4)  Time  Hopping  and  combinations  of  these  four  (2:1). 

In  Direct  Sequence  modulation  or  pseudonoise  modulation  a 
digital  code  sequence'  is  used  to  modulate  a  carrier.  (1:2)  This 
code  sequence  must  have  a  chip  rate  much  faster  than  the  bit  rate 
of  the  information  sequence.  This  operation,  in  effect,  spreads 
the  information  over  a  larger  bandwidth.  The  receiver  uses  the 
same  spreading  code  sequence  to  despread  the  signal  to  obtain  the 
original  information  signal  which  can  then  be  demodulated  as  if 
the  the  direct  sequence  modulation  had  never  existed  on  the 
carrier  (2:2). 

In  Frequency  Hopping,  a  carrier  frequency  is  shifted  through 
a  large  bandwidth  as  dictated  by  a  digital  code  sequence  (1:2). 
At  any  one  instant  of  time,  only  one  narrow  band  signal  is  trans¬ 
mitted  within  the  broad  spectrum  allocated  for  the  system.  This 
contrasts  sharply  with  the  Direct  spread  system  for  it  uses  the 
entire  spectrum  at  any  instant  of  time. 


Worst  Case  Analysis 


The  discrete  periodic  correlations  functions,  9(jl)  and  9(Jl), 
as  defined  by  Massey  and  Uhran  (8:539),  are  undefined  for  non 
integer  values  of  i.  .  If  we  redefine  the  correlation  functions  as 


f o 1 lows : 

e(T)  »  9(Jl)  +  [  9(A)  +  0(A+1)  ]  (t-A) 

9(x)  =  9(A)  +  [  9(A)  +  0(A+1)  ]  (t-A) 

where  A  =  the  integer  part  of  x  chips 
and  0  <  X  <  L  chip  times 

then  the  discrete  correlation  functions  are  defined  for  the 
entire  region  x  e[0,  L]. 

By  the  same  rationale: 

RIGHT(x)  -  RIGHT(Jl)  +  [RIGHT(A)  +  RIGHT  (  A  +  1 )  ]  (  T -A  ) 

LEFT(x)  -  LEFT(A)  +  [LEFT(a)  +  LEFT(A+1 ) ] ( t-A ) 

where  A  ■  the  Integer  part  of  x  chips 
and  0  <  X  <  L  chip  times 

Mow  the  continuous  aperiodic  correlation  function  can  be 

redefined  as: 


R^^^j(x)  -  RIGHT(x) 
R^^,l(x)  -  LEFT(x) 


And  the  periodic  correlation  functions  can  be  redefined  as: 

9(x)  -  RIGHT(x)  +  LEFT(x) 

9(x)  -  RIGHT(x)  -  LEFT(x) 


Subtltuting  equations  2  and  3  into  equation  1,  yields: 

Vj^^(Xj)  «  Ajdj^_^RIGHT(xj)  +  dj^QLEFT(xj)]cos(9j-a)^Xj)  ( 

It  is  obvious  that  v,  ,(x,)  is  maximized  when  cos(9,-{o  x,)  “  1* 

J  » i  J  J  c  J 
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correLattoti  functions.  The  actual  odd  and  even  periodic 
functions  are  continuous  linear  functions  between  integer  values 
of  chip  times  in  the  region  (0,L),  where  L  is  the  length  of  the 
code  sequences  in  chips.  The  only  discontinuities  occur  at  the 
integer  values  of  chips.  It  should  be  noted  that  the  values 
obtained  using  this  procedure  are  unnormalized  values  that  are 
functions  of  the  length  of  the  codes  The  actual  odd  and  even 
continuous  correlation  functions  are  dlplcted  in  Figure  II-4. 


Figure  11-4.  Example  Continuous  Odd  and  Even  Correlation 


then  the  discrete  even  correlation  value  at  that  Jt ,  can  be  deter¬ 


mined  by: 

1)  Determining  the  number  of  aggreements  and  disagreements  in 

Region  1. 

2)  Determine  the  LEFT  Correlation  "  number  of  agreements  - 

number  of  disagreements  in  Region  1. 

3)  Repeat  steps  1  and  2  for  RIGHT  correlation  in  Region  2. 


Then  the  discrete  even  periodic  correlation  is: 

9(JI)  -  RIGHT(JI)  +  LEPT(jl) 

and  the  discrete  odd  periodic  correlation  is  Just: 

e(JL)  -  RIGHT(jt)  -  LEFT(JI) 

To  illustrate  this  procedure  using  p^  and  p^i  let  the  delay  of  d^ 
be  J1  «  3  chip  times.  Then  the  agreements  and  disagreements  for 
the  regions  are  shown  in  Figure  11-3. 
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Figure  II-3.  Discrete  Correlation  Example 


It  can  be  clearly  seen  from  Figure  II-3,  that  the  values  for  the 

procedure  given  a  time  delay  of  X  are  as  follows: 

LEFT  -  -1 
RIGHT  -  -1 
0(3)  -  -2 
9(3)  -  0 

It  can  be  shown  that  both  of  the  discrete  correlation  functions 
define  endpoints  to  line  segments  of  the  continuous  periodic 
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If  the  double  frequency  terms  are  removed  by  filtering  (a  good 


assumption  because  is  always  much  greater  than  1/T  in  spread 
spectrum  multiple  access  systems)  then  the  output  of  the  i th 
correlator  becomes: 


J*1 


T 

+  /  n(t)pj(t)cos(a)^t)  dt 


When  evaluating  code  sets  for  their  use  in  CDMA  systems,  the  term 
inside  the  summation  is  most  important  because  that  is  where  the 
contribution  of  the  J th  signal  to  the  ith  correlator  exists. 
This  term  provides  the  amount  of  Interference  from  the  J th  user 
in  the  1th  user's  communication  channel.  It  will  be  defined  as 
this : 


(tj) 


cos(ej-a)^Tj) 


Discrete  Odd  and  Even  Correlation  Functions 

For  mathematical  purposes,  if  the  spreading  sequences  which 
are  bipolar  NRZ  signals,  are  mapped  from  [->1,1]  into  [1,0],  then 
the  corelation  of  two  codes  can  be  done  using  modulo-2  arithmetic 
with  equlvllent  results.  The  use  of  this  isomorphism  makes  the 
problem  easier  to  present  as  well  as  easier  to  apply  to  any 
algorithm.  As  an  example,  given  the  two  code  sequences  and  2^^ 
represented  in  vector  form  as: 

“  [1,0, 1,1, 0,1, 1,1] 

£2  ■  10.1.0.1,1.0.1.11 

and  given  a  delay,  I,  equal  to  an  Integer  amount  of  chip  times. 


then 


I  6 


'‘j.o  ■  '*3.1 

and  if  d.  »  ^  d.  .  then 
J  f  V  J  f  1 

'’3,i<^3'..x  ■  -  '‘j,o®<^])  <»> 

Normally,  encoded  data  sequences  have  the  property  that  Pr(dj~l) 


Pr(dj  -  -1)  -  .5,  then 


''j  ,i^‘'j^max 


teCtj) 


±e(Tj) 


Concept  of  the  Code  Analysis  Algorithm 

Figure  II-5  displays  the  v,  -(t.)  for  d,  -  ■  1,  using  the 

j , 1  j  max  j , 0 

example  presented  earlier,  with  a  threshold  set  at  1.  If  Is 
considered  to  be  uniformly  random  between  the  Interval  0  -  L 

chips,  then  a  threshold  as  displayed  In  the  figure  can  be  ad¬ 
justed  to  characterize  this  function.  Using  the  previous  example 
at  a  threshold  ~  1,  the  amount  of  the  correlation  function  which 
exceeds  the  threshold  Is  depicted  as  the  shaded  region  on  the 
figure.  It  Is  computed  using  similar  triangle  relationships  to 
be  a  value  of  0.5  when  normalized  by  the  length  of  the  code.  If 
the  threshold  moved  from  0  to  L  In  the  unnormalized  function,  a 
sat  of  data,  Threshold  vs  Amount  Exceeding  Threshold,  can  be 
easily  obtained.  To  Include  the  effect  of  the  sign  of  function 
in  Equations  7  and  8,  a  second  threshold  of  equal  magnitude  and 
opposite  sign  can  be  adjusted  accordingly,  to  obtain  a  similar 
data  set.  If  the  two  data  sets  are  summed,  then  the  Threshold  vs 
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The  discrete  periodic  correlation  functions,  6(X)  and  6(X), 
as  defined  by  Massey  and  Uhran  (8:539),  are  undefined  for  non 
Integer  values  of  Jl .  If  we  redefine  the  correlation  functions  as 
follows : 

0(t)  -  e(X)  +  [  0(1)  +  0(1+1)  ]  (t-1) 

0(t)  -  0(1)  +  [  0(1)  +  0(1+1)  ]  (r-l) 

where  1  >  the  Integer  part  of  x  chips 
and  0  <  X  <  L  chip  times 

then  the  discrete  correlation  functions  are  defined  for  the 
entire  region  x  e[0,  L]. 

By  the  same  rationale: 

RIGHT(x)  -  RIGHT(l)  +  [rIGHT(I)  +  RIGHT(1+1) ] (x-1) 

LEPT(x)  -  LEPT(l)  +  [lEFT(1)  +  LEPT(1+1) ] (x-1) 

where  1  >  the  Integer  part  of  x  chips 
and  0  <  X  <  L  chip  times 


Now  the  continuous  aperiodic  correlation  function  can  be 
redefined  as: 

R^^^j(t)  -  RlGHT(x)  (2) 

Rj^,i(T)  -  LEPT(x)  (3) 

And  the  periodic  correlation  functions  can  be  redefined  as: 

0(x)  -  RIGHT(x)  +  LEFT(x)  (4) 

0(x)  -  RIGHT(x)  -  LEPT(x)  (5) 

Snbtitutlng  equations  2  and  3  Into  equation  1,  yields: 

Vj^j(Xj)  -  [dj^_jRIGHT(Xj)  +  dj^QLEFT(Xj)]co8(0j-a)^Xj)  (6) 


It  is  obvious  that  v,  ,(t.)  Is  maximised  when  cos(0.-(i)  x,) 

J » *  J  j  c  j 
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then  the  discrete  even  correlation  value  at  that  X,  can  be  deter¬ 
mined  by: 

1)  Determining  the  number  of  aggreements  and  disagreements  In 

Region  1. 

2)  Determine  the  LEFT  Correlation  >  number  of  agreements  - 

number  of  disagreements  In  Region  1. 

3)  Repeat  steps  1  and  2  for  RIGHT  correlation  In  Region  2. 

Then  the  discrete  even  periodic  correlation  Is: 

0(2)  -  R1GHT(2)  +  LEFT(2) 

and  the  discrete  odd  periodic  correlation  Is  Just: 

0(2)  -  RIGHT(2)  -  LEFT(2) 

To  Illustrate  this  procedure  using  p^  and  let  the  delay  of 

be  2  ■  3  chip  times.  Then  the  agreements  and  disagreements  for 
the  regions  are  shown  In  Figure  11-3. 
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Figure  11-3.  Discrete  Correlation  Example 

It  can  be  clearly  seen  from  Figure  11-3,  that  the  values  for  the 

procedure  given  a  time  delay  of  2  are  as  follows: 

LEFT  -  -1 
RIGHT  -  -1 
g(3)  -  -2 
0(3)  -  0 

It  can  be  shown  that  both  of  the  discrete  correlation  functions 


define  endpoints  to  line  segments  of  the  continuous  periodic 


If  the  double  frequency  terms  are  removed  by  filtering  (a  good 


assumption  because  Is  always  much  greater  than  1/T  In  spread 
spectrum  multiple  access  systems),  then  the  output  of  the  1th 
correlator  becomes  (letting  *  1,  1  "  1,2, 
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1« 


T 

+  /  n(t)p,(t)cos(<i>  t)  dt 
0  J  c 


When  evaluating  code  sets  for  their  use  In  CDMA  systems,  the  term 
Inside  the  summation  Is  most  Important  because  that  Is  where  the 
contribution  of  the  jth  signal  to  the  1th  correlator  exists. 
This  term  provides  the  amount  of  Interference  from  the  jth  user 
In  the  1th  user's  communication  channel.  It  will  be  defined  as 
this: 

Discrete  Odd  and  Even  Correlation  Functions 

For  mathematical  purposes.  If  the  spreading  sequences  which 
are  bipolar  NRZ  signals,  are  mapped  from  [-1,1]  Into  [1,0],  then 
the  correlation  of  two  codes  can  be  done  using  modulo-2  arith¬ 
metic  with  equlvllent  results.  The  use  of  this  Isomorphism  makes 
the  problem  easier  to  present  as  well  as  easier  to  apply  to  any 
algorithm.  As  an  example,  given  the  two  code  sequences  and  £2 
represented  In  vector  form  as: 

£1  ■  [1,0, 1,1, 0,1, 1,1] 

£2  -  [0.1,0, 1,1. 0,1,1] 

and  given  a  delay,  1,  equal  to  an  Integer  amount  of  chip  times. 


r(t)  •  n(t)  +  Aj^d^(t)pj^(t)coB(«^t)  + 


K 

I 

j-1 

3^1 


Aj  Pj(t  -  dj(t  -  Tj)  C08(o,^t  +  Gj  - 


where  A.  *■  amplitude  of  the  1th  carrier 

d^(t)  ■■  polar  NRZ  message  signal  of  the  1th  user 
p^(t)  ■  polar  NRZ  spreading  sequence  of  the  1th  user 
6^  ■  phase  of  the  1th  carrier  varying  from  0-2*  uniformly 
*  time  delay  of  the  1th  user  varying  from  0  -  T 
uniformly 

T  "  data  symbol  length  In  time  duration 
n(t)  "  additive  white  Gaussian  noise  with  two-sided 
spectral  density  N  /2 


If  the  received  signal  Is  Input  to  the  Integrate  and  dump  cor¬ 
relator  matched  to  the  1th  signal  s^Ct),  then  the  output  of  this 
correlator  at  time  T  la  given  by: 

T 

Z.  •  f  r(t)p,(t)  cosCb)  t)  dt 
1  0  *  c 


To  Illustrate  what  la  going  on  In  the  system  at  this  point  refer 
to  Figure  II-2  where  the  top  diagram  depicts  the  two  data  stream 
bits,  d^  Q  and  d^  All  of  the  other  diagrams  In  this  figure 
are  In  time  reference  to  the  top  diagram  which,  as  mentioned 
earlier.  Is  In  synchronization  with  the  1th  correlator.  Diagram 
5  depicts  the  kth  data  stream  which  Is  delayed  by  which  for 
this  example  Is  equal  to  3  chip  times  of  the  spreading  sequences. 
Diagram  4  depicts  the  product  of  the  1th  and  kth  spreading 
sequences,  P|^(t-X|^)p^(t) ,  and  Is  segmented  In  time  Into  two 
regions:  Region  1,  from  t  ■■  0  to  t  "  x^,  and  Region  2,  from  t  " 
X|^  to  t  "  T.  The  Integrations  of  these  two  regions  In  time 
produce  the  continuous-time  partial  cross-correlation  functions: 


Figure  II-l.  Phase-coded  CDMA  System  Model  (12:797) 


The  data  signal  Is  modulated  on  the  phase  coded  car¬ 

rier  c^(t)  to  produce  the  transmitted  signal  Sj^(t),  each  given 
by:  Cj^(t)  ■  cosCw^t  +  0^^) 

Sj^(t)  ■  Aj^  P|^(t)  d|j(t)  cos(w^t  +  0j^) 

The  received  signal  r(t)  Is  given  by 

)  d.(t  -  T  )  cos(to  t  +  0.  -  «  T  ) 
j  j  c  j  c  3 

If  It  Is  assumed  that  the  1th  correlator  Is  synchronized 
with  the  1th  signal,  l.e.  that  •  0  and  the  phases  are 
referenced  to  signal  1,  l.e.  0^  >  0,  then  the  received  signal 
at  the  1th  correlator  becomes: 


r(t)  ■  n(t)  + 


J-1 


-  ■'j 


situations  where  Interference  and  multipath  problems  exist,  there 
are  significant  degradations  In  both  the  FDMA  and  TDMA  systems. 
For  reasons  such  as  these,  code  division  multiple  accessing 
schemes  have  evolved.  The  CDMA  systems  allow  system  users  to 
transmit  simultaneously  In  time  and  occupy  the  same  RF  bandwidth 
(1:93).  Generally,  the  CDMA  system  operates  asynchronously  to 
eliminate  system  timing  problems.  However,  the  Imperfect  ortho¬ 
gonality  of  the  different  spreading  sequences  reduces  the  capa¬ 
city  of  the  CDMA  system.  An  important  issue  confronting  system 
designers  Is  design  of  sets  of  spreading  sequences  that  operate 
efficiently  In  this  system.  The  correlation  properties  of  these 
sets  of  sequences  have  a  major  Impact  on  the  capacity  of  this 
type  of  system  and  Its  overall  performance. 

The  Code  Division  Multiple  Access  System  Model 

The  basic  model  of  the  Code  Division  Multiple  Access  system 
Is  depicted  In  Figure  II-l.  In  this  system,  each  of  the  K  users 
Is  given  Its  own  spreading  code  sequence  which  has  "orthogonal 
like"  properties  with  the  other  assigned  users.  Because  of  the 
fact  that  any  user  can  Initiate  or  terminate  Its  traffic  at  any 
time  In  an  asynchronous  system.  It  Is  more  difficult  to  design 
and  evaluate  different  code  sets  for  use  In  the  system.  It  Is 
deslreable  to  obtain  code  sets  with  good  cross-correlation  and 
partial-correlation  properties  (12:795).  This  fact  will  become 
evident  as  the  model  Is  explained  In  the  following  pages. 
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Chirp  Modulation  is  most  commonly  used  in  radar,  which 
transmits  swept  frequency  pulses.  The  receiver  uses  a  dispersive 
filter  to  compress  the  signal  into  a  narrower  time  slot,  so  it 
behaves  in  the  same  way  a  high  power  pulse  does  (2:3). 

In  Time  hopping,  the  transmission  time  and  period  is 
controlled  with  the  digital  code  sequence.  When  used  by  itself. 
Time  hopping  can  be  viewed  as  pulse  modulation  under  code  se¬ 
quence  control  (2:3). 

Spread  Spectrum  systems  may  offer  any  or  combinations  of  the 
following  benefits:  l)  Simultaneous  use  of  the  spectrum  by 
multiple  users;  2)  Interference  rejection  caused  by  intentional 
and  unintentional  sources  such  as  Jamming  and  multipath;  3)  Low 
density  output  signals  providing  a  lower  probability  of  intercep¬ 
tion;  4)  High  resolution  ranging;  5)  Accurate  universal  timing; 
and  6)  Inherent  message  privacy  (1,2). 

The  disadvantages  of  using  spread  spectrum  are  the  Increased 
complexity  of  the  system  and  Increased  difficulty  in  allocation 
of  frequencies  (2:4). 

The  effort  of  this  thesis  is  concentrated  in  the  Direct 
Sequence  Spread  Spectrum  area  and  the  advantage  of  multiple 
access  communications.  Multiple  access  communications  systems 
allow  simultaneous  access  to  the  channel  to  many  users.  The 
three  most  common  multiple  access  techniques  are  time  division 
multiple  access  (TDMA),  frequency  division  multiple  access 
(FDMA),  and  code-division  multiple  access  (CDMA).  TDMA  and  FDMA 
separate  the  signals  in  time  and  frequency  domains  respectively. 
CDMA  employs  spread  spectrum  techniques  (SSMA) ,  and  provides 
another  way  to  separate  the  signals  at  the  receiver  (10:868).  In 


1  then 


”  ^J.o  ■  “j,!-  *i  ■ 

and  if  d.  ^  ^  ^4  i>  ^  then 

J  >  ^  J  >  1  ^ 

-  I-EFTCtj)!  -  ij,o9(^j>  (*> 

Normally,  encoded  data  sequences  have  the  property  that  Pr(d^*l) 


Pr(dj  *  -1)  =  .5,  then 


±6( Tj ) 


±0 ( Tj ) 


Concept  of  the  Code  Analysis  Algorithm 

Figure  11-5  displays  the  v,  ,(t,)  for  d,  ^  ■  1,  using  the 
^  j  »i  j  max  J , 0 

example  presented  earlier,  with  a  threshold  set  at  1.  If  Tj  Is 
considered  to  be  uniformly  random  between  the  Interval  0  -  L 
chips,  then  a  threshold  as  displayed  In  the  figure  can  be  ad¬ 
justed  to  characterize  this  function.  Using  the  previous  example 
at  a  threshold  *  1,  the  amount  of  the  correlation  function 
exceeding  the  threshold  Is  depicted  as  the  shaded  region  on  the 
figure.  It  Is  computed  using  similar  triangle  relationships  to 
be  a  value  of  0.5  when  the  length  of  the  code  Is  normalized.  If 
the  threshold  moved  from  0  to  L  In  the  unnormalized  function,  a 
set  of  data.  Threshold  vs  Amount  Exceeding  Threshold,  can  be 
easily  obtained.  To  Include  the  effect  of  the  sign  of  function 
In  Equations  7  and  8,  a  second  threshold  of  equal  magnitude  and 
opposite  sign  can  be  adjusted  accordingly,  to  obtain  a  similar 
data  set.  If  the  two  data  sets  are  summed,  then  the  Threshold  vs 
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Spreading  Code  Sequence  Generation 

Pseudonoise  (PN)  sequences  or  Linear  Feedback  Shift  Register 
(LFSR)  sequences  are  easily  generated  at  the  transmitter  and  the 
receiver.  The  LFSR,  as  shown  in  Figure  II-7  is  composed  of 
binary  storage  elements,  taps  for  reading  storage  contents,  and  a 
device  that  performs  the  modulo-2  addition  of  the  tapped  storage 
elements.  The  binary  storage  elements  must  be  initially  loaded 
with  O's  and  I's  (not  all  0).  As  the  LFSR  is  clocked,  the 
modulo-2  addition  device  generates  the  output  of  the  code  gener¬ 
ator  as  well  as  the  input  to  the  first  binary  storage  location. 


Figure  II-7.  Example  Linear  Feedback  Shift  Register 


This  LFSR  as  also  known  as  a  Simple  Shift  Register  Generator 
(SSRG)  because  all  feedback  signals  are  returned  to  a  single 
input(l:30).  If  the  initial  loading  of  the  generator  is  all  O's 
with  the  exception  of  the  last  position,  then  the  output  sequence 
will  not  have  any  transient  bits.  The  general  form  of  the  LFSR 
is  depicted  in  Figure  11-8  which  has  N  storage  locations  and  taps 
feeding  into  one  summation  device. 
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Figure  II-8.  General  Linear  Feedback  Shift  Register  (1:30). 


The  fundamental  linear  recursion  relation  which  provides  the  code 
sequence  element  Is  written  as: 

"Jj  ‘i\-i 


Maximal  Sequences 

Maximal  length  sequences  are  produced  using  the  LFSR  In 
Figure  11-8  using  a  primitive  polynomial  of  degree  N  to  designate 
which  tap  positions  will  be  used.  They  have  a  maximal  length 

period  of  2  -  1  chips.  The  maximal  sequences  have  the  following 

properties : 

N*  1  N* 1 

1)  Balance  property  -  There  are  2  ones  and  2  -1  zeroes 

2)  Run  property  -  In  any  period,  half  of  the  runs  of  I's  or 

O's  Is  of  length  1,  1/4  are  of  length  2,  1/8  are  of 

length  3  etc. 

3)  The  autocorrelation  function  6(t)  Is 

0(t)  ■  (  ^  r  ■  0,  L,  2L  ...chips 

<  -1  for  T  <  1  and  x  <  L-1  chips 

(  (x-Jl)(L-l)  for  -1  <  X  <  1,  A  - 

Integer  part  of  x  for 
each  period  of  the 
function  (refs  1,4,10). 
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Gold  Code  Generation 


'.N 


Gold  codes  can  be  generated  using  one  period  of  the  output 
sequences  of  two  linear  feedback  shift  registers,  which  have  tap 
assignments  designated  by  preferred  pair  polynomials.  The 
preferred  pair  must  be  of  degree  N,  such  that  N  is  not  divisible 
by  4.  If  the  output  sequences  are  modulo-2  added  to  one  another 
for  each  possible  phase  difference,  then  L  different  Gold  se¬ 
quences  are  produced.  If  the  original  two  sequences  are  included 
then  this  process  yields  a  set  of  L+2  Gold  sequences  (ref 
1,4,5,11).  To  illustrate  by  the  following  example: 

5  2 

f^(x}  "x  -fx  -1-1  Is  a  primitive  polynomial  of  degree  5 

To  obtain  the  preferred  pair.  Gold's  theorem  states  that: 

Given  f.(x)  is  a  primitive  polynomial  of  degree  N  such  that 
N  is  not  divisible  by  4  then  let 


j  1  +  h  .ad 

'  (  1  .  H  ...n 

then  f^Cx)  is  the  minimal  polynomial  of 
In  this  example  t  ■  1  +  2^^  l)/2  _  ^ 

The  minimal  polynomial  of  of  degree  5  is: 
t^ix)  -  x^  +  x^  +  x^  +  X  +  1 

Gold's  theorem  also  states  that  using  a  preferred  pair  of  primi¬ 
tive  polynomials,  of  degree  N,  whose  LFSR  generator  sequences  are 
maximal  length  with  period  L  -  2^-  1,  to  produce  a  set,  C,  of 
L  +  2  sequences  of  length  L  ~  2  -  1.  The  discrete  even  cross 
correlation  function: 


|0^^j(k)|  < 


1  +  N  odd 

1  + 
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and  9.  .(k)  «*  one  of  the  three  [-l,-t,t-2] 

for  k  an  integer  In  [0,Ll 

(refs  1,4,5 ,11) . 

For  this  example,  |0.  ,(k)|  <  1  +  >  9^  and  the  three 

i  *  J 

values  that  it  obtains  are  [-1,-6,  8]. 

Kasami  Code  Generation 

A  set  of  Kasami  codes  can  be  generated  using  the  maximal 

length  sequence  generated  in  the  LFSR  of  even  length.  Using  the 

N/  2 

decimation  of  this  sequence  by  2  '  1  to  produce  a  second 

N  /  2 

sequence  with  length  2  -  1.  A  third  sequence  of  length  L, 

N/  2 

period  2  -  1,  is  produced  by  the  second  sequence.  The 

modulo-2  summation  of  the  first  and  third  sequence  through  the 
N/  2 

2  -  2  different  phase  shifts  of  the  third  sequence  will  pro- 

N/  2 

duce  a  set  of  2  ^  -  1  sequences.  If  the  first  sequence  is 

included  in  the  set,  then  the  entire  Kasami  set  generated  will 
N/  2 

have  2  '  sequences  (11:569).  The  following  example  will  il¬ 
lustrate  this  procedure:  Given  the  maximal  sequence  generated  by 
the  primitive  polynomial  x^  x  -f  1  is 

U1 11 1010^011001 10U1011010£10011100£10111100i01000110£00100000 
the  decimation  by  9  yields  the  second  sequence 

1110010 


i. 


V. 


the  third  sequence  is  made  by  the  second  to  yield 

111001011100101110010111001011100101110010111001011100101110010 


'■j^”Vf>ji '"jiinj*  JM  ' 


sequence  for  a  total  of  8  sequences.  Like  the  Gold  Sequences, 
the  discrete  even  autocorrelation  and  cross-correlation  are  3 
valued.  However  this  set  has  [-1,  +  1),  2^^^  -  1]  and 
for  this  example  have  values  [-1,  -9,  7].  Hence  the  discrete 
even  cross-correlation  function  is 

|Q  ^  2**^^  +  1  N  even  only 
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III.  Software  Overview 


This  chapter  provides  a  general  discussion  of  the  software 
designed  In  the  thesis,  In  the  order  as  the  fully  documented 
programs  appear  In  Appendices  B,  C,  D,  and  E.  This  discussion  Is 
presented  at  the  structure-function  level  for  understandablllty . 
If  the  reader  desires  more  detailed  discussion  of  each  module 
presented  here,  he  should  refer  to  these  Appendices.  It  should 
be  noted  that  the  software  was  designed  and  tested  using  the 
Turbo-Pascal  Compiler,  version  2.0,  from  Borland  International. 
This  compiler  accepts  all  standard  Pascal  types  and  Identifiers 
and  has  many  non-standard  types  and  Identifiers  available  to  the 
user  as  well.  For  transportability  purposes,  the  software  was 
written  such  that  few  nonstandard  types  are  used.  When  they 
were  used,  a  comment  Identifying  its  use,  directly  follows  the 
code  as  It  appears  In  the  appendices.  It  Is  worth  noting  that 
the  EOLN  and  EOF  functions  are  performed  differently  In  many 
compilers,  and  for  this  reason,  any  modules  performing  character 
I/O  may  have  to  be  modified  accordingly.  The  programs  were  run 
using  the  Berkely  Pascal  Compiler,  and  required  additional  proce¬ 
dures  for  string  I/O  because  a  string  Is  a  non-standard  type. 

Program  GenerateCodes 

Figure  III-l  depicts  the  structure  of  the  code  generating 
program,  GenerateCodes.  It  uses  linear  feedback  shift  registers 
(LFSR's)  constructed  from  an  array  of  characters,  I's  and  O's,  to 
generate  Maximal  Length  Codes,  Gold  Codes,  and  Kasaml  Codes.  The 
outputted  sequences  are  stored  In  user  designated  character 
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with  modified  extensions  for  each  code  generated  within  a 


files 

code  set  (i.e.  .001,  .002,  for  the  first  and  second  code  gener¬ 
ated  in  a  code  set).  The  main  module  queries  the  user  to  -deter¬ 
mine  which  type  of  code  to  generate,  and  calls  one  of  three 
modules;  LinearGenera tor ,  Go IdCodeGenera tor  ,  and  KasamlGenerator . 
Each  of  these  performs  most  of  the  user  I/O  required  to  generate 
the  codes.  They  each  make  calls  to  Loadlnput  and  Test,  which 
make  tap  assignments  and  validate  the  inputted  filenames  respect¬ 
ively.  The  GoldCode  Generator  and  the  Kasami  Generator  use  the 
function  GetRidOf  to  free  the  memory  of  linked  lists,  represent¬ 
ing  code-lists,  that  each  generator  no  longer  requires  for  gener¬ 
ation. 

GenLinear .  The  GenLlnear  module  is  called  by  the  LinearGen- 
erator  after  the  input  has  been  validated.  It  loads  an  array, 
representing  the  LFSR,  with  all  O's  and  a  1  in  position  C^,  as 
described  in  Chapter  II.  It  calls  module  ModSumTaps  to  do  the 
modulo-2  addition  of  the  array  elements,  identified  by  the  tap 
assignments,  to  produce  the  output  bit,  CK.  CK  is  written  to  the 
code  file  and  is  passed  to  the  module  Shlf tRegis ter ,  which  shifts 
the  contents  of  the  LFSR  by  one  position,  and  loads  CK  into  the 
C|^_^  position.  This  procedure  is  continued  until  2  -  1  bits 
have  been  generated. 

GenGo Id .  The  GenGold  module  is  called  by  the  GoldCode 
Generator  after  the  tap  assignments  have  been  made  for  two  LFSR's 
and  the  filename  has  been  validated.  The  user  identifies  the 
number  of  sequences  he  wants  generated  from  1  to  2  +1.  It 
produces  two  LFSR  sequences  in  the  same  manner  as  the  GenLlnear 
module,  using  modules  ModSumTaps  and  Shif tRegl s ter ,  but  places 


the  two  output  sequences  into  linked  lists  in  memory,  using 
module  Enque.  The  first  two  sequences  generated  are  written  to 
code  files  with  extensions  .001  and  .002  using  the  Writelist 
module.  The  third  sequence  is  generated  by  calling  module  Mod- 
SumList  which  performs  a  modulo-2  addition  of  the  two  lists  and 
outputs  the  result  into  the  third  file.  Module  ShiftList  is 
called  to  perform  a  cyclic  shift  of  the  first  code  list.  The 
forth  through  the  2  -f  1  sequence  is  generated  by  calling  modules 
ShiftList  and  HodSumList. 

GenKasami .  The  GenKasami  module  is  called  when  the  module 

Kasami-  Generator  has  made  the  tap  assignments  for  one  LFSR  and 

validated  the  inputted  filename.  It  generates  an  LFSR  sequence 

in  the. same  manner  as  the  GenLlnear  module  but  stores  the  output 

sequence  in  a  linked  list  in  memory.  It  calls  module  Decimate  to 

perform  a  decimation  of  the  LFSR  sequence  by  a  factor  N  2 , 

and  produce  a  second  linked  list  in  memory,  using  module  Enque. 

The  second  list  is  then  used  to  produce  a  third  list,  with  period 

of  the  second  list  and  length  of  the  first  list,  using  module 

Enque.  The  first  code  sequence  is  generated  by  writing  the  first 

list  to  file,  with  the  extension  .001,  using  module  Vritelist. 

Successive  sequences  are  generated  using  ModSumList  and  Shift- 

List,  as  described  previously.  Depending  on  the  user's  input, 

N  /  2 

the  GenKasami  module  will  produce  from  1  to  2  sequences. 


Program  PartialCorrelate 

The  PartialCorrelate  program,  as  depicted  in  Figure  I1I-2, 


will  perform  correlations  on  user  designated  code  files  composed 
of  ASCII  I's  and  O's.  The  code  files  are  read  into  linked-lists 


Partial 

Correlate 


Figure  111-2.  Structure 


using  the  module  CreateLlst.  The  program  has  many  options  avail¬ 
able  to  the  user  that  are  presented  in  a  menu  driven  fashion.  In 
all  options  the  actual  discrete  Odd,  and  Even  correlation  func¬ 
tions  are  computed  using  the  Correlate  module,  for  each  integer 
value  of  tau  in  the  interval  [0,L].  The  correlation  function 
values  are  placed  into  an  ordered  linked-list  of  records  using 
the  module  Enque.  If  the  user  desires,  this  list  of  data  can  be 
written  to  file  using  the  module  Writellst.  The  module  ShiftList 
is  only  used  when  the  user  wants  to  evaluate  different  phase 
shifts  of  the  original  codes.  When  called  it  will  shift  both 
code  files  a  designated  number  of  bits  before  correlating.  If 
desired,  the  thresholding  process,  as  described  in  chapter  II, 
can  be  performed  on  the  correlation  functions.  The  user  has  the 
option  of  thresholding  any  or  both  of  the  correlation  functions 
using  a  positive  threshold,  negative  threshold  or  both.  The 
positive  threshold  process  is  performed  in  the  CountUpper  module, 
while  the  negative  threshold  process  is  performed  in  the  Count- 
Lower  module.  The  main  module  calls  the  modules,  CountUpper  and 
CountLower,  and  passes  a  specific  threshold  value  to  obtain  a 
value  for  the  amount  of  the  correlation  function  exceeding  the 
threshold.  The  threshold  values  passed  in,  vary  from  0  to  the 
length  of  the  codes,  and  are  incremented  a  fixed  amount  until  the 
threshold  value  has  reached  the  length  of  the  code  or  the  amount 
exceeding  the  threshold  goes  to  0.  During  each  cycle  of  this 
process,  the  normalized  threshold  value  and  the  amount  of  the 
correlation  function  exceeding  the  threshold  are  written  to  a 
user  designated  file.  The  modules  DisposeOf  and  GetRidOf  are 
used  to  free  the  memory  of  the  code  lists  and  the  correlation 


Lists,  when  their  use  is  no  longer  needed. 


Program  PhaseCorreLate 

The  PhaseCorreLate  program  determines  the  best  and  worst 
phases  of  the  two  inputted  codes  in  terms  of  the  Odd  correLation 
function.  The  structure  of  the  program  is  shown  in  Figure  III-3. 
The  program  reads  the  code  sequences  into  L inked- L i s ts  using  the 
moduLe  CreateList.  The  program  wlLL  perform  correLations  of  the 
user  designated  code  fiLes  using  moduLe  CorreLate,  and  determine 
the  minimum  and  maximum  values  of  the  Odd  function  for  a  certain 
phase  in  moduLe  Mlnlmax.  The  program  repeats  this  process  for 
every  phase  shift  of  the  code  files,  and  uses  module  Enque  to 
place  these  values  into  another  linked-llst.  The  ShiftLlst  mod¬ 
ule  cyclically  shifts  the  code  lists  appropriately  for  each  phase 
prior  to  commencing  the  correlation  process.  After  correlations 
have  been  performed  on  all  phase  shifts  of  the  code  files,  then 
the  module  Maxlmlnlmlze  Is  called.  This  module  will  traverse  the 
llnked-llst  of  minimum  and  maximum  Odd  values  to  determine  the 
maximum  and  minimum  spread  for  all  phases.  The  main  program  will 
then  traverse  the  linked  list  and  write  the  phase,  maximum  and 
minimum  Odd  values  for  each  phase  shift  that  has  a  Maximum  or 
Minimum  spread.  Just  as  before,  the  DisposeOf  and  GetRidOf 
modules  free  the  memory  used  for  the  Data  linked- lis ts ,  and  code 
lists,  when  no  Longer  required  in  the  program.  The  user  has  the 
menu-driven  option  of  performing  off-peak  autocorrelations  or 
total  cross  correLations.  This  program  can  easily  be  modified  to 
find  the  min  and  max  of  the  even  or  odd  correlation  functions  for 
each  phase,  if  desired. 
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detected. 


Comparison  of  Maximal ,  Gold  and  Kasaml  Codes 

Figures  A-73  and  A-74  show  the  maximum  even  and  odd  thresh¬ 
olds  obtained  from  the  Mass-correlation  algorithm  for  all  three 
un-optlmlzed  classes  of  codes  of  length  63.  In  the  threshold 
range  0.15  to  0.4,  the  Kasaml  set  clearly  has  fewer  peaks  and 
sub-peaks  than  either  of  the  other  sets  of  even  correlations.  In 
terms  of  the  even  correlations,  the  Maximal  set  shows  better 
performance  than  the  Gold  set  within  this  range.  The  odd  correl¬ 
ations  are  very  similar  In  this  range,  with  the  Kasaml  set  show¬ 
ing  a  slightly  better  perf oribance. 

Figures  A-75  and  A-76  show  the  maximum  even  and  odd  thresh¬ 
olds  obtained  from  the  Mass-correlation  algorithm  for  the  Linear 
and  Gold  un-optlmlzed  classes  of  codes  of  length  127.  In  terms 
of  the  even  thresholds  in  the  range  0.13  to  0.20,  the  Gold  codes 
show  fewer  peaks  than  the  Maximal  code  set.  The  odd  thresholds 
show  the  odd  correlation  functions  to  be  quite  similar  for  both 
Maximal  and  Gold  code  set  over  the  entire  range  of  threshold 
values . 

Figures  A-77  and  A-78  show  the  maximum  even  and  odd  thresh¬ 
olds  obtained  from  the  Mass-correlation  algorithm  for  the  Linear 
and  Kasaml  un-optlmlzed  classes  of  codes  of  length  255.  In  the 
threshold  range  of  0.07  to  0.2  of  the  even  threshold  functions, 
the  Kasaml  code  set  has  fewer  peaks  than  the  Maximal  set.'  In  the 
odd  threshold  functions,  both  the  sets  are  similar  for  thresholds 
above  0.1.  Because  of  the  equal  importance  of  the  even  and  odd 
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number  of  sub-peaks  in  the  threshold  range  of  0.1  and  0.13.  The 
odd  correlations  however,  show  variations  below  0.32. 

Length  255  Kasami .  The  odd  and  even  cross-correlation  func¬ 
tions  for  Kasami  code  2  and  3,  of  length  255,  are  shown  in 
Figures  A-64  and  A-65.  The  even  correlation  function  is  three 
valued  [15,  -1,  -17]  and  somewhat  irregular  while  the  odd  correl¬ 
ation  is  vary  irregular  with  maximum  peaks  at  +  46.  The  thresh¬ 
old  function  for  both  is  depicted  in  Figure  A-66.  Neither 
threshold  plot  has  significant  values  above  0.13,  but  the  even 
function  shows  more  peaks  below  this  threshold.  The  Mass-correl¬ 
ation  output  of  entire  subset  of  10  un-optlmlzed  Kasami  code 
sequences  is  plotted  in  Figures  A-67  and  A-68.  There  is  some 
variation  in  the  number  of  sub-peaks  region  between  0,06  and  0,07 
for  the  even  correlations  and  the  odd  correlations  show  increas¬ 
ing  variations  below  the  threshold  of  0.15, 

Length  1023  Kasami .  Figures  A-69  and  A-70  show  the  output 
of  the  Mass-correlate  algorithm  for  the  8  un-optlmlzed  Kasami 
codes  of  length  1023.  There  is  no  variation  in  the  even  correla¬ 
tions  above  the  threshold  of  0.03  and  minor  fluctuations  below 
this  value.  The  odd  correlations  show  slight  variations  below 
the  threshold  value  of  0.06. 

Comparison  of  Kasami  Code  Sets ♦  Figures  A-71  and  A-72  show 
the  maximum  even  and  odd  thresholds  obtained  from  the  Mass- 
correlate  algorithm  for  the  Kasami  codes  generated  in  Table  3. 
As  the  length  of  the  code  increases,  the  performance  dramatically 
increases  for  both  functions.  The  even  functions  clearly  out¬ 
performing  the  odd  in  terms  of  the  number  of  peaks  and  sub-peaks 
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Length  63  Kasaml .  Figures  A-54  and  A-55  show  the  even  and 
odd  cross-correlation  functions  codes  2  and  3.  The  even  function 
is  three  valued  [7,  -1,  -9],  for  Integer  values  of  tau.  The  odd 
cross-correlation  function  is  irregular  with  maximum  peaks  of 
+  17.  Figure  A-56  shows  the  threshold  plots  of  this  cross- 
correlation  for  both  functions,  and  obviously  the  odd  function 
has  more  peaks  in  a  range  above  0.11.  It  is  worth  noting  that 
the  apparent  discontinuities  of  the  even  threshold  occur  as  the 
threshold  passes  the  three  values  mentioned  earlier. 

The  off-peak  odd  auto-correlation  of  code  2  and  the  odd 
cross-correlation  of  codes  2  and  3  were  optimized  in  terms  of 
phase  and  maximum  peaks  and  the  resulting  correlation  functions 
are  depicted  in  Figures  A-57  and  A-58.  The  worst  and  best  phase 
of  the  odd  auto  and  cross-correlation  threshold  functions  are 


plotted  in  Figures  A-59  and  A-60.  In  both  cases,  the  odd  correl¬ 
ation  functions  show  a  significant  difference  in  the  number  of 
peaks  and  sub-peaks  in  the  threshold  range  of  0.05  to  0.25.  When 
the  even  and  odd  cross-correlation  are  analyzed  in  Figure  A-60, 


it  is  apparent  that  the  odd  function  has  more  peaks  and  sub-peaks 
in  the  range  between  0.15  and  0.25.  When  the  minimum  odd  and 
even  cross-correlation  functions  are  compared  in  Figure  A-61  the 
even  function  appears  to  provide  more  Interference  in  the  thresh¬ 
old  range  of  0.15  to  0.21. 

Figures  A-62  and  A-63  show  the  output  of  the  Mass-correlate 
algorithm  for  all  8  of  the  un-optimlzed  codes  of  length  63.  The 
even  correlations  show  little  or  no  variation  above  a  threshold 
of  0.13  for  the  entire  set.  and  substantial  variation  in  the 


threshold  below  a  value  of  0.15  which  Indicates  that  the  even 


w' 


correlation  functions  have  greater  variation  In  the  number  of 
sub-peaks  In  this  range. 

Length  1023  Gold .  Figures  A-50  and  A-51  show  the  output  of 
the  Mass-correlate  algorithm  for  the  8  un-optlmlzed  sequences  In 
the  Gold  subset  of  length  1023.  There  Is  little  deviation  across 
the  subset  for  both  the  even  and  odd  threshold  functions. 

Comparison  of  Gold  Code  Sets ♦  Figures  A-52  and  A-53  show 
the  maximum  even  and  odd  threshold  functions  obtained  from  the 


Mass-correlate  algorithm  for  all  of  the  code  sets  defined  In 
Table  2.  As  the  length  of  the  code  Increases,  the  performance 
Increases  significantly  for  both  the  even  and  odd  threshold 
functions . 

Kasaml  Codes 

The  tap  assignments  used  for  the  generation  of  Kasaml  code 

sets  are  given  by  the  polynomials  In  Table  3.  Each  polynomial 

N/  2 

yields  a  set  of  2  code  sequences,  of  which  the  number  used 

for  analysis  Is  Indicated  In  the  table.  Three  different  lengths 
of  Kasaml  code  sets  are  analyzed;  63,  255,  and  1023. 

Table  3 

Code  Generation  Polynomials  for  Kasaml  Code  Sets 


I 

k 

0 

o 

o 


<>• 


Number 

Generated 

8 

10 

8 


Length  Polynomial 

63  [1,6] 

255  [2,3,4,81 

1023  [3,10] 
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and  odd  functions  within  the  set  significantly  vary  in  terms  of 
sub-peaks  below  a  threshold  of  0.28. 


Length  127  Gold. 

Figures  A-40 

and  A- 41 

show  the  even 

and 

cross-correlation 

functions  for 

codes  3 

and  4,  length 

127 

Gold  code  set.  The  even  correlation  function  is  three  valued 
[15,  -1,  -17]  at  integer  values  of  tau.  The  odd  correlation 
function  is  Irregular  with  max  peaks  at  +  28.  The  threshold 
algorithm  plot  in  Figure  A-42,  of  these  two  functions  shows  the 
even  function  having  fewer  peaks  above  a  threshold  of  0.14  and 
more  sub-peaks  below  this  threshold. 

The  off-peak  odd  auto-correlatlon  of  code  3  and  the  odd 
cross-correlation  of  codes  3  and  4  were  optimized  In  terms  of 
phase  and  maximum  peaks  and  and  yielded  the  correlation  functions 
In  Figures  A-43  and  A-44.  The  worst  and  best  phase  of  the  odd 
auto-correlatlon  function  was  thresholded  to  provide  the  plot  In 
Figure  A-45.  It  Is  apparent  that  there  Is  a  significant  differ¬ 
ence  In  the  number  of  peaks  and  sub-peaks  In  threshold  values 
above  0.08.  Figure  A-46  shows  similar  results  for  the  worst  and 
best  phase  of  the  odd  cross-correlation  function  above  a  thresh¬ 
old  of  0.75.  When  the  minimum  odd  and  even  cross-correlation 
functions  are  compared  In  Figure  A-47  the  even  function  appears 
to  provide  more  interference  below  a  threshold  of  0.14. 

Figures  A-48  and  A-49  show  the  output  of  the  Mass-correla¬ 
tion  algorithm  for  the  10  un-optlmized  Gold  codes  In  the  subset 
of  length  127.  The  even  functions  show  no  variation  above  a  0.14 
threshold  while  the  odd  functions  begin  varying  at  a  threshold 
value  of  0.22.  The  even  threshold  varies  much  more  than  the  odd 


Table  2 


Code  Generation  Polynomials  for  Gold  Code  Sets 


Number 

Length  Polynomial  Generated 


63 

[1,6] 

10 

127 

[3,71 

10 

1023 

[3,10] 

8 

plot  of  this  cross  correlation  for  both  even  and  odd  functions. 
Both  functions  perform  similarly  between  threshold  values  of  0.30 
and  0.25.  The  odd  function  performs  slightly  better  than  the 
even  between  'threshold  values  of  0.15  and  0.25.  This  can  be 
attributed  to  fewer  sub-peaks  In  the  odd  function. 

The  off-peak  odd  auto-correlatlon  of  code  3  and  the  odd 
cross-correlation  of  codes  3  and  4  were  optimized  by  varying  the 
phase  and  obtaining  the  minimum  peak  functions  depicted  In  Fig¬ 
ures  A-33  and  A-34.  The  worst  and  best  phase  of  the  odd  auto  and 
cross-correlation  functions  Is  shown  in  Figures  A-35  and  A-36. 
The  difference  between  the  worst  and  best  phase  Is  significant  In 
both  cases,  and  Indicates  that  optimization  of  phase  relation¬ 
ships  will  realize  significant  gains.  Figure  A-37  depicts  the 
optimized  odd,  and  even  cross-correlation  functions  In  terms  of 
the  threshold  algorithm.  This  figure  Indicates  that  the  even 
function  has  more  peaks  In  the  region  above  0.18. 

Figures  A-38  and  A-39  depict  the  output  of  the  Mass-correla¬ 
tion  algorithm  for  ten  of  the  un-optlmlzed  Gold  codes  In  the  set 
of  codes  of  length  63.  These  figures  Indicate  that  both  the  even 


codes  Is  plotted  in  Figures  A~24  and  A-25.  There  Is  little 
variation  above  a  threshold  value  of  0.15  for  both  even  and  odd 
threshold  plots.  Below  0.15,  the  even  threshold  shows  much  more 
variation  than  the  odd  threshold. 

Length  1023  Maximal .  Figures  A-26  and  A-27  show  the  output 
of  the  Mass-correlate  algorithm  for  all  8  un-optlmlzed  Maximal 
codes  of  length  1023.  There  Is  little  variation  above  a  thresh¬ 
old  of  0.075  for  both  the  even  and  odd  threshold  plots.  Below 
0.075,  the  even  threshold  shows  much  more  variation  again. 

Comparison  of  Maximal  Code  Sets .  Figures  A-28  and  A-29 
show  the  maximum  even  and  odd  thresholds  obtained  from  the  Mass- 
correlate  algorithm  for  all  the  Maximal  code  sets  defined  herein. 
As  the  length  of  the  code  Increases,  the  performance  dramatically 
Increases  for  both  the  even  and  odd  threshold  functions.  •  The 
gains  made  from  shorter  to  longer  code  sets  appear  to  gradually 
decrease  In  magnitude  as  the  order  Increases. 

Gold  Codes 

The  tap  assignments  used  for  the  generation  of  Gold  code 
sets  are  given  by  the  polynomials  In  Table  2.  Each  pair  of 
polynomials  yields  a  set  of  2  +  2  code  sequences,  of  which  the 
number  used  for  analysis  Is  indicated  In  the  table.  Three  dif¬ 
ferent  lengths  of  Gold  code  sets  are  analyzed;  63,  127,  and  1023. 

Length  63  Gold .  Figures  A-30  and  A-31  depict  the  even  and 
odd  cross-correlation  functions  for  codes  3  and  4.  The  even 
correlation  function  Is  three  valued  [15,  -1,  -17]  at  Integer 
values  of  tau.  The  odd  cross-correlation  function  Is  Irregular, 
and  has  peaks  of  19  and  -18.  Figure  A-32  shows  the  threshold 


threshold  functions  of  each  Is  plotted  in  Figure  A-13. 
case  the  odd  cross-correlation  has  more  Impact. 


For  this 


The  optimization  of  the  odd  auto  and  cross-correlation  of 
code  1  and  2  produce  the  correlation  functions  shown  In  Figure 
A-14  and  15.  Figure  A-16  shows  the  worst  and  best  phase  of  the 
auto-correlatlon  function  of  code  1.  Again,  It  shows  a  signifi¬ 
cant  difference  between  the  worst  and  best  odd  auto-correlatlon. 
Figure  A-17  show  similar  results  for  the  odd  cross-correlation  of 
codes  1  and  2.  The  even  threshold  function  In  this  figure  as 
well  as  Figure  A-18  Indicates  that  It  has  slightly  less  Impact 
than  either  the  worst  or  best  odd  cross-correlation  function  for 
thresholds  above  0.14.  As  the  threshold  Is  lowered  below  0.14, 
the  even  cross-correlation  has  greater  Impact. 

Figures  A-19  and  A-20  show  the  output  of  the  Mass-correlate 
algorithm  for  all  10  un-.optlmlzed  codes  of  length  127.  There  is 
little  variation  with  threshold  above  0.15.  The  variation  In¬ 
creases  slightly  below  that  threshold. 

Length  255  Maximal.  The  odd  and  even  cross  correlation 
functions  for  Maximal  code  1  and  2  of  Length  255,  are  shown  in 
Figures  A-21  and  A-22.  The  even  function  Is  5  valued  and  ap¬ 
pears  to  be  more  Irregular  than  previous  even  cross-correlation 
functions.  The  odd  function  Is  Irregular  and  has  more  peaks 
around  40  than  the  even  function.  Figure  A-23  depicts  the 
threshold  functions  of  both  the  even  and  odd  cross-correlation. 
Neither  threshold  plot  has  significant  values  above  0.13,  but  the 
even  function  shows  more  peaks  below  this  threshold.  The  Mass- 
correlation  algorithm  for  the  entire  code  set  of  8  un-optlmlzed 
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Important.  This  portion  of  the  function  characterizes  the  peak 
as  well  as  all  of  the  significant  sub-peaks.  For  this  cross- 
correlation,  the  even  and  odd  functions  are  similar  with  the  even 
function  .providing  slightly  more  Interference  on  the  1th  correl¬ 
ators  received  signal. 

The  off-peak  odd  auto-correlatlon  of  code  1  and  the  odd 
cross-correlation  of  codes  1  and  2  were  optimized  In  terms  of 
phase  and  maximum  peaks  and  are  depicted  In  Figures  A-4  and  A-5. 
The  worst  and  best  phase  of  the  odd  auto-correlatlon  threshold 
functions  Is  depicted  In  Figure  A-6.  All  of  the  other  phase 
shifts  of  code  1  have  threshold  plots  between  these  two.  The 
difference  between  the  worst  and  best  phase  is  significant  In 
terms  of  the  odd  auto-correlatlon.  Figure  A-7  shows  the  minimum 
and  maximum  odd  cross-correlation.  The  minimum  threshold  func¬ 
tion  reaches  zero  quicker  than  the  maximum  function  thereby 
providing  slight  Improvement.  When  the  minimum  odd  cross-correl¬ 
ation  Is  compared  to  the  even  cross-correlation  function  In 
Figure  A-8,  the  even  function  appears  to  provide  more  Interfer¬ 


ence. 


Figures  A-9  and  A-10  show  the  output  of  the  Mass-correlate 
algorithm  for  all  six  of  the  un-optlmlzed  Maximal  codes  of  length 
63.  These  figures  show  little  variation  In  the  code  set  In  terms 
of  peaks  and  sub-peaks  over  all  the  cross-correlations,  but  show 
significant  variations  for  lower  thresholds  (below  0.3). 

Length  127  Maximal .  Figures  A-11  and  A-12  show  the  odd  and 
even  cross-correlation  functions  for  code  1  and  2,  length  127 
Maximal  code  set.  Again  the  even  function  Is  three  valued  at 
Integer  value  of  tau  and  the  odd  function  Is  Irregular.  The 
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Table  1 

Code  Generation  Polynomials  for  Maximal  Length  Code  Sets 


Length 


Name 


Polynomial 

[1.6] 

[1.2. 5. 6] 

[2.3,5,61 

[5.6] 

[1,4,5,61 

[1.2,3,61 

[3.71 

[1.2,3,71 

[1,2,4,5,6,71 

[2.3,4.71 

[1,2,3,4,5,71 

[2,4,6,71 

[1.7] 

[1,3.6,71 

[2.5,6,71 

[4.7] 

[2.3,4,81 

[3,5,6,81 

[1,2,5,6,7,81 

[1,3,5,81 

[2.5,6,81 

[1,5,6,81 

[1,2,3,4,6,81 

[1,6,7,81 

[3,101 

[2,3,8,101 

[3.4.5.6.7.8.9.10] 

[1.2.3.5.6.10] 
[2,3,6,8,9,101 
[1,3,4,5,6,7,8,101 
[7,101 

[2.7.8.10] 


threshold  function  Is  due  to  the  threshold  passing  through  the 
value  +  1,  and  Is  common  for  most  even  cross-correlation  func¬ 
tions  of  linear  codes.  The  odd  cross-correlation  threshold  func¬ 
tion  has  a  more  general  slope  and  this  can  be  attributed  to  the 
Irregularity  of  the  odd  function.  In  terms  of  the  system  model 
In  Chapter  II,  the  lower  part  of  the  threshold  functions  are  most 
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IV.  Code  Performance  Analysis 


In  this  chapter,  the  performance  of  three  classes  of  codes, 
as  measured  by  the  algorithms  presented  in  Chapter  III,  are 
analyzed.  The  three  classes  analyzed  are  Maximal  Length  Codes, 
Gold  Codes  and  Kasami  Codes.  The  performance  of  individual 
codes  within  a  code  set  are  analyzed  first,  followed  by  a  mea¬ 
sure  of  the  code  set  Itself.  This  analysis  is  repeated  for 
varying  length  code  sets  for  each  code  class.  Then  the  different 
length  code  sets  are  compared  within  each  class.  Finally,  the 
three  classes  are  compared  for  varying  lengths  of  code  sets. 

Appendix  A  contains  all  the  figures  referenced  in  this 
chapter.  This  appendix  contains  actual  even  and  odd  correlation 
functions  for  code  sets,  as  well  as  plots  from  the  thresholding 
algorithm  (threshold  plots). 

Maximal  Length  Codes 

The  tap  assignments  used  for  the  generation  of  Maximal 
length  code  sets  are  given  by  the  polynomials  in  Table  1.  Four 
different  lengths  of  Maximal  length  codes  are  analyzed;  63,  127, 
255  and  1023. 

Length  63  Maximal.  Figures  A-1  and  A-2  depict  the  even  and 
odd  cross-correlation  functions  for  codes  1  and  2.  It  should  be 
noted  that  the  even  correlation  function  Is  three  valued  [15,  -1 , 
-17]  at  Integer  values  of  tau.  The  odd  cross-correlation  func¬ 
tion  Is  more  erratic  however,  and  has  maximum  peaks  at  17  and 
19.  Figure  A-3  shows  the  threshold  plot  of  this  cross-correla¬ 
tion  for  both  functions.  The  apparent  discontinuity  of  the  even 
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Figure  111-4.  Structure  Chart  for  Mass  Correlate  Algorithm 


with  the  exception  of  the  modules  De termlneName  and  Update 
DataLlst.  The  program  Is  much  more  limited  in  user  options  than 
the  PartlalCorrelate  program.  This  program  is  designed  to  calcu¬ 
late  both  the  Even  and  Odd  correlation  functions,  and  performs 
both  the  positive  and  negative  thresholding  process  on  each  to 
determine  the  amount  of  the  correlation  function  exceeding  these 
thresholds.  However,  this  program  will  perform  this  process  on 
an  entire  code  set  rather  than  just  two  code  sequences.  All  of 
modules  in  this  program  that  have  the  same  names  as  those  in  the 
PartlalCorrelate  program  perform  the  same  function.  The  Module 
DetermineHame  is  called  to  convert  the  entered  filename  extens¬ 
ions  to  Integer  values,  In  order  to  set  up  looping  controls 
within  the  program,  which  perform  correlations  of  pairs  of  code 
sequences.  Instead  of  writing  the  threshold  and  amount  exceeding 
threshold  directly  to  an  output  file,  this  program  will  enter 
this  information  into  a  datalist  in  memory  using  module  Update- 
DataList.  Each  record  in  the  list  contains  fields  for  the 
threshold  value,  the  minimum,  maximum  and  average  value  for  both 
even  and  odd  functions.  As  each  pair  of  code  sequences  Is  pro¬ 
cessed,  the  fields  of  the  data  list  are  updated  accordingly. 
When  all  code  sequences  in  the  user  defined  code  set  have  been 
cross-correlated  and  thresholded,  the  program  will  write  the  data 
list  to  a  user  designated  file. 


Phase 


cross-correlation  functions,  the  Kasami  set  appears  to  show  bet¬ 
ter  performance  in  the  range  above  0.1. 

Figures  A-79  and  A-80  show  the  maximum  even  and  odd  thresh¬ 
olds  obtained  from  the  Mass-correlation  algorithm  for  all  three 
un-optimlzed  classes  of  codes  of  length  1023.  In  the  threshold 
range  0.03  to  0.10  the  Kasami  set  is  clearly  the  best  performer 
in  terms  of  the  number  of  peaks  in  the  even  threshold  functions. 
In  the  same  range  the  Maximal  set  slightly  out-performs  the  Gold 
Set.  In  the  odd  threshold  functions  of  the  same  range,  the  three 
sets  are  very  similar.  The  Kasami  set  slightly  out-performs  the 
other  two  followed  by  the  Maximal  and  Gold  sets. 
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times  listed  In  Table  5.  Based  on  this  data,  the  following 
equation  can  be  used  to  approximate  the  required  runtime  as  a 
function  of  the  length  of  the  code  sequences  to  be  correlated: 

Runtime  “  (1.326  x  10  +  (1.2  x  10  ^)L  (secs) 

Typical  runtimes  for  the  correlation  process  are  listed  In  Table 

5. 


Table  5 

Correlation  Software  Performance  Times  (Secs) 


Code  Length 

63 

127 

255 

511 

1023 

2047 


Runtimes 

2 

3.5 

14 

48 

148 

570 


Phase-Correlate  Program  Performance 

The  Phase-Correlate  program  must  cycle  through  L  correlation 
processes  to  find  the  best  phase  performance  of  each  pair  of 
codes.  For  this  reason,  the  runtime  required  for  each  pair  of 
sequences  baloons  rapidly.  The  runtime  function  can  be 
approximated  as  follows: 

Runtime  ■  (1.326  x  10  ^)L^  +  (1.2  x  10  ^)L^  +  OH  (secs) 

OH  is  the  overhead  time  associated  with  storing  the  maximum  peaks 
for  each  phase,  and  Is  small  relative  to  the  other  factors  In  the 
equation.  It  Is  obvious  that  attempting  to  optimize  an  entire 
code  set  for  the  best  phase  cross-correlation  relationships  would 
be  futile  for  any  reasonable  length.  For  this  reason,  heuristic 
approaches  must  be  used  to  "optimize”  the  performance  of  each 
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code  set. 


Thresholding  Process 

The  thresholding  algorithm  runtime  varies  with  the  the 
length  of  the  code,  correlation  functions  themselves,  and  the 
number  of  threshold  values  required.  All  runs  used  In  this 

work,  varied  the  threshold  values  from  0  to  the  highest  peak  In 
the  correlation  function  being  thresholded.  Generally  the 
threshold  process  (using  approximately  150  thresholds)  required 
less  time  than  the  correlation  process  for  all  work  In  this 
thesis.  In  order  to  shorten  required  runtimes,  the  threshold 
region  can  be  reduced  to  a  smaller  region  of  Interest.  At 
worst,  the  runtime  of  this  process  Is  a  linear  function  of  length 
of  the  code. 

Mass-correlate  performance 

2 

The  Mass-correlate  program  performs  (N  -  N)/2  correlations, 

where  N  Is  the  number  of  code  sequences  within  the  user  defined 

code  set.  The  runtime  for  this  program  Is  a  function  of  the 

cross-correlations  required,  the  length  of  the  code  sequences, 

and  the  number  of  thresholds  of  interest,  and  Is  In  the  order  of 
2 

L  .  This  software  can  be  Improved  by  limiting  the  thresholds 
used  to  a  desired  region  of  Interest. 
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VI 


Conclusions  and  Recommendations 


Conclusions 

From  the  analysis,  of  Chapter  IV  It  can  be  concluded  that  the 
thresholding  algorithm  can  be  used  to  effectively  evaluate  the 
performance  of  binary  code  sequences  In  the  phase-coded  code- 
division  multiple-access  (CDMA)  model.  The  MassCorrelate  program 
will  perform  the  thresholding  algorithm  on  an  entire  code  set  of 
binary  sequences  and  can  be  used  to  evaluate  the  performance  of 
the  set.  The  computing  time  required  to  run  the  MassCorrelate 
program  can  be  reduced  significantly  by  reducing  the  region  of 
threshold  Interest  as  dictated  by  particular  receiver  specifica¬ 
tions  . 

Recommendations 

The  following  recommendations  are  proposed  for  further  study 
In  this  area: 

1.  Using  algorithm?  for  finding  suboptlmal  sequence  phases, 
as  outlined  by  Gahutu  (25:22-26),  and  evaluating  the  performance 
of  these  and  any  other  optimization  techniques,  using  the 
thresholding  algorithm. 

2.  Optimize  the  software  in  terms  of  memory  requirements 
and/or  runtime  to  Improve  performance  of  the  algorithms. 

3.  In  order  to  reduce  the  runtime  of  the  correlation 
algorithm  for  long  code  sequences,  It  seems  feasible  to  Implement 
this  algorithm  using  discrete  Fourier  techniques  on  an  array 
processor. 

4.  Incorporating  additional  code  generation  capabilities 


In  the  GenerateCodes  program,  to  Include  other  linear  and  non¬ 
linear  generators,  so  that  these  code  sets  can  be  evaluated  by 
the  thresholding  algorithm  and  compared  with  the  Maximal  length. 
Gold  and  Kasaml  code  sets. 
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Appendix  A 

Performance  Plots 

This  appendix  contains  all  the  peformance  plots  of  the  data 
obtained  from  the  PartlalCorrelate  program  and  the  MassCorrelate 
program.  Its  order  Is  parrallel  to  that  of  the  order  of 
presentation  In  Chapter  IV.  All  plots  were  obtained  using  the  S 
plotting  package  available  on  the  VAX  11-780  SSC,  and  the  HP  7220 
plotter.  Each  plot  that  Is  obtained  from  the  threshold  algorithm 
Is  normalized  In  terms  of  length  of  the  code  sequence.  All  plots 
of  actual  correlation  functions  are  unnormallzed . 
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Fig.  A-25.  Threshold  Plot  of  Minimum,  Maximum  and  Average 
Odd  Correlation  Functions  from  Mass-Correlation 
Algorithm  for  Maximal  Codes  1-8,  Length  255 
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Fig.  A-26.  Threshold  Plot  of  Minimum,  Maximum  and  Average 
Even  Correlation  Functions  from  Mass-Correlation 
Algorithm  for  Maximal  Codes  1-8,  Length  1023 
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Fig.  A-19.  Threshold  Plot  of  Mlalmua,  Maximum  and  Average 
Even  Correlation  Functions  from  Mass-Correlation 
Algorithm  for  Maximal  Codes  1-10,  Length  127 
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A-20.  Threshold  Plot  of  Minimum,  Maximum  and  Average 
Odd  Correlation  Functions  from  Mass-Correlation 
Algorithm  for  Maximal  Codes  1-10,  Length  127 
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Fig.  A-17.  Threshold  Plot  of  Minimum  and  Maximum  Odd  Cross 
Correlation  and  Even  Cro88'"Correlatlon  of  Maximal 
Codea  1  &  2,  length  127 
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Pig.  A-18.  Threshold  Plot  of  Minimum  Odd  Auto-Correlatlon 
Minimum  Odd  Cross-Correlation,  and  Even  Cross-Cor¬ 
relation  of  Maximal  Codes  1  &  2,  Length  127 


A-16 


-  •  ' ' 


ami;  •xo4t«idln3  ^hrMahold 


Fig,  A-15.  Optimal  Odd  Cross-Correlation  Function  of  Maximal 

Codes  1  &  2»  Length  127 


Fig.  A-16.  Threshold  Plot  of  Minimum  and  Maximum  Odd  Auto 
Correlation  and  Even  Auto-Correlatlon  of  Maximal 
Codes  I  &  2,  Length  127 
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Fig.  A-3.  Typical  Amount  Exceeding  Threshold  vs  Threshold  of  Even 
and  Odd  Correlation  Functions  Maximal  Length  Code  of  Length  63 


Fig.  A-4.  Optimal  Odd  Auto*~Correlatlon  Function  of  Maximal 

Codes  1  &  2.  Length  63 
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A-78  Threshold  Plot  of  Maximum  Odd  Correlation  Func¬ 
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Threshold  Plot  of  Maximum  Odd  Correlation  Func¬ 
tions  from  Mass-Correlation  Algorithm  for  Maximal 
Gold  and  Kasami  Codes  of  Length  1023 . 
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Fig.  A-27.  Threshold  Plot  of  Minimum,  Maximum  and  Average 
Odd  Correlation  Functions  from  Mass-Correlation 
Algorithm  for  Maximal  Codes  1-8,  Length  1023 
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Fig.  A-28.  Threshold  Plot  of  Maximum  Even  Correlation  Func 
tlons  from  Mass-Correlation  Algorithm  for  Maximal 
Codes  of  Length  63,  127,  255,  1023 
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A-33.  Optimal  Odd  Auto^Correlatlon  Function  of  Gold 
Codes  3  &  4,  Length  63 


12  Z2  32  42  52  52 


Fig.  A-34.  Optimal  Odd  Cross^Correlatlon  Function  of  Gold 

Codes  3  &  4,  Length  63 
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Fig.  A-38.  Threshold  Plot  of  Minimum,  Maximum  and  Average 
Even  Correlation  Functions  from  Mass*-Correlatlon 
Algorithm  for  Gold  Codes  l-IO,  Length  63 
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Fig.  A-39.  Threshold  Plot  of  Minimum,  Maximum  and  Average 
Odd  Correlation  Functions  from  Mass-Correlation 
Algorithm  for  Gold  Codes  1-10,  Length  63 
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Fig.  A-41.  Typical  Odd  Cross  Correlation  Function  of 
Gold  Code  of  Length  127 
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Fig.  A-42.  Typical  Amount  Exceeding  Threshold  vs  Threshold  of 
Even  and  Odd  Correlation  Functions  Gold  Code  of  Length  127 
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Pig.  A-43,  Optimal  Odd  Auto-Correlatlon  Function  of  Gold 

Codes  3  &  4,  Length  127 
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Fig.  A-49.  Threshold  Plot  of  Minimum,  Maximum  and  Average 
Odd  Correlation  Functions  from  Mass-Correlation 
Algorithm  for  Gold  Codes  1-10,  Length  127 
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Fig.  A-50.  Threshold  Plot  of  Minimum,  Maximum  and  Average 
Even  Correlation  Functions  from  Mass-Correlation 
Algorithm  for  Gold  Codes  1-8,  Length  1023 
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Fig.  A-51.  Threshold  Plot  of  Mlalmum,  Maximum  and  Average 
Odd  Correlation  Functions  from  Mass-Correlation 
Algorithm  for  Gold  Codes  1-8,  Length  1023 
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.-57.  Optimal  Odd  Auto-Correlation  Function  of  Kasaml 
Codes  2  &  3,  Length  63 
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Fig.  A-58.  Optimal  Odd  Cross-Correlation  Function  of  Kasaml 

Codes  2  &  3,  Length  63 


Fig.  A-59.  Threshold  Plot  of  Mlaimum  and  Maximum  Odd  Auto* 
Correlation  and  Even  Auto-Correlatlon  of  Kasaml 
Codes  2  &  3,  Length  63 


Fig.  A-60.  Threshold  Plot  of  Minimum  and  Maximum  Odd  Cross- 
Correlation  and  Even  Cross-Correlation  of  Kasaml 
Codes  2  &  3,  Length  63 
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Fig.  A-61.  Threshold  Plot  of  Minimum  Odd  Auto-Correlatlon , 
Minimum  Odd  Cross-Correlation,  and  Even  Cross-Correlation 
of  Kasaml  Codes  2  &  3,  Length  63 


Pig.  A-62.  Threshold  Plot  of  Minimum,  Maximum  and  Average 
Even  Correlation  Functions  from  Mass-Correlation 
Algorithm  for  Kasaml  Codes  1-8,  Length  63 


A-38 


Min  Odd  Cross 
Max  Odd  Cross 
Ave  Odd  Cross 


\ 


0 


0.  1  0.  2  0.  3  0.  4 

t.hr«*hold  a*  f  CL) 


Fig.  A-63.  Threshold  Plot  of  Minimum,  Maximum  and  Average 
Odd  Correlation  Functions  from  Mass-Correlation 
Algorithm  for  Kasaml  Codes  1~8,  Length  63 


Fig.  A-64.  Typical  Even  Cross  Correlation  Function  of 
Kasaml  Code  of  Length  255 
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Pig.  A-69.  Threshold  Plot  of  Minimum,  Maximum  and  Average 
Even  Correlation  Functions  from  Mass~Correlatlon 
Algorithm  for  K'asaml  Codes  1-8,  Length  1023 


Fig.  A-70.  Threshold  Plot  of  Minimum,  Maximum  and  Average 
Odd  Correlation  Functions  from  Mass-Correlation 
Algorithm  for  Kasaml  Codes  1-8,  Length  1023 
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Fig.  A-71.  Threshold  Plot  of  Maximum  Even  Correlation  Func 
tions  from  Mass-Correlation  Algorithm  for  Kasami 
Codes  of  Length  63,  255,  1023 
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Fig.  A-72.  Threshold  Plot  of  Maximum  Odd  Correlation  Func 
tions  from  Mass-Correlation  Algorithm  for  Kasami 
Codes  of  Length  63,  255,  1023 
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Fig.  A-73.  Threshold  Plot  of  Maximum  Even  Correlation  Func¬ 
tions  from  Mass-Correlation  Algorithm  for  Maximal, 

Gold  and  Kasami  Codes  of  Length  63 
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Fig.  A-74.  Threshold  Plot  of  Maximum  Odd  Correlation  Func¬ 
tions  from  Mass-Correlation  Algorithm  for  Maximal, 

Gold  and  Kasami  Codes  of  Length  63 
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begin 

M  0; 

L  :*  1; 

while  (M  <  N)  and  (L<16383)  do  begin 
L  2  *  L; 

M  M  +  1; 

end ; ( *while  loop* ) 
for  !;■  1  to  N  do 
Reglfl];-  'O'; 

Regl[N]:- 
M;-  1; 

while  M  <  L  do  begin 

ModSumTapsdapl  ,N,Regl  ,CK) ; 

Enque(Listl,Backl,CK) ; 

ShlftRegister(Regi,N,CK) ; 

M:-  M  +  1; 
end;  C*while  loop*) 
for  !;■  1  to  N  do 
Regl[I]:-  'O'; 

RegllN]:-  '1'; 

M:-  1; 

while  M  <  L  do  begin 

ModSumXaps(Tap2 ,N,Regl ,CR) ; 

Enque (Li8t2,Back2,CK); 

ShiftRegister(Regi ,N,CK) ; 

M:-  M  +  1;  ' 

end;  (*whlle  loop*) 

wrlteln('How  many  codes  do  you  want  generated?  You  may  choose') 
writelnC 'from  1  to  ',L+1,'  Please  input  number’); 
readln(Limlt) ; 
if  Limit  >■  1  then  begin 
Filename;"  Name  +  ’.001’; 
writeln(Fllename) ; 

WriteListCLlstl,  Filename); 
end;  (*lf*) 

if  Limit  >■  2  then  begin 

Filename;"  Name  +  '.002'; 
wrlteln(Fllename) ; 

WriteLlst(List2, Filename); 
end;(*if*) 

Limit;"  Limit  -  2; 

X;"  1; 

while  (X  <  L)  and  (X  <"  Limit)  do  begin 
str(X+2;3,st); 

if  8t[l]  ■  '  '  then  at[l];"  'O'; 
if  st[2j  "  '  '  then  st[2];"  'O'; 

Filename ; "  Name  +  ' . '  +  st ; 
writeln(Filename) ; 

ModSimList (Li  St l,Llst2, Filename) ; 

Shif tLlst(Listl) ; 

X;"  X+1; 

end  (*whlle  loop*) 
end ; ( *procedure  GenGold* ) 


(*  *) 

(*  GenGold  *) 

(*  *) 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  A. 21  *) 

(*  DESIGN  DATE:  Dec  84  LAST  UPDATE:  4  Jan  85  *) 

(* - *) 

(*  DESCRIPTION:  This  procedure  constructs  and  manipulates  two  linked  *) 

(*  lists  representing  the  outputs  of  two  LFSR's.  The  shifting  and  *) 

(*  modulo-2  addition  of  the  lists  are  required  to  produce  Gold  codes.  *) 

(*  The  user  has  the  option  of  producing  1  to  the  maximum  number  of  *) 

(*  unique  gold  codes  the  registers  will  produce  and  stores  the  codes  *) 

(*  In  text  files  with  the  same  name  and  different  extensions.  *) 

(  * - *  ) 

(*  CALLING  ARGUMENTS:  Name — string  of  characters  representing  the  file  *) 

(*  the  user  desires  to  store  output  *) 

(*  Tapi — Tap  assignments  for  LFSR  *) 

(*  N — Length  of  LFSR  used  to  generate  code  *) 

(* - *) 

(*  INTERNAL  VARIABLES:  I,M,X — counters  used  for  looping  *) 

(*  L — Length  of  LFSR  sequences  *) 

(*  Limit — Max  number  of  Gold  Codes  to  be  produced  *) 

(*  CK — output  of  LFSR  after  addition  of  taps  assigned*) 

(*  Coded ta — Text  file  variable  for  storage  of  output  *) 

(*  st — string  of  character  representing  filename  ex-  *) 

(*  tensions  for  storage  of  output  *) 

(*  Regl — LFSR  array  contents  *) 

(*  Backl,Back2 — pointers  to  end  of  lists  *) 

(*  Filename — string  of  12  characters  w/extenslon  ♦) 

(* - *) 

(♦  CONSTANTS:  None  *) 

(  * - : - *  ) 

(*  FILES  USED:  Codedta — storage  of  outputted  code  *) 

(* - *) 

(*  MODULES  CALLED:  ModSumTaps  *) 

(*  Enque  *) 

(*  ShiftRegister  *) 

(*  ModSumLlst  *) 

(*  ShiftLlst  *) 

(*  WrlteLlst  *) 

(♦ - *) 

procedure  GenGold (Ncune: Word;  N:lnteger;  Tapl,Tap2:RegStructure; 

var  Llstl,  Llst2:dataptr) ; 


var  Llmit,I,M,L,X:lnteger; 
CK: Binary; 

Codedata:  text; 

St:  strlng[3]; 

Regl:  RegStructure ; 
Backl , Back2 : dataptr ; 
Filename : NameArray ; 
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(*  *) 

(*  Decimate  *) 

(*  *) 

(*«***********************************************************************) 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  A.311  *) 

(*  DESIGN  DATE;  Feb  85  LAST  UPDATE;  14  Mar  85  *) 

(* - *) 

(*  DESCRIPTION;  This  procedure  decimates  a  linked  list  representing  a  *) 

(*  maximal  sequence  of  length  Length,  at  decimates  by  a  factor  of  *) 

(*  N+2  to  produce  a  sequence  of  length  L.  It  then  uses  the  second  *) 

(*  sequence  to  produce  a  third  sequence  of  period  L,  length  Length.  *) 

(* - *) 

(*  CALLING  ARGUMENTS;  Llstl,Llst2,Llst3 — pointers  to  beginning  of  lists  *) 

(*  Backl,Back2,Back3 — pointers  to  ends  of  lists  *) 

(*  L — Length  of  sequence  two  *) 

(*  Length — length  of  sequence  three  *) 

(* - *) 

(*  INTEBINAL  VARIABLES;  P,ptr — pointer  used  to  traverse  the  list  *) 

(*  Count ,R — counters  used  to  do  process  *) 

(* - *) 

(*  CONSTANTS;  None  *) 

(* - *) 

(*  PILES  USED;  None  *) 

(* - ; - - *) 

^  (*  MODULES  CALLED;  Enque  *) 

(* - *) 


procedure  DeclmateCvar  Llstl,Backl,Llst2,Back2,Llst3,Back3;  dataptr; 

L, Length ; Integer) ; 

var  P,Ptr; dataptr; 

Count ,K; Integer ; 

begin 

P;-  Listl; 

Enque(Llst2  ,Back2  jP'^  .Info) ; 

Count ;■!; 
repeat 

for  K;*  1  to  (L+2)  do 
P;“  P'‘.ncrt; 

Enque(Llst2  ,Back2  ,P'' .  Info) ; 

Count;*  Count  +  1; 
until  Count  ■  L+1; 

Count;*  0; 

whlleCCount  <  Length)  do  begin 
P;*  List2; 
repeat 

Enque(Llst3,Back3,P'' .Info) ; 

Count;*  Count  +  1; 

P;*  P^.next; 
until  (P'‘.next  *  nil); 
end(*whlle  Count*) 
end;  (*procedure  Decimate*) 
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(^Ic************************************************************************'^ 

(*  *) 

(*  GetRidOf  *) 

(*  *) 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  A.  03  *) 

(*  DESIGN  DATE;  Jan  85  LAST  UPDATE;  14  Jan  85  *) 

(* - *) 

(*  DESCRIPTION:  This  procedure  releases  a  linked  list  from  memory  alio-  *) 

(*  cation  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  List — pointer  to  beginning  of  linked  list  *) 

(  * - *  ) 

(*  INTERNAL  VARIABLES:  P,ptr — pointer  used  to  traverse  the  list  *) 

(* - *) 

(*  CONSTANTS;  None  *) 

(* - *) 

(*  PILES  USED:  None  *) 

(* - *) 

(*  MODULES  CALLED:  None  *) 

(* - *) 

procedure  GetRidOf (var  List:dataptr); 

var  P,ptr:dataptr; • 

begin 

P;-  List; 
ptr:"P; 

while  P''.next  O  nil  do  begin 
ptr:*  P''.nextj 
dlsposeCP) ; 

P;»ptr; 

end  (*whlle  loop*); 
dlsposeCP); 

List:*  nil; 

end ; ( *procedure  GetRidOf* ) 


(******«******************************************«**************★********) 


(*  *) 

(*  WrlteList  *) 

(*  *) 

(4t***********************4k****************4r*4r*****4r4r**************4t**4t^/***) 

(*  DESIGNER:  Richard  C.  Gender  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  A.07  *) 

(*  DESIGN  DATE:  Dec  84  LAST  UPDATE:  10  Dec  84  *) 

(* - *) 

(*  DESCRIPTION:  This  procedure  traverses  a  linked  list  of  records  and  *) 

(*  outputs  the  Information  In  the  records  In  file  Codedta.  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  List— pointer  to  beginning  of  Unked  list  *) 

(*  Filename — string  of  char  Identifying  file  to  be  *) 

(*  written  to  *) 

(* - *) 

(*  INTERNAL  VARIABLES:  P — pointer  used  to  traverse  the  list  *) 

(* - *) 

(*  CONSTANTS:  None  *) 

(* - *) 

(*  FILES  USED:  Writes  to  Codedta  *) 

(• - ^ - *) 

(*  MODULES  CALLED:  None  *) 

(* - *) 

procedure  WrlteLlst(Llst:dataptr:Fllename:NameArray) ; 


var  P:dataptr; 

Codedta: text; 

begin 

asslgnCCodedta, Filename) ; 
rewrite (Codedta) ; 

P:"Ll8t; 

while  P  O  nil  do  begin 
write  (  Codedta ,  P'" .  Inf  o  ) ; 
P:“  P''.next; 
end;  (*whlle  loop*) 
close(Codedta) ; 
end; 


I  TS 


t 


(*******************'*******'******************************4^****************) 
(*  *) 

(*  ModSumLlst  *) 
(*  *) 
(*************************************************************************) 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 
(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  A.  06  *) 
(♦DESIGN  DATE;  Dec  84  LAST  UPDATE;  4  Jan  1985  *) 
* _ 


( 


) 


(*  DESCRIPTION;  This  procedure  traverses  two  Linked  Lists  of  binary  data  *) 
(*  and  performs  a  bit  by  bit  moduIo-2  sumatlon  of  the  lists  and  writes  ♦) 
(*  the  result  Into  the  text  file,  Codedta.  *) 


( 


) 


(*  CALLING  ARGUMENTS:  Llstl,Llst2 — pointer  values  Indicating  the  first  *) 
(*  record  of  each  linked  list  or  records  *) 
(*  Filename — string  of  characters  Indicating  name  of  *) 
(*  file  output  should  be  stored  In  *) 


0 


(*  INTEEINAL  VARIABLES:  Pl,P2 — pointers  used  for  traversing  the  lists  *) 
(*  bit — used  for  result  of  modulo-2  addition  of  two  *) 
(*  records  *) 


(*  CONSTANTS;  None 
(* - 


0 


(*  FILES  USED:  writes  to  Codedta 

*  ■ 


*) 

-*) 

*) 


( 


(*  MODULES  CALLED;  None 

(* 


0 


procedure  ModSufflLlst(Llstl:dataptr;Llst2:dataptr;Fllename:NaffleArray) ; 


*) 

*) 


var  Pl,P2  ;  dataptr; 
Bit: binary; 
Codedta:  text; 


begin 

asslgnC Codedta, Filename) ; 
rewrlte(Codedta) ; 

PI;-  Llstl; 

P2;-  Llat2; 

while  Pl'‘.next  O  nil  do  begin 
If  PI ''.Info  -  P2''.lnfo  then 
Bit;-  'O' 

else 

Bit;-  '1'; 
wrlte(Codedta,Blt) ; 

PI;-  PI*. next; 

P2;-  P2*.next;  ' 
end(*whlle  loop*); 

If  PI''. Info  -  P2''.lnfo  then 
Bit;-  'O' 
else 

Bit;-  '1'; 
wrlte(Codedta,Blt) ; 
close(Codedta); 
end  (*procedure  ModSumLlst*); 
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(*  *) 

(*  ShlftLiat  *) 

(*  *) 

(*************************************************************************) 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  A. 08  *) 

(*  DESIGN  DATE:  Dec  84  LAST  UPDATE:  12  Dec  84  *) 


(*  DESCRIPTION:  This  procedure  takes  a  linked  list  of  records  and  places  *) 
(*  the  first  record  at  the  end  of  the  list.  *) 


(*  CALLING  ARGUMENTS:  List — pointer  to  the  first  record  In  the  linked 


(*  INTEEINAL  VARIABLES:  First, P — pointers  used  for  traversing  the  list 
(*  and  assigning  a  new  beginning  of  the  list 


(*  CONSTANTS:  None 


C*  FILES  USED:  None 


(*  MODULES  CALLED:  None 


procedure  ShlftLlst(var  Llst:dataptr) 
var  Fir s t , P : dataptr ; 


P:“  List} 

while  P^.nezt  O  nil  do 
P;"  P'^.next; 

First:"  First''. next; 
P'*.next:"  List; 

List''. next:"  nil; 

List:"  First; 

end;  (*procedure  ShlftLlst*) 


(^mcitit*********************************************************************) 

(*  *) 

(*  E  n  q  u  e  *) 

(*  *) 

(A************************************************************************) 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  A. 09  *) 

(*  DESIGN  DATE;  Dec  84  LAST  UPDATE;  12  Dec  1984  *) 


(*  DESCRIPTION;  This  procedure  takes  an  Inputted  value  (Bit)  and  places  *) 
(.*  It  In  a  record,  and  then  places  the  record  at  the  bottom  of  the  list  *) 
i*  Identified  by  the  calling  argument  List.  It  updates  the  calling  mod-*  *) 
(*  ule  of  the  pointer  values  for  the  front  and  back  of  the  list.  *) 


(*  CALLING  ARGUMENTS;  List — pointer  Indicating  the  beginning  of  the  list  *) 
(*  Back — pointer  Indicating  the  end  of  the  list  *) 
(*  Bit — the  data  value  to  be  placed  In  the  linked  list*) 


(*  INTERNAL  VARIABLES;  P — temporary  pointer  used  for  traversing  the  list  *) 


(*  CONSTANTS;  None 


(*  FILES  USED;  None 


(*  MODULES  CALLED;  None 


procedure  EnqueCvar  Llst;dataptr;var  Back:dataptr;  Blt;blnary); 


var  P;  dataptr; 


begin 

if  List  ■  nil  then  begin 
New(P); 

List;*  P; 

P''.lnfo;-  Bit; 

Back;*  P; 

Back''. next;*  nil; 
end 

else  begin 
New(PJ ; 

P'.info;-  Bit; 

Back'. next:*  P; 

Back;*  P; 

Back'. Next;*  nil; 
end  (*lf  then  else*) 
end;  (*procedure  Enque*) 


(»  *) 

(*  ShiftRegister  *) 

(*  *) 

(*4t***************'****************«*ii*«*»*.*«*ii*ir*ir'*x*«*K*ir*«*«*»***it*****) 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  A. 04  *) 

(*  DESIGN  date;  Nov  84  LAST  UPDATE:  15  Nov  1984  *) 


(*  DESCRIPTION:  This  procedure  takes  the  contents  of  an  LFSR  and  shifts  *) 
(*  It  to  the  right  one  position.  It  Loads  the  first  position  with  the  *) 
(*  ved.ue  CK  (input  paraaeter).  *) 


(*  CALLING  ARGUMENTS:  C — ^Array  representing  the  LFSR  values 
(*  N — Length  or  the  LFSR 

(*  CK — value  to  be  input  to  first  position  of  LFSR 

(*  after  shifting  the  register  one  position 


(*  INTERNAL  VARIABLES:  I— used  as  a  counter 


(*  CONSTANTS;  None 


(*  FILES  USED;  None 


(*  MODULES  CALLED;  None 


procedure  ShlftReglsterCvar  C:RegStructure;  N: Integer;  CK: binary) 


var  I: Integer; 


begin 
I.—  N; 
repeat 

C[I];-  C(I-1]; 

I;-  I  -  1; 
until  I  <  2; 

CllJ.—  CK; 

end  (*procedure  ShiftRegister*); 


(*  *) 

(*  ModSumTaps  *) 
(*  *) 
(********************************************************************«****) 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 
(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  A. 05  *) 
(*  DESIGN  DATE:  Nov  84  LAST  UPDATE;  15  Nov  84  *) 


(*  DESCRIPTION:  This  procedure  takes  an  LFSR  and  Its  designated  taps  and  *) 
(*  performs  a  modulo  2  addition  of  the  contents  of  the  register  as  des-*) 
(*  Ignated  by  the  taps  and  produces  the  result — CK  *) 


(*  CALLING  ARGUMENTS:  T — Tap  assignments 
(*  N—Length  of  LFSR 
(*  R — ^Reglater  array 

(*  CK — output  character  of  the  operation 


(*  INTERltAL  VARIABLES:  I — used  as  a  counter 


(*  CONSTANTS;  None 


k*  rlLES  USED:  None 


(*  Modules  CALLED:  None 


procedure  ModSumTaps (T:RegStructure;  N: Integer; 

R:RegStructure;var  CK: binary) 


var  I  :  Integer; 


begin 

CK  'O'; 

I;-0; 

repeat 

I:-  I  +  1; 

If  T[I]  ■  '1'  then  begin 
If  CK  -  R[I]  then 
CK  :-  'O' 
else  CK  ;■  '1' 
end;  (*lf*) 
until  I  "  N; 

end; (*procedure  ModSumTaps*) 


(**********4r***ik*4t****4r  ;***************************************************  ^ 

(*  *) 

(*  Loadinput  *) 

(*  *) 

(***********4r***4r*****4t**4r*4tiir**1^******4t*****4t**il!4;**4t***4r*4r******4r******4t**) 

(*  DESIGNER;  Richard  C.  Gonder  .  MASTER  ThESIb  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  ’  MODULE  NUMBER  A. 01  *) 

(*  DESIGN  DATE;  Nov  84  LAST  UPDATE;  15  Nov  84  *) 

(* - *) 

(*  DESCRIPTION:  This  procedure  loads  the  tap  assignments  for  an  LFSR  and  *) 

(*  returns  the  loading  to  the  calling  module.  *) 

(*  *) 

(* - *) 

(*  CALLING  ARGUMENTS;  N—Length  of  the  register  *; 

i*  T — array  of  taps  that  are  assigned  in  module  *) 

(*  Identifier — String  of  characters  used  to  communl-  *) 

(*  cate  with  the  user  *) 

(* - *) 

(*  INTERNAL  VARIABLES:  I—Countlng  variable  *) 

(* - *) 

(*  CONSTANTS;  None  *) 

(* - : - *) 

K*  ifiLES  USED;  None  *) 

(* - *) 

(*  MODULES  CALLED:  None  *) 

(* - * - *) 

procedure  Loadinput (N: Integer;  var  T;RegStructure;  Identifier: StringSO) 

var  I:  Integer; 

begin 

I;"l; 

while  I  <"'N  do  begin 
T[Il;-  'O’; 

I;-I+l; 

end;(*wfalle  loop*) 

wrltelnC'Now  Select  the  positions  for  the  Identifier, 

'  . . .when  complete ' ) ; 
writeln( ' Input  value  greater  than  15'); 
readln(l); 

while  I  <  16  do  begin 
T[I]  ;-  '!'; 
readln(I) 

end  (*while  loop*); 
end;  (*procedure  Loadinput*) 
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(*************************************************************************) 


(*  *) 

(*  Test  *) 

(*  *) 

(********1k**«*****4t***************llt**************************4r************^ 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  A.02  *) 

(*  DESIGN  DATE:  Nov  84  LAST  UPDATE:  15  Nov  84  *) 

(* - *) 

(*  DESCRIPTION:  This  procedure  Tests  to  Insure  that  the  fllenaffle  Is  a  *) 

(*  valid  one  for  an  8-bit  CPM  machine.  If  correct  it  will  return  a  *) 

(*  flag  ■  true.  With  minor  errors  it  will  also  correct  them.  If  the  *) 

(*  error  is  not  recoverable  it  will  return  a  flag  ■  false.  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  Name-  8-letter  string  *) 

(*  flag-  Boolean  indicating  validity  of  Name  *) 

(* - *) 

(*  INTERNAL  VARIABLES;  I— counting  variable  *) 

(*  Len — Length  of  the  string  Name  *) 

(*  ch — used  to  change  the  case  of  Name  *) 

(* - *) 

(*  CONSTANTS;  None 

(* - ^ - *) 

(*  PILES  USED;  None  *) 

(* - *) 

(*  MODULES  CALLED;  None  *) 

(♦ - *) 

procedure  Test(var  Name:Word;  var  flag: boolean); 


var  I, Len:  Integer; 
ch:  char; 


begin 

l:-  1; 

Len:**  Length(Name); 

while  (I  <**  Len)  and  flag  do  begin 

if  Name! I]  in  [ 'A' . . ’Z* , 'a' . . 'z* , ’0’ . . *9* , ’ : ' ]  then  begin 
flag;"  true; 

if  Name[I]  in  ['a'..'z']  then  begin 
ch;"  chr(ord(Name[l3)  -32); 

Name[I]:"  ch; 
end;(*if*) 

end  (*front  clause  of  if*) 
else  flag:"  false; 

I;"  I  +  1; 
end;  (*whlle  loop*) 
end;  (*procedure  Test*) 


(^t'kifit'kic'kltieitititititititicicititicitigicifisic'kitltifiticisiticitititicicic’kli'kicicicicicicit'kicicicitisiticititifititifititis’kitifit'^ 

(*  *) 

(*  GenerateCodes  *) 

(*  *) 

(^Ic************************************************************************) 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODDLE  NUMBER  A.O  *) 

(*  DESIGN  DATE;  Nov  1984  LAST  UPDATE;  4  Mar  1985  *) 


(*  DESCRIPTION:  This  program  will  generate  codes  using  linear  feedback  *) 
(*  shift  registers  (LFSR)  and  store  the  coues  In  ASCII  text  files  as  *) 
(*  designated  by  the  user  for  further  manipulation.  It  c£ui  be  used  for  *) 
(*  generating  Maximal  Length  Codes,  Gold  Codes  and  Kasaml  codes  provided*) 
(*  the  user  Identifies  the  appropriate  tap  assignments  for  each.  This  *) 
(*  version  Is  limited  to  register  lengths  of  15  because  It  was  written  *) 
t*  tor  an  8-blt  machine.  *) 


(*  CALLING  ARGUMENTS:  Std  Input,  Output 


(*  INTERNAL  VARIABLES; 

(*  ch — character  Input  to  the  program  from  user 


(*  CONSTANTS;  None 


(*  FILES  USED:  User  designated  Filename  for  code  storage 


C*  MODULES  CALLED:  LlnearGenerator 
(*  GoldCodeGenerator 

(*  KasamlGenerator 


program  GenerateCodes (Input , output) ; 


Binary  - 

RegStructure  ■  array  [1..15]  of  Binary 
NameArray  ■  strlngtlZ]; 

Word  ■  strlng[8]; 

StrlngSO  ■  strlng[30]; 
dataptr  ■  '‘data; 
daca  ■  record 

Info:  Binary; 
next:  dataptr; 
end; (*record*) 


var  ch:  char; 


Fig.  A-80.  Threshold  Plot  of  Maximum  Odd  Correlation  Func¬ 
tions  from  Mass-Correlation  Algorithm  for  Maximal, 

Gold  and  Kasaml  Codes  of  Length  1023 
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(*  *) 

(*  GenLinear  *) 

(*  *) 

(****************************«********************************************) 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  A. 11  *) 

(*  DESIGN  DATE;  Nov  84  LAST  UPDATE;  4  Jan  85  *) 


(*  DESCRIPTION;  This  procedure  produces  a  sequence  of  bits  from  an  LFSR  *) 
(*  and  writes  the  sequence  Into  file  Codedta.  Maximum  length  of  LFSR  *) 
(*  Is  15  positions.  *) 


(*  CALLING  ARGUMENTS;  Filename — string  representing  storage  file 
(*  N— Length  of  LFSR 

(*  Tapi — array  of  tap  assignments 


(*  INTERNAL  VARIABLES;  M — counter  for  looping 
(*  L.Leng — Length  of  output  sequence 

(*  Maxi — ^Number  of  ones  In  output  sequence 

(*  CK — Output  of  generator  at  each  clock  cycle 

(*  Maxi — nimber  of  ones  generated  In  code 


(*  CONSTANTS;  None 


(*  FILES  USED;  Writes  to  Codedta 


(*  MODULES  CALLED:  ModSunilaps 
(*  ShlftReglster 


procedure  GenLinear ( Filename ; Namearray ; N ; Integer ; var  Tapi ;RegStructure) 

var  Maxi, M,L,Leng; Integer; 

Codedta;  text; 

CK;  Binary; 

Regl ;RegStructure ; 

begin 

Leng:“  0; 

M  ;-  0; 

L  1; 

while  (M  <  N)  do  begin 
L  2  *  L; 

M  ;-  M  +  1; 
end;  (*  while  loop*) 
for  M;*  1  to  N  do 
Regl[M];-  'O'; 

RegllN];- 
L  ;-  L  -  1; 

Maxi;*  0; 

M  ;-  0; 

asslgnC  Codedta , Filename ) ; 
rewrite (Codedta) ; 


■  •  )p  ■  w  »  " 


.VW'.'-'V** 


while  (M  <  L)  do  begin 

Mods umTapsC Tapi ,N,Regl ,CK) ; 

Leng  ;■  Leng  +  1; 

Shlf tReglsterCRegl ,N,CK) ; 
wrlte(Codedta,CK) ; 

If  CK  -  '1’  then 
Maxi;*  Mazl  +1; 

M  M  +  1 
end;  (*while  loop*) 
close(Codedta); 

wrltelnC 'Length  of  code  Is  ',Leng,'  #  of  ones  Is  Maxi) 
end( *procedure  GenLlnear* ) ; 


(ir 
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(*  *) 

(*  GenKasaml  *) 

(*  *) 

(********************************************jk****************************) 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  A. 3.1  *) 

C*  DESIGN  DATE:  Nov  84  LAST  UPDATE:  4  Jan  85  *) 

(* - *) 


(*  DESCRIPTION:  This  procedure  produces  a  sequence  of  bits  from  an  LFSR  *) 
(*  and  places  the  sequence  Into  a  linked  list  In  memory.  It  then  dec-  *) 
(*  Imates  this  sequence  to  produce  a  second  sequence.  The  second  se-  *) 
(*  quence  Is  used  to  produce  a  third  sequence  of  the  period  of  the  se-  *) 
(*  cond  that  Is  the  same  length  of  the  first.  Then  the  procedure  modu-*) 
(*  lo-2  sums  the  first  sequence  with  phase  shifts  of  the  third  to  pro-  *) 
(*  duce  the  set  of  Kasaml  Sequences.  The  first  sequence  must  by  a  max-*) 


(*  sequence  for  the  generator  to  produce  Kasaml  sequences.  *) 

(  * - *  ) 

(*  CALLING  ARGUMENTS:  Name —  representing  storage  file  *) 

(*  N— Length  of  LFSR  *) 

(*  Tapi — ^Tap  assignments  for  LFSR  *) 

(*  Llstl,Llst2, lists — pointers  to  fronts  of  lists  *) 

(*  Backl,Back2,Back3 — pointer  to  ends  of  lists  of  se-  *) 

(*  quences  used  to  generate  Kasaml  Sequences  *) 

(* - *) 

(*  INTERNAL  VARIABLES:  Length — Length  of  sequence  one  *) 

(*  L — Decimating  factor  *) 

(*  Limit— maximum  number  of  Kasaml  Sequences  that  *) 

(*  can  be  produced  In  this  coriflguration  *) 

(*  Regl — array  representing  LFSR  producing  seq  1  *) 

(*  Filename — Filename  w/ extension  added  for  output  *) 

(*  sequences  *) 

(*  CK — output  of  LFSR  produces  seq  1  *) 

(*  I,K,M,X — counters  used  In  procedure  *) 

(*  st — string  of  length  3  used  as  filename  extension  *) 

(*  LjLeng — ^Length  of  output  sequence  *) 

(* - *) 

(*  CONSTANTS;  None  *) 

(* - *) 

(*  FILES  USED:  Writes  to  Codedta  *) 

(* - *) 

(*  MODULES  CALLED;  ModSumTaps  *) 

(*  ShlftReglster  *) 

(*  Enque  *) 

(*  Decimate  *) 

(*  Wrltellst  *) 

(*  ShlftLlst  *) 

(*  ModSumLlst  *) 

(♦ - *) 

procedure  GenKasaml (Name: Word;  N:lnteger;  Tapl:RegStructure; 


var  Llstly  Llst2,Llst3,Backi,Back2,Back3:dataptr); 


var  Length, K, Limit, I, M,L,X:lnteger; 

CK; Binary; 

Codedata: file  of  char; 

St:  strlng[3]; 

Regl:  RegStructure; 

Filename : NameArray ; 

begin 

M  :-  0; 

Length  :■  1; 

while  (M  <  N)  and  (Length<16383)  do  begin 
Length  :■  2  *  Length; 

M  M  +  1; 
end;(*\rtille  loop*) 

Length: “Length- 1 ; 

M:“  0; 

L:“  1; 

while  (M<  N/2)  do  begin 
L:“  2  *  L; 

end(*whlle  loop*); 

L:“  L  -1; 
for  !:■  1  to  N  do 
Regl[Il;“  ’O’; 

Regl[N]:“  ’1’; 

M:-  1; 

while  M  <“  Length  do  begin 
ModSumTapsClapl ,N,Regl ,CK) ; 

Enque(Llstl , Backl , CK) ; 

Shlf  tReglsterCRegl  ,N,CK) ; 

M:-  M  +  1; 
end;  (*whlle  loop*) 

Declmate(Llstl, Backl, Llst2,Back2,Llst3,Back3,L, Length) ; 

M:“  1; 

wrltelnC’How  many  codes  do  you  want  generated?  You  may  choose’) 
wrltelnC’from  1  to  ’,L,’  Please  Input  ntimber’); 
readln(Llmlt) ; 

If  Limit  >“  1  then  begin 
Filename:*  Name  +  ’.001’; 
wrltelnCFUename) ; 

WrlteLlstCLlstl,  Filename); 
end;  (*lf*) 

Limit ;■  Limit  -  1; 

X:“  1; 

while  (X  <  L)  and  (X  <“  Limit)  do  begin 
str(X+l:3,st); 

If  at[l]  “  ’  ’  then  st[l]:“  ’O’; 

If  st[2]  “  ’  ’  then  st[2]:“  ’O’; 

Filename:*  Name  +  ’.’  +  st; 
wrltelnCFUename  ) ; 

ModSumLlstCLlstl ,Llst3, Filename) ; 

ShlftLlst(Llstl); 

X:*  X+1; 

end;  (*whlie  loop*) 
end;(*procedure  GenKasaml*) 


(*  ■  *) 
(*  GoldCodeGenerator  *) 

(*  *) 
(*************************************************************************) 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  A. 2  *) 

(*  DESIGN  DATE:  Dec  84  LAST  UPDATE:  4  Jan  85  *) 


(*  DESCRIPTION:  This  procedure  does  most  of  the  program/ user  Interaction  *) 
(*  to  generate  Gold  codes.  Reads  In  Length  of  LFSR,  Tap  assignments  and  *) 
(*  Filename  for  storage.  *) 


(*  CALLING  ARGUMENTS:  None 


(*  INTERNAL  VARIABLES: 

Regl — array  representing  LFSR 

.mj 

*) 

(* 

Tapl,Tap2 — Tap  assignments  for  two  LFSR's  used  to 

*) 

(* 

produce  two  linear  Sequences 

*) 

(* 

N— length  If  LFSR's 

*) 

(* 

Name — user  Inputted  filename  for  storage  of  Gold 

*) 

(* 

sequences . . . w/ o  extension 

*) 

(* 

Llstl,Llst2 — pointers  to  front  of  LFSR  sequences 

*) 

(* 

Backl,Back2 — pointers  to  end  of  LFSR  sequences 

*) 

(* 

Flag — boolean  used  for  control 

*) 

(* 

(* - 

ch — character  used  for  user  I/O 

*) 

(*  CONSTANTS:  None 


(*  FILES  USED:  None 


procedure  GoldCodeGenerator; 


var  Regl,Tapl,Tap2  :  RegStructure ; 

N  :  Integer; 

Name:  Word; 
ch:  char; 

Flag:  boolean; 

Llstl,Llst2,Backl,Back2:  dataptr; 

begin 

repeat 

Llstl:"  nil; 

Llst2:~  nil; 

wrltelnC Input  the  length  of  the  two  registers, 
..2  -  15  please  ’); 
readln(N); 

while  (N  <  2)  or  (N  >  15)  do  begin 

wrltelnC 'range  must  be  2  to  15  please'); 


readln(N) 

end;  (*\diile  loop*) 

LoadInputCN, Tapi, 'Taps  for  register  nvmber  1'); 
LoadInput(N,Tap2, 'Taps  for  register  number  2'); 
wrltelnC 'What  file  do  you  want  to  store  the  sequence 
generated ' ) ; 

wrltelnC 'No  more  than  8  characters  In  filename'); 
repeat 

Flag;"  true; 
readln(Name) ; 

Test(Name,Flag) ; 

If  not  Flag  then  begin 

wrltelnC 'Invalid  Input  Please  try  again'); 
wrltelnC 'No  Special  Characters  on  Inputted  Name'); 
end;  C*lf  statement*) 
until  Flag; 

wrltelnC ' Code  will  be  stored  In  file  ' ,Name) ; 

GenGoldCName ,N,Tapl ,Tap2 ,Llstl ,Llst2) ; 
wrlteln; 

wrltelnC 'Do  you  wish  to  generate  any  other  Gold  codeflles?') 

wrltelnC 'If  so.  Input  Y  If  not.  Input  N'); 

readlnCch); 

GetRldOf CLlatl) ; 

GetRldOf  CUst2) ; 
tmtll  Cch  ■  'N')  or  Cch  *  'n'); 
end;  C*of  procedure  GoldCodeGenerator*) 


(*  *) 

(*  LlnearGenerator  *) 

(*  *) 

(^titititit********************************************************************') 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  A.l  *) 

(*  DESIGN  DATE;  Nov  84  LAST  UPDATE;  4  Dec  84  *) 


(*  DESCRIPTION;  This  procedure  does  most  of  the  user/program  Interaction  *) 
(*  required  to  produce  a  sequence  from  a  single  LFSR.  The  Length  of  *) 
(*  The  LFSR,  tap  assignments  and  Filename  for  storage  are  obtained  *) 


(*  CALLING  ARGUMENTS:  None 


(*  INTEEINAL  VARIABLES;  Tapi — tap  assignment  for  LFSR 


N—length  of  LFSR 

Filename — string  for  output  filename  w/ext 
Name — string  for  output  filename  w/o  ezt 
Flag — boolean  used  for  control 
ch — character  used  for  user  I/O 


(*  CONSTANTS;  None 


(*  PILES  USED:  None 


(*  MODULES  CALLED:  Loadinput 
(*  Test 

(*  Genlilnear 


procedure  LlnearGenerator; 

var  Tapi  :  RegStructure; 
N  :  Integer; 
Filename;  NameArray; 
Name:  Word; 
ch:  char; 

Flag:  boolean; 


begin 

repeat 

wrlteln( ' Input  the  length  of  the  linear  shift  register'); 
wrltelnC 'Length  should  be  between  2-15  please'); 
readln(N) ; 

while  (N  <  2)  or- (N  >  15)  do  begin 

wrltelnC 'range  must  be  between  2  to  15.. please  enter  again'); 
readln(N); 
end(*whlle  loop*); 

LoadInputCN, Tapi, 'Taps  for  the  shlftreglster' ) ; 
wrltelnC 'What  file  do  you  want  to  store  the  sequence 
generated? ' ) ; 

wrltelnC 'No  more  than  8  characters  In  the  filename'); 


repeat 

Flag;*  true; 
readln(Name) ; 

Test  (Neune,  Flag) ; 

If  not  Flag  then  begin 

writelnC ' Invalid  Input. . .Please  try  again' ) ; 
wrltelnC 'No  special  characters  In  filename ' ) ; 
end; (*lf*) 
until  Flag; 

Filename  ;■  Name  +  *.DTA'; 

wrltelnC 'Code  will  be  stored  In  file  ', Filename); 
GenLlnear(Fllename,N,Tapl) ; 

wrltelnC 'Do  you  wish  to  generate  any  Linear  shift  codeflles?') 
wrltelnC 'If  so.  Input  Y...lf  not.  Input  N'); 
readlnCch); 

until  Cch  ■  'N')  or  Cch  “  'n'); 
end;  C*procedure  Linear  Generator*) 
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(A************************************************************************) 


(*  *) 

(*  KasamlGenerator  *) 

(*  *) 

(*4t*****************4r****4k****it*4t*********1k************4t*4r1(t****4:******4r4r4r*^ 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  A. 3  *) 

(*  DESIGN  DATE:  Jan  85  LAST  UPDATE:  9  Mar  85  *) 

(* - *) 

(*  DESCRIPTION:  This  procedure  does  most  of  the  user /program  Interaction  *) 
(*  required  to  produce  a  set  of  Kasaml  codes.  The  Length  of  the  LFSR,  *) 
(*  tap  assignments  and  Filename  for  storage  are  obtained  here.  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  None  *) 

(* - *) 

(*  INTEEINAL  VARIABLES:  Tapi — tap  assignment  for  LFSR  *) 

(*  N— length  of  LFSR  *) 

(*  Name — string  for  output  filename  w/o  ext  *) 

(*  Flag — boolean  used  for  control  *) 

(*  ch — character  used  for  user  I/O  *) 

(*  Regl — ^Reglster  array  for  LFSR  *) 

(*  Llstl,Llst2,Llst3 — pointers  to  sequences  used  to  *) 

(*  generate  Kasaml  Sequences  *) 

(*  Backl,Back2,Back3 — pointers  to  ends  of  sequences  *) 

(* - *) 

(*  CONSTANTS:  None  *) 

(* - *) 

(*  FILES  USED;  None  *) 

(* - *) 

(*  MODULES  CALLED;  LoadInput  *) 

(*  Test  *) 

(*  GenKasaml,  GetRldOf  *) 

(* - *) 

procedure  KasamlGenerator; 


var  Regl, Tapi  :  RegStructure; 

N  :  Integer; 

Name:  Word; 
ch:  char; 

Flag:  boolean; 

Llstl ,Llst2,Llst3,Backl ,Back2,Back3:  dataptr ; 

begin 

repeat 

Llstl:*  nil; 

Ll8t2:"  nil; 

Backl:*  nil; 

Back2:"  nil; 

Ll8t3:"  nil; 

Back3:*  nil; 

wrltelnC ’This  generation  uses  one  linear  sequence  and  Its 
decimation  to  produce'); 
irrlteln( 'the  desired  code  set'); 


writeln( 'Input  the  length  of  the  register... 2  -  15  please  '); 
readln(N); 

while  (N  <  2)  or  (N  >  15)  or  odd(N)  do  begin 

writelnC ' range  must  be  2  to  14  and  an  even  Integer  please'); 
readln(N) 

end;  (*while  loop*) 

LoadInput(N, Tapi, 'Taps  for  register  number  1'); 

writelnC 'What  file  do  you  want  to  store  the  sequence  generated ' ) ; 

writelnC 'No  more  than  8  characters  in  filename'); 

repeat 

Flag;"  true; 
readlnCName) ; 

TestCName,Flag) ; 
if  not  Flag  then  begin 

writelnC ' Invalid  Input  Please  try  again'); 
writelnC 'No  Special  Characters  on  Inputted  Name'); 
end;  C*if  statement*) 
tmtll  Flag; 

writelnC 'Code  will  be  stored  in  file  ',Name); 

6enKasamiCNaiiie,N, Tapi, Llstl,l<l8t2, Lists,  Backl,Back2,Back3) ; 
wrlteln; 

writelnC 'Do  you  wish  to  generate  any  other  Kasaml  codefiles?'); 
writelnC 'If  so,  input  7  if  not,  input  N'); 
readlnCch) ; 

GetRldOf CLlstl) ; 

GetRldOf CLlat2) ; 

GetRldOf  CUstS); 
until  Cch  ■  'N')  or  Cch  ■  'n'); 
end;  C*of  procedure  KasamlGenerator*) 


(*  *) 

(*  MAINPROGRAM  *) 

(*  *) 

(^*1c1i1t1t*1i1t1t1t1t1t1t1tit1tiHi1c1i1t1c1i*1t1ft1i1i1t1i1t1t*1t1t1t1i1c1t1t1i1t*it1i1t1t1Hi1fk1t1c1t1t1t1t1t1t1t1Ht1t1i1fk*1c*1t1t1i**') 

begin  (*Main  program*) 

wrltelnC  "This  program  generates  code  sequences  of  Linear  shift 
registers ' ) ; 

wrltelnC 'and  stores  them  In  designated  flies  for  further 
manipulation* ) ; 

wrltelnC 'You  have  the  option  to  generate  Hazlmal  codes  using  the 
Linear*); 

wrltelnC 'Shift  register  option  with  appropriate  Input  for  tap 
aslgnments ' ) ; 

wrltelnC 'or  you  can  produce  gold  codes  using  two  linear  shift 
registers ' ) ; 

wrltelnC 'and  Inputting  the  appropriate  taps  as  dictated  by  the 
preferred ' ) ; 

wrltelnC ' pair  polynomljLLs ' ) ; 
wrlteln; 

wrltelnC 'DNF0RTDNATEL7  the  software  does  not  catch  Improper  tap 
assignments ' ) ; 

wrltelnC '  for  ^^^l■r^nla1  or  gold  code  sequences ,  Therefore  you  must 
Insure  that'); 

wrltelnC 'your  assignments  axe  correct. . .Check  a  book  If  necessaryl ' ) ; 
wrlteln; 

wrltelnC 'Would  you  like  to  generate  any  codeflles  at  this  time?'); 

wrltelnC ' Input  Y  for  yes  or  N  for  no’); 

readlnCch); 

while  Cch  “  *Y*)  or  Cch  "  *y')  do  begin 

wrltelnC 'Which  type  of  code  would  you  like  to  generate  at  this 
time? ' ) ; 

wrltelnC ' Input  L  for  Linear  shift  register  code  generation'); 

wrltelnC ' Input  G  for  Gold  code  generation'); 

wrltelnC 'Input  K  for  Kasaml  Generation'); 

readlnCch); 

case  ch  of 

'L*  ,  '!'  ;  begin 

ClrScr; 

LlnearGenerator ; 
end; 

*G'  ,  'g*  ;  begin 

ClrScr; 

GoldCodeGenerator ; 
end; 

'K*  ,  'k*  ;  begin 

ClrScr; 

KasamlGenerator ; 
end; 


end;C*case*) 


mmmmmm 


vrriteInC 'Would  you  like  to  generate  any  more  codefiles  Y/N  ?'); 
readln(ch); 
end(*while  loop*); 
end.(*of  program  GenerateCodea*) 


1 

« 

I 
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(^*1t1c*1t1t1c1t*1e**it*1i1c1i*1t1t1t1c1i1i1c*1tit1t1c**1i1i1e1c1e1t1t**1i1t1clc1c1c1fk1s1t1t1c*1i1fk1c**1fk1c1c1t1t*1c*1t1t1c*') 


(*  *) 

(*  PartialCorrelate  *) 

(*  *) 

I*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  B.O  *) 

(*  DESIGN  DATE;  Jan  1985  LAST  UPDATE;  1  Apr  1985  *) 

(* - *) 

(*  DESCRIPTION;  This  program  will  perform  correlations  on  user  designs-  *) 
(*  code  flies  composed  of  ASCII  I's  and  O's.  The  options  offered  by  *) 


(*  the  porgram  In  a  menu  driven  fashion,  are  to  produce  the  periodic  *) 
(*  odd  and  even  correlation  functions.  Other  options  performed  are  the  *) 
(*  thresholding  of  the  Odd,  Even  functions.  This  process  will  produce  *) 
(*  data  that  can  be  plotted,  threshold  vs  amount  exceeding  threshold.  *) 
(*  The  user  is  given  the  option  of  i^ch  thresholding  process  positive,  *) 
(*  negative  or  both.  All  manlpxilatlon  of  data  Is  done  dynamically  with  *) 
(*  pointers,  so  that  the  program  needs  no  modification  when  larger  mem-  *) 
(*  ory  devices  are  used.  The  size  of  the  memory  of  the  machine  used  Is  *) 


(*  the  only  limitation  on  the  length  of  the  Input  codes..  *) 

(* - ^*) 

(*  CALLING  ARGUMENTS;  Std  Input,  Output,  Designated  Text  File  ’  *) 

(* - - *) 

(*  INTERNAL  VARIABLES;  *) 

(*  Phase — the  phase  shift  of  desired  process  *) 

<*  PctRlght,PctLeft— used  for  output  of  thresholding  values  *) 

(*  Sumleven,Sum2even,Sumlodd,Sum2odd — used  temp  storage  values  *) 

(*  for  thresholding  process  *) 

(*  Threshold — actual  threshold  passed  Into  threshold  process  *) 

(*  DataList,Back — pointers  used  to  create/manipulate  output  data  *) 

(*  Resultflle,Codedta — File  variable  for  file  I/O  *) 

(*  Result, Fllenal,Fllena2;  string  for  filename  with  extension  *) 

(*  Name — User-inputted  filename  for  storage  of  the  code  *) 

(*  ch — character  Input  to  the  program  from  user  *) 

(*  FlagA,FlagB,Flwg£,stop — booleans  used  to  control  processes  *) 

(*  performed  as  set  by  user's  response  to  menus  *) 

(*  Llstl,Llst2  used  to  identify  the  beginning  of  linked  lists  *) 

(*  representing  the  codes  to  be  correlated  *) 

(  * - *  ) 

(*  CONSTANTS;  None  *) 

(* - *) 

(*  FILES  USED;  User  designated  Filename  for  code  readln  and  data  output  *) 

(* - *) 

(*  MODULES  CALLED;  CreateList  *) 

(*  ShlftList  *) 

(*  Correlate  *) 

(*  Writellst  *) 

(*  CountUpper  *) 

(*  CountLower  *) 

(*  DlsposeOf  *) 

(*  GetRldOf  *) 

(* - *) 


program  PartialCorrelateC Input,  output); 


type 

Name  ■  string[12]; 
pointer  “  ^  datarecord; 
datarecord  = 
record 

data:  char; 
next;  pointer 
end;  (*of  record*) 
dataptr  ■  '*  info; 
info  ■ 
record 

Even;  real; 

Odd:  real; 
nextl:  dataptr 
end;  (*record*) 


var 

Phase,  Delay,  J,  N:  Integer; 

Listl,  Llst2:  pointer; 

DataLlst,  Back:  dataptr; 

Result,  Filenal,  Fllena2:  Name; 

Resultflle,  Codedta:  text; 

PctRight,  PctLeft,  PctEven,  PctOdd, 

Sufflleven,  Sum2even,  Sumlodd,  Sum2odd,  Threshold;  real 
ch:  char; 

FlagA,  FlagB,  FlagE,  stop:  boolean; 
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wrlteln( 'A — Even  function  only'); 
wrltelnC ' B — Odd  function  only ' ) ; 
writelnC'E — both  of  the  above'); 
writeln; 

writelnC ' Select  one  please'); 

FlagA  :■  false; 

FlagB  ;*  false; 
readln(ch); 
case  ch  of 

'A',  'a'; 
begin 

FlagA  :■  true; 

writelnC 'Even  function  selected') 
end; 

'B',  'b': 
begin 

FlagB  true; 

writelnC 'Odd  function  selected') 
end; 

'E',  'e'; 
begin 

FlagA  ;■  true; 

FlagB  :■  true; 

FlagE  ;■  true; 
writelnC 'All  selected*) 
end 

end;  C*ca8e*) 

writelnC 'You  have  a  choice  of  — '); 
writelnC  U — using  upper  Threshold  only*); 
writelnC  L — using  lower  Threshold  only'); 
writelnC  B — using  both  Threshold'); 
writelnC  please  select  one'); 
readlnCch) ; 

Threshold  :■  0; 
stop  f«ilse; 

writelnC  WORKING  ON  Threshold  ALGORITHM'); 
Sumleven  ;■  0; 

Sum2even  ;■  0; 

Suinlodd  0; 

Sum2odd  0; 

while  CThreshold  <■■  N)  and  not  stop  do  begin 
case  ch  of 
'U',  'u'; 
begin 

CountUpperCDataLlst,  Threshold, 
Sumleven,  Sumlodd) 


(^isititigisiiititticiticit'kit’kieiticiticititicicic’k'kieifit’k’kitificiciticltititit’kieitificic’k’kiiiiiticicisicitltit’kit'kicicicicitieititis'kicit) 

(*  *) 

(*  MAINPROGRAM  *) 

(*  *) 

(A***************************************************************************) 

begin 

wrltelnC "This  program  allows  you  to  do  correlations  of  ASCII  code  files'); 
wrltelnC 'representing  the  pseudorandom  sequences.  You  may  get  output  of); 
wrltelnCthe  actual  Odd  and  Even  correlation  function,  Right  and  Left  '); 
writeln( 'Partial  correlations  or  the  threshold  functions  of  any  or  all  of); 
wrltelnC ' of  these  four. ' ) ; 

wrltelnC 'The  Phase  of  the  codefiles  may  be  varied  from  0  to  L'); 
wrlteln; 

wrltelnC 'Do  you  want  to  do  any  correlations  at  this  time?. .Y/N'); 
readlnCch); 

If  Cch  ■  'y')  or  Cch  ■  'Y')  then  begin 
repeat 

wrltelnC ' enter  the  first  file  you  wish  to  correlate'); 
readlnCFUexial ) ; 

wrltelnC 'now  enter  the  second  Filename'); 
readlnCFUenal) ; 

wrltelnC ' Enter  phaseshlft  desired  for  run...O  to  length'); 
readlnC Phase ) ; 

wrltelnC 'What  of  the  following  options  do  you  desire?'); 
wrltelnC 'A  —  a  listing  of  actual  correlation  functions ' ) ; 
wrltelnC 'B  —  a  listing  of  Threshold  functions'); 
readlnCch); 

CreateLlstCFllenal,  Llstl,  N); 

CreateListCFllena2,  Llst2,  N); 
for  J  1  to  Phase  do  begin 
Shlf tLlst Cllst 1 ) ; 

ShlftLlstCLlst2) 
end;  C*for  loop*) 

DataLlst  nil; 

Back  nil; 

CorrelateCDataLlst,  Back,  Delay,  Llstl,  Llst2,  N); 
repeat 

wrltelnC 'Where  do  you  want  to  store  the  output?'); 

readlnCResult) ; 

asslgnCResultf lie, Result) ; 

rewrlteCResultflle) ; 

case  ch  of 

'A',  'a': 
begin 

WrltellstCDataLlst,  Resultflle) 
end; 

'B',  'b'; 
begin 


ClrScr; 

wrltelnC 'You  have  the  choice  of  the  following  listings 
for  Threshold  functions : ' ) ; 


(iticitisicit’k'kitis'kicitii'kiciciciticic'kiclciciticicic'kit'kieisicitisittis'kit'kicitftlc'kiticitieicif'kicic'kiciticiciciglclclcic’k’kicit*^ 

(*  *) 

(*  Writelist  *) 

(*  *) 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  B.4  *) 

(*  DESIGN  DATE;  Nov  1984  LAST  UPDATE;  14  Mar  1985  *) 


(*  DESCRIPTION:  This  procedure  takes  a  linked  list  of  data  and  writes  *) 
(*  the  data  fields,  even,  odd,  right  and  left  to  the  designated  file,  as*) 
(*  well  as  the  chip  location  from  0  to  L,  the  length  of  the  list.  *) 


(*  CALLING  ARGUMENTS;  List — pointer  to  front  of  datallst 
(*  OutFlle-strlng  for  output  filename 


(*  INTERNAL  VARIABLES: 

(*  P — temporary  pointer  used  to  traverse  list 

(*  K— Counter  used  as  chip  position  in  list 


(*  CONSTANTS;  None 


(*  FILES  USED;  Writes  to  OutFlle 


(*  MODULES  CALLED;  None 


procedure  Writelist (List:  dataptr;  var  OutFlle:  text); 


K:  Integer; 
P:  dataptr; 


begin 

K  :-  0; 

P  ;■  List; 

while  P  O  nil  do  begin 
wrltelnCOutFlle,  K, 

K  ;-  K  +  1; 

P  ;■  P''.nextl 
end  (*while  loop*) 
end;  (*procedure  Writelist*) 


P^.Even;  5;  2,  P^.Odd;  5;  2) 


R  :■  (-Ptrl''.Even  -  Thresh)  /  (-Ptrl'‘ .Even  +  Ptr2'' .Even) ; 
Countl  :■  Countl  +  R 

end 

end  (*lf*);  (*if  FlagA*) 
if  Flags  then  begin 

if  (-Ptrl''.0dd  <  Thresh)  and  (-Ptr2'‘.0dd  >■  Thresh)  then  begin 
R  :■  (-Ptr2''.0dd  -  Thresh)  /  (-Ptr2'‘.0dd  +  Ptrl^'.Odd); 

Count2  Count2  +  R 

end;  (*if*) 

if  (-Ptrl''.Odd  >“  Thresh)  and  (-Ptr2'‘.0dd  >■  Thresh)  then  begin 
Count2  Count2  +  1 

end;  (*if*) 

if  (-Ptrl'^.Odd  >“  Thresh)  and  (-Ptr2''.0dd  <  Thresh)  then  begin 
R  (-Ptrr.Odd  -  Thresh)  /  (-Ptrl''.Odd  +  Ptr2'*.0dd); 

Count2  Count2  +  R 

end 

end  (*if*);  (*if  FlagB*) 

Ptrl  Ptr2; 
until  Ptrl^.nextl  “  nil 
end;  (^procedure  CountLower*) 


(*  *) 

(*  CountLower  *) 

(*  *) 

(^Itlilililili*******************************************-^^:*********************'^ 

(*  DESIGNER:  Richard  C.  Gender  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  B.6  *) 

(*  DESIGN  DATE:  Jan  1985  LAST  UPDATE:  26  Jan  1985  *) 

* - * 


( 


) 


(*  DESCRIPTION:  This  procedure  t£dces  a  data  linked  list  representing  the*) 
(*  the  even  an  odd  correlation  functions,  and  determines  what  amount  of  *) 
of  each  curve  lies  below  the  negative  threshold.  It  utilizes  the  *) 

fact  that  from  one  data  record  to  the  next,  1  chip  time  elapses.  *) 

Essentially  a  similar  triangle  operation  using  the  correlation  values*) 
from  one  point  to  the  next,  and  the  threshold.  *) 


(* 

(* 

(* 

(* 


0 


(*  CALLING  ARGUMENTS:  List — pointer  to  front  of  dat5illst  *) 
(*  Thresh — threshold  value  *) 
(*  Count l,Count2 — Counters  used  to  add  amount  of  curve*) 
(*  exceeds  the  threshold  for  each  function  *) 

* _ 


( 


) 


(*  INTERNAL  VARIABLES: 

(*  Ptrl,Ptr2 — temp  pointer  used  to  do  the  process 

(*  R — temp  storage  location  used  with  cotinter 


*) 

*) 

*) 


0 


(*  CONSTANTS:  None 
(* - 


(*  FILES  USED; 
_ _ _ 


None 


( 


*) 

-*) 

*) 


(*  MODULES  CALLED; 

(* 


0 


None 


procedure  CountLowerCLlst:  dataptr;  Thresh:  real; 

var  Countl,  Count2:  real); 


*) 

-*) 


var 


Ptrl,  Ptr2;  dataptr; 
R:  real; 


begin 

Countl  ;■  0; 

Count2  0; 

Ptrl  ;■  List; 
repeat 

Ptr2  ;■  Ptrl''.nextl; 

If  FlagA  then  begin 

If  (“Ptrl'*. Even  <  Thresh)  and  (-Ptr2''.Even  >■  Thresh)  then  begin 
R  ;■  (-Ptr2''.Even  -  Thresh)  /  (-Ptr2''.Even  +  Ptrl''. Even); 
Covintl  ;■  Countl  +  R 
end;  (*lf*) 

If  (“Ptrl''. Even  >■  Thresh)  and  (-Ptr2''.Even  >■  Thresh)  then  begin 
Countl  Countl  +  1 
end;  (*if*) 

If  (“Ptrl''. Even  >“  Thresh)  and  (-Ptr2''.Even  <  Thresh)  then  begin 
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R  ;■  (Ptrl''.Eveii  -  Thresh)  /  (Ptrl''.Eveii  -  Ptr2^.Even) ; 
Countl  Countl  +  R 
end;  (*lf*) 
end;  (*if  Flag  A*) 
if  FlagB  then  begin 

if  (Ptrl^'.Odd  <  Thresh)  and  (Ptr2'‘.0dd  >“  Thresh)  then  begin 
R  (Ptr2'‘.0dd  -  Thresh)  /  (Ptr2'‘.0dd  -  Ptrl'‘.Odd); 

Count2  ;■  Count2  +  R 
end;  (*if*) 

if  (Ptrl''.Odd  >*■  Thresh)  and  (Ptr2'‘.0dd  >“  Thresh)  then  begin 
Count2  :■  Count2  +  1 
end;  (*if*) 

if  (Ptrl''.Odd  >“  Thresh)  and  (Ptr2'*'.0dd  <  Thresh)  then  begin 
R  (Ptrl^.Odd  -  Thresh)  /  (Ptrl^.Odd  -  Ptr2''.0dd); 

Cotmt2  ;■  Count2  +  R 

end 

end  (*if*);  (*if  FlagB*) 

Ptrl  Ptr2; 
until  Ptrl'“.nextl  ■  nil 
end;  (*procedure  CountUpper*) 


(*  *) 

(*  CountUpper  *) 

C*  *) 

(^*1tit1t1t1t1fk1t1t1fit1t1t1c1c1t1tit1i1t1iifk1t1i1e1t1c1t1t1i1t1cic1t1c*1fkit1t1t1fk1c1i'k1i1i1t1t1i1i1t1i1fk1i1fk1c1i1t1i1t1c1t1tic1i1i1i'^ 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  B.5  *) 

(*  DESIGN  DATE:  Jan  1985  LAST  UPDATE;  26  Jan  1985  *) 


(*  DESCRIPTION:  This  procedure  takes  a  data  linked  list  representing  the*) 
(*  the  even  and  odd  correlation  fimctlons,  and  determines  what  amount  of*) 
(*  each  function  lies  above  the  Inputted  threshold.  It  utilizes  the  *) 
(*  fact  that  from  one  data  record  to  the  next,  1  chip  time  elapses.  *) 
(*  Essentially  a  similar  triangle  operation  using  the  correlation  values*) 
(*  from  one  point  to  the  next,  and  the  threshold.  *) 


(*  CALLING  ARGUMENTS:  List — pointer  to  front  of  datallst  *) 
(*  Thresh — threshold  value  *) 
(*  Countl,Count2 — Counters  used  to  add  amount  of  curve*) 
(*  exceeds  the  threshold  for  each  function  *) 


(*  INTERNAL  VARIABLES: 

(*  Ptrl,Ptr2 — temp  pointer  used  to  do  the  process 

(*  R— temp  storage  location  used  with  counter 


(*  CONSTANTS:  None 


(*  FILES  USED;  None 


(*  MODULES  CALLED:  None 


procedure  CountUpper (List:  dataptr;  Thresh:  real; 

var  Countl,  Count2:  real); 


var 

Ptrl,  Ptr2;  dataptr; 

R:  real; 

begin 

Countl  :*  0; 

Count2  ;■  0; 

Ptrl  :■  List; 
repeat 

Ptr2  Ptrl^'.nextl; 

If  FlagA  then  begin 

if  (Ptrl^.Even  <  Thresh)  and  (Ptr2''.Even  >■  Thresh)  then  begin 
R  :■  (Ptr2''.Even  -  Thresh)  /  (Ptr2''.Even  -  Ptrl''. Even); 
Countl  ;■  Cotmtl  +  R 
end;  (*lf*) 

If  (Ptrl''. Even  >■  Thresh)  and  (Ptr2''.Even  >■  Thresh)  then  begin 
Countl  :■  Countl  +  1 
end;  (*lf*) 

If  (Ptrl''. Even  >■  Thresh)  and  (Ptr2*.Even  <  Thresh)  then  begin 


if  Ptrl‘'.data  ■  Ptr2''.data  then 
Agree  ;■  Agree  +  1 
else 

DlsAgree  DlsAgree  +  1; 
Ptrl  Ptrl^.next; 

Ptr2  ;•  Ptr2'“.next 
end;  (*whlle  Loop*) 

RightCor  :■  Agree  -  DlsAgree; 

Agree  :■  0; 

DlsAgree  0; 

Ptrl  :■  Llstl; 

while  Ptr2  O  nil  do  begin 

if  Ptrl''. data  "  Ptr2''.data  then 
Agree  ;■  Agree  +  1 
else 

DlsAgree  DlsAgree  ■¥  1; 
Ptrl  :■  Ptrl''. next; 

Ptr2  ;■  Ptr2^.next 
end;  (*whlle  Loop*) 

LeftCor  :■  Agree  -*  DlsAgree; 

Even  ;■  LeftCor  +  RightCor; 

Odd  :■  RightCor  -  LeftCor; 
EnqueCDatalilst,  Back,  Odd,  Even); 

K  K  +  1; 

Delay  ;■  Delay  +  1 
end;  (*whlle  loop*) 

Agree  ;•  0; 

DlsAgree  :■  0; 

Ptrl  ;■  Llstl; 

Ptr2  :*  Llst2; 

while  Ptr2  O  nil  do  begin 

If  Ptrl''. data  "  Ptr2''.data  then 
Agree  ;■  Agree  +  1 
else 

DlsAgree  DlsAgree  +  1; 

Ptrl  ;■  Ptrl". next; 

Ptr2  Ptr2‘'.next 
end;  (*whlle  Loop*) 

LeftCor  Agree  -  DlsAgree; 

Even  :*  LeftCor; 

Odd  :■  -LeftCor; 

RightCor  :■  0; 

EnqueCDataLlst,  Back,  Odd,  Even) 
end;  (*procedure  correlate*) 


(^mt***********************************************************************) 

(*  *) 

(*  Correlate  *) 

(*  *) 

(*************************************************************************) 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  B.l  *) 

(*  DESIGN  DATE;  Jan  1985  LAST  UPDATE;  14  Jan  1985  *) 


(*  DESCRIPTION;  This  procedure  takes  the  2  Linked  Lists  of  characters  *) 
(*  representing  codesets  and  determines  the  even,  odd,  right  and  left  *) 
(.*  correlation  functions  at  discrete  values.  It  places  the  determined  *) 
(*  values  at  corresponding  positions  In  a  third  linked  list  *) 


(*  CALLING  ARGUMENTS;  DataLlst ,Back — pointers  to  front  and  rear  of  output*) 
(*  Datallst  *) 
(*  Dela3r-offset  of  lists  during  correlation  cycle(TAU)*) 
C*  Llstl,Llat2 — pointers  for  front  of  codellsts  *) 
(*  N — ^Length  of  List  *) 


(*  INTERNAL  VARIABLES;  *) 
(*  Ptrl,Ptr2 — temp  pointer  used  to  do  the  process  *) 
(*  LeftCor,  RlghtCor — partial  corelatlons  values  used  to  de-  *) 
(*  termlne  the  even  and  odd  correlation  functions  *) 
<*  Even, Odd— discrete  correlation  values  *) 
(*  Agree, Disagree — used  as  counters  during  correlation  process  ♦) 
(*  Count ,K— counters  used  In  process  *) 


(*  CONSTANTS;  None 


(*  FILES  USED; 


(*  MODULES  CALLED;  Enque 


procedure  Correlate (var  DataLlst,  Back;  dataptr;  Delay:  Integer; 

Llstl:  pointer;  Llst2:  pointer;  N:  Integer); 


LeftCor,  RlghtCor,  Even,  Odd:  real; 
Count,  Agree,  DlsAgree,  K:  Integer; 
Ptrl,  Ptr2:  pointer; 

begin 

K  ;-  1; 

while  K  <■  N  do  begin 
Ptrl  ;■  Llstl; 

Ptr2  :■  Llst2; 

Agree  :■  0; 

DlsAgree  :■  0; 
for  Count  :*  2  to  K  do 
Ptrl  ;-  Ptrl^.next; 
while  Ptrl  <>  nil  do  begin 
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(*  *) 

(*  ShiftLlat  *) 

(*  *) 


(*  DESIGNER:  Richard  C.  Gonder 
(*  ADVISOR:  MAJ  Ren  Castor 
(*  DESIGN  DATE:  Jan  1985 


MASTER  THESIS  WORK 
MODULE  NUMBER  B.3 
LAST  UPDATE;  14  Jan  1985 


(*  DESCRIPTION:  This  procedure  merely  takes  the  first  record  of  the 
(*  linked  list  and  places  it  at  the  end  of  the  list. 


(*  CALLING  ARGUMENTS:  List — pointers  to  front  of  list 


(*  INTERNAL  VARIABLES; 

(*  PjPtr — temp  pointer  used  to  do  the  process 


(*  CONSTANTS:  None 


(*  FILES  USED 


(*  MODULES  CALLED 


procedure  ShlftLlstCvar  List:  pointer) 


First,  P:  pointer; 


begin 

First  List; 

P  ;■  List; 

while  P'^.next  O  nil  do 
P  :■  P^'.next; 

First  ;■  First'*. next; 


P  .next  ;■  List; 

List'*. next  :•  nil; 

List  :■  First 

end;  (*procedure  ShlftLlst*) 


(*  *) 

(*  GetRidOf  *) 

(*  *) 

(*************************************************************************) 


(*  DESIGNER:  Richard  C.  Gonder 
(*  ADVISOR:  MAJ  Ren  Castor 
(*  DESIGN  DATE:  Jan  1985 


MASTER  THESIS  WORK 
MODULE  NUMBER  B.8 
LAST  UPDATE:  14  Jan  1985 


(*  DESCRIPTION:  This  procedure  releases  a  linked  list  from  memory 
(*  allocation.  Specifically  used  on  output  data  lists  destruction 


(*  CALLING  ARGUMENTS:  List — pointers  to  front  of  list 


(*  INTERNAL  VARIABLES: 

(*  PjPtr — temp  pointer  used  to  do  the  process 


(*  CONSTANTS:  None 


(*  PILES  USED 


(*  MODULES  CALLED 


procedure  GetRldOf(var  List:  dataptr) 


P,  ptr:  dataptr; 

begin 

P  :■  List; 
ptr  :■  P; 

while  P'^.nextl  O  nil  do  begin 
ptr  :■  P''.nextl; 
dlspose(P); 

P  :■  ptr 

end  (*^rtille  loop*); 
dlspose(P); 

List  :■  nil 

end;  (*procedure  GetRidOf*) 


(*  *) 

(*  DlaposeOf  *) 

(*  *) 

^1t1t1i1clt1t1c1e1t1t1t1i1t1t1i1t1tielc1i1cic1c1i1t1fk1c1t1t1iii1fk1t±1t1cic1i1t1i*1i1t1i1tit1tit1i1t1c1iiiii1t1i1i1i1i1i1t1t1i1c1i1t1i1t1i1i1i'^ 


(*  DESIGNER:  Richard  C.  Gonder 
(*  ADVISOR;  MAJ  Ken  Castor 
(*  DESIGN  DATE;  Jan  1985 


MASTER  THESIS  WORK 
MODULE  NUMBER  B.7 
LAST  UPDATE:  14  Jan  1985 


(*  DESCRIPTION:  This  procedure  releases  a  linked  list  from  memory 
(*  allocation.  Specifically  used  on  Input  data  lists  destruction 


(*  CALLING  ARGUMENTS;  List — pointers  to  front  of  list 


(*  INTERNAL  VARIABLES: 

(*  PjPtr — temp  pointer  used  to  do  the  process 


(*  CONSTANTS:  None 


(*  FILES  USED 


(*  MODULES  CAT.T.m 


procedure  DlaposeOf (var  List:  pointer) 


P,  ptr:  pointer; 

begin 

P  :■  List; 
ptr  :■  P; 

while  P'^.next  O  nil  do  begin 
ptr  :•  P'‘.next; 
dlspose(P); 

P  ;■  ptr 

end  (*whlle  loop*); 
dlspose(P); 

List  ;■  nil 

end;  (*procedure  DlsposeOf*) 


^*****:<t1lr*:(lr*ir*4r***1^ii^4r*i^llr*1lr***1lr1k******1^***!(^*******^4r*4r*1^****4r1k4r*1lr*Kr1ir1lri(^*:/^*4r*il^^ 

(*  *) 

(*  E  n  q  u  e  *) 

(*  *) 

(^1i1t1c1t1c1i1i1t1t1t1t1t1c1s1i1t1i1t1t1cit1i1t1t1t1c1t1tifk1t1i1i1i1t1t1t1t1t1t1c1t*1t1t1t1i1i1t1c1c1c1t1t1t1c1t1i1c1t1c1i1t1c1t*1i1fk1t1i1t1t'^ 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  B.1-1  *) 

(*  DESIGN  DATE;  Nov  1984  LAST  UPDATE;  14  Dec  1984  *) 


(*  DESCRIPTION:  This  procedure  enque  a  data  record  representing  the 
(*  values  of  the  even,  odd,  right  and  left  correlation  functions  at  a 
(*  discrete  V£d.ue  In  an  ordered  linked  list. 


(.*  CALLING  ARGUMENTS:  List, Back — pointers  to  front  and  rear  of  list 
(*  Odd, Even — values  of  correlation  functions 


(*  INTERNAL  VARIABLES; 

(*  P — tenp  pointer  used  to  do  the  process 


(*  CONSTANTS:  None 


(*  FILES  USED 


(*  MODULES  CALLED 


procedure  Enque(var  List,  Back:  dataptr;  Odd,  Even:  real) 


P:  dataptr; 


begin 

If  List  *  nil  then  begin 
new(P) ; 

List  ;■  P; 

P^'.Even  :«  Even; 
P^'.Odd  Odd; 

Back  ;■  P; 

Back''. next  1  nil 

end  else  begin 
new(P); 

P'.Even  ;■  Even; 
P'.Odd  Odd; 

Back'.neztl  ;■  P; 

Back  ;■  P; 

Back'. next 1  ;■  nil 
end  (*if  then  else*) 
end;  (*procedure  Enque*) 


(*  *) 

(*  CreateLlst  *) 

(*  *) 

(^*it1t1c1t1t1t1c1t1c*1t1t*1t1t***1c1t'?!-****1tic1t1iit1c1t1tifk1t*1i1t1c1c1i1i1i1t1i1i*1i1c1i1i1t1i1t1i*1i**1c1t1t*1t1fk1t1c1i1i*') 

(*  DESIGNER:  Richard  C*  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  B.2  *) 

(*  DESIGN  DATE:  Jan  1985  LAST  UPDATE:  4  Jan  1985  *) 

(* - *) 

(*  DESCRIPTION:  This  procedure  will  read  in  a  file  of  char,  placing  each*) 

(*  character  into  a  linked  list.  It  terminates  at  the  EOF  of  the  file  *) 
(*  read.  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  Filename — string  of  char  *) 

(*  List — pointer  to  the  beginning  of  linked  list  *) 

(*  K — number  of  characters  read  into  list  *) 

( * - *  ) 

(*  INTERNAL  VARIABLES:  *) 

(*  Pjptr — used  to  construct  list  *) 

(*  bit — temp  location  for  character  read  *) 

(* - *) 

(*  CONSTANTS:  None  *) 

(* - *) 

(*  FILES  USED:  User  designated  Filename  for  code  storage  *) 

(  * - *  ) 

(*  MODULES  CALLED;  None  *) 

- *) 


procedure  CreateLlst (Filename:  Name;  var  List:  pointer;  var  K:  Integer); 


P,  ptr;  pointer; 
bit:  char; 

begin 

asslgnCCodedta, Filename) ; 
reset ( Coded ta ) ; 
new(P) ; 

List  P; 
ptr  P; 

K  0; 

while  not  eof(Codedta)  do  begin 
readCCodedta,  bit); 

K  K  +  1; 

P'‘.data  ;■  bit; 
ptr '‘.next  ;*  P; 
ptr  ;■  P; 
new(P) 

end;  (*whlle  loop*) 
close (Codedta) ; 
ptr''. next  ;■  nil; 
end;  (*procedure  CreateLlst*) 


•L',  ’1’; 
begin 

CountLower (DataLlst ,  Threshold , 

Suffl2even,  Sum2odd) 

end; 

else 

CountUpper(DataLlst,  Threshold, 

Suffileven,  Sumlodd) ; 
CountLower(DataLlst,  Threshold, 

Suffl2even,  Sum2odd) 

end;  (*case*) 

PctEven  ;■  (Sumleven  +  Su]n2even)  /  N; 

PctOdd  :■  (Sumlodd  +  Sum2odd)  /  N; 
wrlte(Resultflle,  Threshold  /  N:  5:  4); 

If  FlagA  then 

write (Resultf lie,  PctEven:  8:  7); 

If  Flags  then 

write (Resultf lie,  PctOdd:  8:  7); 

write (Threshold  /  N:  5:  4); 

If  FlagA.  then 

wrlteC,*,  PctEven:  8:  7); 

If  Flags  then 

wrlteC,',  PctOdd:  8:  7); 
wrlteln; 

wrlteln(Resultflle) ; 

If  Threshold  /  N  >  0.10  then 

Threshold  :■  Threshold  +  N  /  100 
else 

Threshold  :■  Threshold  +  N  /  500; 

If  not  FlagE  and  FlagA.  and  (PctEven  ■  0.0)  then 
atop  :■  true; 

If  not  FlagE  and  FlagS  and  (PctOdd  *■  0.0)  then 
stop  :■  true; 

If  FlagE  and  (PctEven  ■  0.0)  and  (PctOdd  ■  0.0)  then 
stop  :■  true 


end 

end 

end  (*whlle  loop*)  (*E  -option*);  (*case*) 
close(Resultflle) ; 

wrltelnCDo  you  want  any  other  files  listed  with  this 
dataset  of  codes?'); 

wrlteln( 'What  of  the  following  options  do  you  desire? ' ) ; 
wrlteln( 'A  —  a  listing  of  all  actual  correlation  f\mctlons ' ) ; 
wrltelnCS  —  a  listing  of  Threshold  functions'); 
wrlteln( 'N  —  None,  You  want  to  leave  this  mode'); 
readln(ch) 

until  (ch  ■  'n')  or  (ch  ■  'N'); 


DlsposeOf (Llst2) ; 

DlsposeOf (Llstl) ; 

GetRldOf (DataLlst) ; 

wrltelnCdo  you  wish  to  do  any  more  correlations? .  .Y/N' ) ; 
readln(ch) 

until  (ch  ■  'N')  or  (ch  ■  'n') 
end  (*lf  y/n*) 
end.  [  PartlalCorrelate  } 


(*  *) 

(*  PhaseCorrelate  *) 

(*  *) 

(**************4t4t*****ik***********************1(r****1k4r*******4t*it****4r****4r*^ 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  D.l  *) 

(*  DESIGN  DATE:  Jan  1985  LAST  UPDATE:  1  Apr  1985  *) 

(* - *) 


(*  DESCRIPTION:  This  program  will  perform  correlations  on  user  designs-  *) 
(*  code  files  composed  of  ASCII  I's  and  O's,  and  determine  the  minimum  *) 
(.*  and  maximum  values  of  the  Odd  correlation  function  for  L  different  *) 
(*  phases  of  the  code  file.  For  each  phase,  this  program  will  shift  *) 
(*  both  code  files  and  compute  the  Odd  function.  It  then  determines  the*) 
(*  minimum  and  maximum  values  of  Odd  fimctlon  for  that  phase,  and  plases*) 
(*  this  Information  Into  a  list.  Once  all  the  phases  have  been  cycled  *) 
(*  through,  a  listing  of  all  of  the  phases  where  the  maximum  and  minimum*) 


(*  spreads  of  the  Odd  function,  as  well  as  their  mar  and  min  values.  Is  *) 
(*  written  to  a  user  designated  file.  It  will  perform  off-peak  autocor-  *) 
(*  relation  and  total  cross-correlation  of  codes  as  determined  by  user  *) 

(*  Input.  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  Std  Input,  Output,  Designated  Text  File  *) 

(* - *) 

(*  INTERNAL  VARIABLES:  *) 

(*  J~uaed  as  counter  for  control  *) 

(*  N— Length  of  codefiles  in  characters  *) 

(*  Phasemax,Pha8emln — Value  of  Maximum  and  Minimum  spread  for  *) 

(*  for  Odd  fxinction  for  all  phase  shifts  *) 

(*  Look, DataLlst, Back — pointers  used  to  create/ manipulate  data  *) 

(*  lists  to  obtain  outputted  data  file  *) 

(*  Resultflle,Codedta — File  variable  for  file  I/O  *) . 

(*  Result, Fllenal,Fllena2:  string  for  filename  with  extension  *) 

(*  ch — character  Input  to  the  program  from  user  *) 

(*  FlagA — boolean  used  to  control  process  *) 

(*  Llstl,Llst2  used  to  Identify  the  beginning  of  linked  lists  *) 

(*  representing  the  codes  to  be  correlated  *) 

(*  K,S— used  to  set  correlation  process  for  off-peak-autocorrel-  *) 

(*  atlon  or  complete-cross-correlation  *) 

(* - *) 

(*  CONSTANTS:  None  *) 

(* - *) 

(*  FILES  USED:  User  designated  Filename  for  code  readln  and  data  output  *) 

(* - *) 

(*  MODULES  CALLED:  CreateLlst  *) 

(*  ShiftList  *) 

(*  Correlate  *) 

(*  DisposeOf  *) 

(*  GetRidOf  *) 

(*  Maxlmlnlmlze  *) 

(* - *) 
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program  PhaseCorrelateClnput, output) ; 


type  Name  ■  strliig[12]; 

pointer  ■  '^datarecord ; 

datarecord  ■  record 

data  :  char; 
next  :  pointer; 
end;  (*of  record*) 
dataptr  ■  '“Info; 

Info  ■  record 

MlnOdd  :  real; 
MazOdd  :  real; 
neztl  :  dataptr; 
end  (*record*); 

var  K,J,N,S:  Integer; 

Llstl,Ll8t2:  pointer; 

Look, DataLlst, Back:  dataptr; 
Restilt,Fllenai,Fllena2:  Name; 
Resultflle, Coded ta:  Text; 
PhasemaXjP^semln:  real; 
ch:  char; 

FlagA:  boolean; 


twji 


Vo-'-' 


(*  *) 

(*  CreateLis-t  *) 

(*  *) 

^***4k**4rlk*****4rik**4k*****1lt*it***lt******************A4t*4r4t*******ilt****4r*******^ 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  D.2  *) 

(*  DESIGN  DATE;  Jan  1985  LAST  UPDATE;  4  Jan  1985  *) 

(* - *) 

(*  DESCRIPTION;  This  procedure  will  read  In  a  file  of  char,  placing  each*) 
(*  character  Into  a  linked  list.  It  terminates  at  the  EOF  of  the  file  *) 

(*  read.  *) 

(* - *) 

(*  CALLING  ARGUMENTS;  Filename — string  of  char  *) 

(*  List — pointer  to  the  beginning  of  linked  list  *) 

(*  K — number  of  characters  read  Into  list  *) 

(* - *) 

(*  INTERNAL  VARIABLES;  *) 

(*  PjPtr — used  to  construct  list  *) 

(*  bit — temp  location  for  character  read  *) 

(* - *) 

(*  CONSTANTS;  None  *) 

(* - *) 

(*  FILES  USED;  User  designated  Filename  for*  code  storage  *) 

(* - *) 

(*  MODULES  CALLED;  None  *) 

(  * - *  ) 

procedure  CreateLlst( Filename:  Name;  var  List;  pointer;  var  K:  Integer); 

var 

P,  Ptr;  pointer; 

bit:  char; 

begin 

asslgn( Codedta , Filename ) ; 

reset(Codedta); 

new(P); 

List  ;■  P; 

Ptr  ;-  P; 

K  ;-  0; 

while  not  eof (Codedta)  do  begin 
readCCodedta,  bit); 

K  ;-  K  +  1; 

P^'.data  ;■  bit; 

Ptr'*. next  ;■  P; 

Ptr  ;■  P; 
new(P) 

end;  (*whlle  loop*) 

close(Codedta) ; 

Ptr‘'.next;"  nil; 

end;  (*procedure  CreateLlst*) 
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(*  *) 

(*  E  n  q  u  e  *) 

(*  *) 

(**4t******4t**4r************4t4r***ir**ik**********4t*4r******it*******1lt***********^ 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  D.4.1  *) 

(*  DESIGN  DATE;  Nov  1984  LAST  UPDATE;  14  Dec  1984  *) 


(*  DESCRIPTION;  This  procedure  enques  a  data  record  representing  the  *) 
(*  values  of  the  Minimum  and  Maximum  Odd  values  for  a  phase  of  the  Odd  *) 
(*  correlation  function.  The  list  Is  ordered  by  phase,  from  0  to  L.  *) 


(*  CALLING  ARGUMENTS;  List, Back — pointers  to  front  and  rear  of  list  *) 
(*  MlnOdd,MaxOdd — min  and  max  values  of  Odd  function  *) 
(*  for  a  particular  phase  *) 


(*  INTERNAL  VARIABLES: 

(*  p — temp  pointer  used  to  do  the  process 


(*  CONSTANTS:  None 


(*  PILES  USED 


(*  MODULES  CALLED 


procedure  Enque(var  List,  Back:dataptr;  MlnOdd,  MaxOdd:real) 
var  P:  dataptr; 


begin 

If  List  ■  nil  then  begin 
New(P); 

List:"  P; 

P'‘.MaxOdd;-  MaxOdd; 
P'‘.MlnOdd;-  MlnOdd; 
Back;-  P; 

Back'‘,nextl;"  nil; 

end 

else  begin 
New(P); 

P''. MaxOdd:"  MaxOdd; 

P" .MlnOdd;"  MlnOdd; 
Back^'.nextl;"  P; 
Back:"  P; 

Back''.nextl;"  nil; 
end  (*lf  then  else*) 
end;  (*procedure  Enque*) 


(*  *) 

(*  DisposeOf  *) 

(*  *) 

(**********************************«**************************************) 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  D.5  *) 

(*  DESIOJ  DATE:  Jan  1985  LAST  UPDATE:  14  Jan  1985  *) 

(  * - *  ) 

(*  DESCRIPTION:  This  procedure  releases  a  linked  list  from  memory  *) 

(*  allocation.  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  List — pointers  to  front  of  list  *) 

(* - *) 

(*  INTERNAL  VARIABLES:  *) 

(*  PjPtr — temp  pointer  used  to  do  the  process  *) 

(* - *) 

(*  CONSTANTS:  None  *) 

(* - *) 

(*  FILES  USED:  None  *) 

(* - *) 

(*  MODULES  CALLED:  None  *) 

(* - *) 


procedure  DisposeOf (var  List: pointer); 

var  PjPtr: pointer; 

begin 

P:-  List; 

Ptr;“P; 

while  P^'.nezt  O  nil  do  begin 
Ptr:-  P''.next; 
dlspose(P) ; 

P:“Ptr; 

end  (*whlle  loop*); 
dlspose(P); 

List:**  nil; 

end ; ( *procedure  DisposeOf* ) 
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(*  *) 

(*  GetRldOf  *) 

(*  *) 


(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  D.6 

(*  DESIGN  DATE:  Jan  1985  LAST  UPDATE:  14  Jan  1985 

(* - 

(*  DESCRIPTION:  This  procedure  releases  a  linked  list  from  memory 
(*  allocation. 

(* - 

(*  CALLING  ARGUMENTS:  List — pointers  to  front  of  list 

(* - 

(*  INTERNAL  VARIABLES: 

(*  P,Ptr — temp  pointer  used  to  do  the  process 

(* - 

(*  CONSTANTS;  None 

(* - - 

(*  FILES  USED;  None 

(* - 

(*  MODULES  CAUSED;  None 

(* - 


I  tp 


procedure  GetRldOf(Llst:dataptr): 

var  P,Ptr:dataptr; 

begin 

P;"  List; 

Ptr:“P; 

while  P^.nextl  O  nil  do  begin 
Ptr;*  P''.nextl; 
dlsposeCP) ; 

P;“Ptr; 

end  (*whlle  loop*); 
dlsposeCP); 

List:"  nil; 

end; (*procedure  GetRldOf*) 


•l**  ' 
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(*  *) 

(*  ShlftLlst  *) 

(*  *) 

(^lililtitltltltltlclilileltltlilutliltlnslcltltltltltlclelt*******************************************'^ 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  D.3  *) 

(*  DESIGN  DATE:  Jan  1985  LAST  UPDATE:  14  Jan  1985  *) 

(* - *) 

(*  DESCRIPTION:  This  procedure  merely  takes  the  first  record  of  the  *) 

(*  linked  list  and  places  it  at  the  end  of  the  list.  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  List — pointers  to  front  of  list  *) 

(* - *) 

(*  INTERNAL  VARIABLES:  *) 

(*  P, First — temp  pointer  used  to  do  the  process  *) 

(* - *) 

(*  CONSTANTS:  None  *) 

(* - - » - *) 

(*  FILES  USED:  None  *) 

(* - *) 

(*  MODULES  CALLED:  None  *) 

(* - *) 


procedure  ShiftLlstCvat  List: pointer); 
var  FirstjP: pointer; 
begin 

First:*  List; 

P:*  List; 

while  P^.next  O  nil  do 
P:*  P'^.next; 

First:*  First''. next; 

P'*.next:*  List; 

List''. next:*  nil; 

List:*  First; 

end;  (*procedure  ShlftLlst*) 


■w*.' 
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(*  *) 

(*  MlniMax  *) 

<*  *) 

(*ik****4t*****4t******4r******************************4t****4r********4r4r4r****4r*) 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  D.4.2  *) 

(*  DESIGN  DATE:  Jan  1985  LAST  UPDATE:  14  Jan  1985  *) 

(  * - *  ) 

(*  DESCRIPTION:  This  procedure  merely  takes  In  the  values  min,  max  and  *) 

(*  Temp,  and  sets  Min  or  Max  to  Temp  when  It  Is  less  than  Min  or  greater*) 

(*  than  Mm.  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  Temp — Input  value  *) 

(*  Min — Input  set  to  temp  when  temp  Is  less  than  MLn  *) 

(*  Max — Input  set  to  temp  when  temp  Is  greater  than  *) 

(* - *) 

(*  INTERNAL  VARIABLES:  None  *) 

(* - *) 

(*  CONSTANTS:  None  *) 

(* - *) 

(*  PILES  USED:  None  *) 

(* - ^ - *) 

(*  MODULES  CALLED:  None  *) 

(* - *) 

procedure  Mlnlmax(Temp:real;var  Mln,Max:real); 


begin 

If  Temp  <  Min  then 
Mln:“  Temp; 

If  Temp  >  Max  then 
Max:"  Temp; 

end;  (*procedure  Mlnlmax*) 


(*  *) 

(*  Correlate  *) 

(*  *) 

(*  DESIGNER:  Richard  C.  Gonder  MSTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  D.4  *) 

(*  DESIGN  DATE:  Jan  1985  LAST  UPDATE:  14  Jan  1985  *) 

(* - *) 

(*  DESCRIPTION:  This  procedure  takes  the  2  Linked  Lists  of  characters  *) 

(*  representing  codesets  and  determines  the  odd  correlation  function  at  *) 

(*  discrete  values.  It  determines  the  maximum  and  minimum  Odd  fxmctlon  *) 

(*  values  and  places  these  values  In  a  linked  list  of  records.  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  Start, Stop — values  setting  region  over  which  cor-  *) 

(*  relation  Is  perfoirmed.  *) 

(*  Llstl,Llst2 — pointers  for  front  of  codellsts  *) 

(*  N — Length  of  List  *) 

(* - *) 

(*  INTERNAL  VARIABLES:  *) 

(*  Ptrl,Ptr2 — temp  pointer  used  to  do  the  process  *) 

(*  LeftCor,  RlghtCor — partial  corelatlons  values  used  to  de-  *) 

(*  termlne  odd  correlation  value  *) 

(*  Odd— discrete  correlation  values  *) 

(*  Agree, Disagree — used  as  counters  during  correlation  process  *) 

(*  Count.,K — counters  used  In  process  *) 

(*  MlnOdd,MaxOdd — represent  mln/max  Odd  value  for  each  correl-  *) 

(*  atlon  phase  *) 

(* - *) 

(*  CONSTANTS:  None  *) 

(* - *) 

(*  FILES  USED;  None  *) 

(* - *) 

(*  MODULES  CALLED:  Enque  *) 

(*  Mini  max  *) 

(* - *) 


procedure  Correlate(Llstl:polnter;  Llst2:polnter;  N:lnteger; 

var  Start , Stop : Integer) ; 


var  MlnOdd,MaxOdd,Leftcor,Rlghtcor,Odd, :real; 
Count , Agree , D1 sAgr ee , K : Integer ; 

Ptr 1 , Ptr2 : pointer ; 

begin 

MlnOdd;-  N; 

MaxOdd:"  -N; 

while  Start  <~  Stop  do  begin 
Ptrl:"  Llstl; 

Ptr2:-  Ll8t2; 

Agree:"  0; 

DlsAgree:"  0; 


for  Count:*  2  to  Start  do 
Ptrl;*  Ptrl''.next; 
while  Ptrl  O  nil  do  begin 

if  Ptrl*. data  ■  Ptr2*.data  then 
Agree:*  Agree  +  1 
else  DlsAgree:*  DisAgree  +  1; 
Pt r 1 : *Ptr 1 * . next ; 
Ptr2:*Ptr2*.next; 
end;(*whlle  Loop*) 

RlghtCor:*  (Agree-DlsAgree); 
Agree:*  0; 

DlsAgree:*  0; 

Ptrl:*  Llstl; 

while  Ptr2  O  nil  do  begin 

if  Ptrl*. data  *  Ptr2*.data  then 
Agree:*  Agree  +  1 
else  DlsAgree:*  DlsAgree  -t-  1; 
Ptrl :*Ptrl*. next; 

Ptr2 : *Ptr2* . next ; 
end;(*whlle  Loop*) 

LeftCor:*  (Agree-DlsAgree); 

Odd:*  RlghtCor  -  LeftCor  ; 
MlnlnaxCOdd ,MlnOdd ,MaxOdd) ; 

Start :*S tart  +  1; 
end;  (*whlle  loop*) 

Enque(DataLlst  >Back ,MlnOdd ,MaxOdd) ; 
end ; ( *procedure  correlate ) 
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(*  *) 

(*  Maxlmialmize  *) 

(*  *) 

(*******************  ****4t*4r****************jlr**:k******4t****1lr****4t****4r4t***ik^ 

(*  DESIGNER;  Richard  C.  Gonder  M/^TER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  D.7  *) 

(*  DESIGN  DATE;  Jan  1985  LAST  UPDATE:  12  Feb  1985  *) 

(* - *) 

(*  DESCRIPTION:  This  procedure  takes  a  linked  list  of  min  and  max  values*) 
(*  representing  each  phase  shifted  correlation  function,  and  determines  *) 
(*  the  maximum  and  minimum  difference  for  «n  phases.  It  will  write  *) 

(*  these  values  to  the  screen.  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  DataLlst — pointer  to  front  of  list  of  records  *) 

(*  PhasemaXjPhasemln — output  value  of  Mlntumm  and  *) 

(*  maximum  spread.  *) 

(*  N — Length  of  list  *) 

(* - *) 

(*  INTERNAL  VARIABLES:  *) 

(*  Ptr — temp  pointer  used  to  do  the  process  .  *) 

(* - *) 

(*  CONSTANTS:  None  *) 

(* - : - *) 

(*  FILES  USED:  None  *) 

(* - *) 

(♦  MODULES  CALLED:  None  *) 

(* - *) 


procedure  Maxlmlnlmlze(Datallst:dataptr;  var  Phasemax,  Phasemln:real; 

N: Integer) ; 


var  Ptr:  dataptr; 
begin 

Phasemax  :"^N; 

Phasemln:*  N; 

Ptr:“datalist; 

while  Ptr next 1  O  nil  do  begin 

if  (Ptr'^.MaxOdd  -Ptr'' .MlnOdd)  >  Phasemax  then 
Phasemax:"  Ptr''.MaxOdd-Ptr''.MlnOdd; 

If  (Ptr''.MaxOdd  -Ptr'* .MlnOdd)  <  Phasemln  then 
Phasemin;"  Ptr''.MaxOdd-Ptr''.MinOdd; 

Ptr:"  Ptr''.nextl; 
end;(*whlle  loop*) 

wrltelnCmln  spread  Is  Phasemln: 5:3, '  max  spread  Is  ', 
Phasemax: 5: 3) ; 

end;  (*procedure  MaxlMlnlmlze*) 
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(^•k1t1t1fk1i*1t1tit1i*1i1c1c1t1tit1i1c1i1i1t1fk*1c1tit±*1t1c1c1fk1e1i‘k1c1t1t1t1t1e1c1c1c1c1i1c1c1i1t1i1t1i1c1i1t1i1tii1c1i1i1t1t1t1iic1i1t'^ 

(*  *) 

(*  MAINPROGRAM  *) 

(*  *) 

begin 

FlagA:*  false; 

wrltelnC 'Select  A — for  off  peak  minimax  of  autocorrelation'); 
wrltelnC 'Select  B — for  crosscorrelation  mlnimax  procedure'); 
readln(ch); 

DataList!"  nil; 

Back:*  nil; 
case  ch  of 

'A'  ,  'a'  :begln 

FlagA:*  true; 

wrltelnC 'Autocorrelation  selected — enter  filename  of 
code'); 

readlnC  Fllena 1 ) ; 

CreateLlstCf llenal ^Llstl ,N) ; 

CreateLlstCfllenal ,Llst2,N) ; 

wrltelnC ' enter  filename  where  you  would  like  to  store 
results ' ) ; 
readlnCResult) ; 

ClrScr; 

^  GoToXYC22,14); 

f.  wrlteC 'WORKING  ON  PHASE  SHIFT  OF  '); 

for  J:*  0  to  N-1  do  begin 
GoToXyC48,14); 
wrlteCJ); 

K:-  2; 

S:-  N-1; 

CorrelateCLlstl ,List2,N,K,N) ; 

ShlftLlstCListl); 

ShlftLi8tCLlat2); 
end;C*for  loop*) 
end; 

'B'  ,  'b'  :begin 

FlagA;*  true; 

wrltelnC ' Cross  correlation  selected. .Enter  the  first 
code  filename'); 
readlnC f llenal) ; 

wrltelnC 'Now  enter  the  second  filename'); 
readlnC fllena2) ; 

wrltelnC ' enter  filename  where  you  would  like  to  store 
results'); 
readlnCResult) ; 

CreateLlstCf llenal ,Ll8tl ,N) ; 

CreateLis t  C  f llena2 , List2 ,N) ; 

ClrScr; 

GoToXYC22,14); 

WRITEC 'WORKING  ON  PHASE  SHIFT  OF  '); 
for  J;*  0  to  N-1  do  begin 
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GoToXY(48,14); 
write(J) ; 

K:-  1; 

S:-  N; 

Correlate(Llstl ,List2 ,N,K,  S ) ; 

ShiftListCLlstl) ; 

ShiftList(List2) ; 
end;C*for  loop*) 
end;  (*B  select*) 

end;(*  case*) 
wri t eln ; wrl teln ; 

If  flagA  then  begin 
DlsposeOf (Llst2) ; 

DlsposeOf (Llstl) ; 
asslgn(Res\iltflle, Result) ; 
rewrite (Resultf lie) ; 

Hax1  minimize  (Datallst ,  Phasemar  ,Phaseiiiln,N) ; 

Look : "Datallst ; 

J:*0; 

i^le  Look  O  nil  do  begin 

If  Look'^.MaxOdd  -  Look"' .MlnOdd  ■  Phasemax  then 

wrltelnCResultflle,  'phasemar  *,j,'  Look'*. MlnOdd: 5:2, '  to 
Look'*  .MaxOdd :  5 : 2  )  5 

If  Look'* .MaxOdd  -  Look'*.MlnOdd  ■  Phasemln  then 

wrlteln(Resultflle,  'phasemln  ’,j,'  Look'*. MlnOdd; 5:2, '  to  ', 
Look'*  .MaxOdd :  5 : 2) ; 

Look :  ■Look'* .  nextl ; 

J.—  J+1; 

end(*whlle  loop*); 
close(Resultflie) ; 
end  (*lf  flagA.*) 
end. 


rv 


(^if^iciciticltit’kiticifititititiiitititif’kiiititlcit’kiciticititlc’kicicitlcitlticicitlciciticiticiticifitisiciticitit’kiflcicicitisitlciticicic') 

(*  *) 

(*  CountLower  *) 

(*  *) 

(A************************************************************************) 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  E.6  *) 

(*  DESIGN  DATE;  Jan  1985  LAST  UPDATE:  26  Jan  1985  *) 


(*  DESCRIPTION:  This  procedure  takes  a  data  linked  list  representing  the*) 
(*  the  Odd  and  Even  correlation  functions,  and  determines  what  percent-  *) 
(*  age  of  these  functions  lie  below  the  Inputted  threshold.  *) 


(*  CALLING  ARGUMENTS:  List — pointer  to  front  of  datallst  *) 
(*  Thresh — threshold  value  *) 
(*  Count l,Count2 — used  to  add  amount  of  each  function  *) 
(*  that  exceeds  the  threshold  *) 


(*  INTERNAL  VARIABLES: 

(*  Ptrl,Ptr2 — temp  pointer  used  to  do  the  process 

(*  R — temp  storage  location  used  with  counter 


(*  CONSTANTS;  None 


(*  FILES  USED 


(*  MODULES  CALLED:  None 


procedure  CountLowerCLlst:  dataptr;  Thresh:  real;  var  Countl, 

Count 2:  real); 

var 

Ptrl,  Ptr2:  dataptr; 

R:  real; 

begin 

Countl  :*  0; 

Count2  :■  0; 

Ptrl  List; 
repeat 

Ptr2  :■  Ptrl^'.Nextl; 

If  (-Ptrl'*. Even  <  Thresh)  and  (-Ptr2''.Even  >“  Thresh)  then  begin 
R  (-Ptr2''.Even  -  Thresh)  /  (-Ptr2'‘.Even  +  Ptrl''. Even); 
Countl  ;■  Countl  +  R 
end;  (*lf*) 

If  (-Ptrl''. Even  >■  Thresh)  and  (-Ptr2''.Even  >■  Thresh)  then  begin 
Countl  ;■  Countl  +  1 
end;  (*lf*) 

if  (-Ptrl*. Even  >■  Thresh)  and  (-Ptr2''.Even  <  Thresh)  then  begin 
R  :■  (-Ptrl'*. Even  -  Thresh)  /  (-Ptrl*. Even  +  Ptr2*.Even); 
Countl  :■  Countl  +  R 
end; 


if  CPtrl'‘.Odd  <  Thresh.)  and  (Ptr2'‘.0dd  >“  Thresh)  then  begin 
R  (PtrZ^'.Odd  -  Thresh)  /  (Ptr2''.Odd  -  Ptrl^'.Odd); 

Count2  :■  Count2  +  R 
end;  (*if*) 

if  (Ptrl'‘.Odd  >■  Thresh)  and  (Ptr2'‘.Odd  >“  Thresh)  then  begin 
Count2  ;■  Count2  +  1 
end;  (*if*) 

if  (Ptrl'^.Odd  >*  Thresh)  and  (Ptr2^.0dd  <  Thresh)  then  begin 
R  (Ptrl''.Odd  -  Thresh)  /  <Ptrl''.0dd  -  Ptr2".Odd); 

Count2  :■  Count2  +  R 
end; 

Ptrl  :■  Ptr2; 
until  Ptrl'‘.Nextl  ■  nil 
end;  (*procedure  CountUpper*) 


(*  *) 

(*  CountUpper  *) 

(*  *) 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  E.5  *) 

(*  DESIGN  DATE:  Jan  1985  LAST  UPDATE:  26  Jan  1985  *) 

(* - *) 

(*  DESCILEPTION:  This  procedure  takes  a  data  linked  list  representing  the*) 

(*  the  Even  and  Odd  correlation  functions,  and  determines  what  amount  *) 

(*  of  these  functions  lie  above  the  Inputted  threshold.  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  List — pointer  to  front  of  datallst  *) 

(*  Thresh — threshold  value  *) 

(*  Cotmtl,Count2 — used  to  add  amount  of  curve  that  *) 

(*  exceeds  the  threshold  *) 

(* - *) 

(*  INTERNAL  VARIABLES:  *) 

(*  Ptrl,Ptr2 — temp  pointer  used  to  do  the  process  *) 

(*  R~temp  storage  location  used  with  counter  *) 

(* - *) 

(*  CONSTANTS:  None  *) 

(* - *) 

(*  FILES  USED:  None  *) 

(* - *) 

(*  MODULES  CALLED:  None  *) 

(* - ^ - *) 


procedure  CountUpperCLlst:  dataptr;  Thresh:  real;  var  Countl, 

Count2:  real); 


var 

Ptrl,  Ftr2:  dataptr; 

R:  real; 

begin 

Countl  :■  0; 

Cotmt2  :*  0; 

Ptrl  :■  List; 
repeat 

Ptr2  :■  Ptrl'*.Nextl; 

if  (Ptrl'*. Even  <  Thresh)  and  (Ptr2'‘.Even  >■  Thresh)  then  begin 
R  :■  (Ptr2''.Even  -  Thresh)  /  (Ptr2^.Even  -  Ptrl''.Even); 
Countl  :*  Countl  +  R 
end;  (*if*) 

if  (Ptrl'*. Even  >■  Thresh)  and  (Ptr2'‘.Even  >■  Thresh)  then  begin 
Countl  :•  Countl  +  1 
end;  (*if*) 

If  (Ptrl'*. Even  >“  Thresh)  and  (Ptr2'*.Even  <  Thresh)  then  begin 
R  :■  (Ptrl'*. Even  -  Thresh)  /  (Ptrl''.Even  -  Ptr2''.Even); 
Countl  :•  Countl  +  R 
end; 
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DisAgree  ;■  DisAgree  +1;  j 

Ptrl  :■  Ptrl^.Neact;  i 

Ptr2  ;-  Ptr2'‘.Next 
end;  (*while  Loop*) 

RightCor  :■  Agree  -  DisAgree; 

Agree  ;■  0; 

DisAgree  :■  0; 

Ptrl  ;■  Listl; 
while  Ptr2  O  nil  do  begin 

if  Ptrl^.data  “  Ptr2'‘.data  then 
Agree  ;■  Agree  +  1 
else 

DisAgree  DisAgree  +  1; 

Ptrl  ;■  Ptrl''. Next; 

Ptr2  Ptr2''.Next 
end;  (*iriille  Loop*) 

LeftCor  :*  Agree  -  DisAgree; 

Even  :■  LeftCor  +  RightCor; 

Odd  :■  RightCor  -  LeftCor; 

Enque(DataLlst,  Back,  Odd,  Even); 

K  K  +  1; 
end;  (*whlle  loop*) 

Agree:"  0; 

DisAgree:"  0; 

Ptrl  :"  Listl; 

Ptr2  :"  List2; 

QlT  while  Ptr2  O  nil  do  begin 

If  Ptrl*. data  ■  Ptr2*.data  then 
Agree  :■  Agree  +  1 
else 

DisAgree  :"  DisAgree  +  1; 

Ptrl  :"  Ptrl'*. next; 

Ptr2  :"  Ptr2*.next; 
end;  (*whlle  loop*) 

LeftCor  :"  Agree  -  DisAgree; 

Even  :"  LeftCor; 

Odd  :"  -LeftCor; 

Enque(DataLlst ,  Back,  Odd,  Even) 
end;  (*procedure  correlate*) 
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(*  *) 

(*  Correlate  *) 

(*  *) 

^*************************************1^***********************************) 


(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  E.l  *) 

(*  DESIGN  DATE;  Jan  1985  LAST  UPDATE;  14  Jan  1985  *) 

(* - *) 

(*  DESCRIPTION:  This  procedure  takes  the  2  Linked  Lists  of  characters  *) 

(*  representing  codesets  and  determines  the  Even  and  Odd  correlation  *) 

(*  functions  at  discrete  values.  It  places  the  determined  values  at  *) 

(*  corresponding  positions  In  a  third  linked  list.  *) 

(* - *) 

(*  CAT.TiTNG  ARGUMENTS:  DataLlst ,Back — pointers  to  front  and  rear  of  output*) 
(*  Datallst  *) 

(*  Llatl,Llat2 — pointers  for  front  of  codellsts  *) 

(*  N — Length  of  List  *) 

(* - - - *) 

(*  INTERNAL  VARIABLES:  *) 

(*  Ptrl,Ptr2 — tenqt  pointer  used  to  do  the  process  *) 

(*  LeftCor,  RlghtCor — partial  corelatlons  values  used  to  de-  *) 

(*  termlne  the  Even  and  Odd  values  *) 

(*  Even, Odd, — discrete  correlation  values  *) 

(*  Agree, Disagree — used  as  counters  during  correlation  process  *) 

(*  Count  ,K~counters  used  In  process  *) 

(* - *) 

(*  CONSTANTS:  None  *) 

(* - *) 

(*  FILES  USED:  None  *) 

(* - - *) 

(*  MODULES  CALLED:  Enque  *) 

(* - *) 

procedure  CorrelateCvar  DataLlst,Back:dataptr;Llstl,Llst2:polnter: 


N: Integer); 

var 

LeftCor,  RlghtCor,  Even,  Odd:  real; 

Count,  Agree,  DlsAgree,  K:  Integer; 

Ptrl,  Ptr2;  pointer; 


begin 

K  1; 

while  K  <"  N  do  begin 
Ptrl  ;•  Llstl; 

Ptr2  Llst2; 

Agree  ;■  0; 

DlsAgree  ;■  0; 
for  Count  :■  2  to  K  do 
Ptrl  ;■  Ptrl*. Next; 
while  Ptrl  O  nil  do  begin 

If  Ptrl*. data  ■  Ptr2*.data  then 
Agree  :■  Agree  +  1 
else 


(****************«********************************************************) 
(*  *) 

(*  ShlftLlst  *) 

(*  *) 

(*************************************************************************) 
(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  E.3  *) 

(*  DESIGN  DATE;  Jan  1985  LAST  UPDATE:  14  Jan  1985  *) 

_ _ _ _ * 


(* - 

(*  DESCRIPTION:  This  procedure  merely  takes  the  first  record  of  the 
(*  linked  list  and  places  It  at  the  end  of  the  list. 

(* - 

(*  CALLING  ARGUMENTS;  List — pointers  to  front  of  list 

(* - 

(*  INTERNAL  VARIABLES: 

(*  P,Ptr — temp  pointer 

(* - 

(*  CONSTANTS:  None 


<?ed  to  do  the  process 


(* - 

(*  FILES  USED; 
(*- 


) 

*) 

*) 

-*) 

*) 

-*) 

*) 

*) 

-*) 

*) 

-*) 

*) 

-*) 

*) 

-*) 


None 


(*  MODULES  CALLED; 
(* - 


None 


procedure  ShlftLlstCvar  List:  pointer); 


var 


First,  P:  pointer; 


begin 

First  ;■  List; 

P  List; 

while  P^.Next  O  nil  do 
P  P*.Next; 

First  First''. Next; 
P'‘.Next  ;•  List; 

List''. Next  :■■  nil; 

List  ;■  First 

end;  (^procedure  ShlftLlst*) 


(*  *) 

(*  GetRidOf  *) 

(*  *) 

(******ir***********ik********:k**************ir*********4lt*:k4t***4r***1k*********^ 


(*  DESIGNER:  Btlchard  C.  Gonder 
(*  ADVISOR:  MAJ  Ren  Castor 
(*  DESIGN  DATE:  Jan  1985 


MASTER  THESIS  WORK 
MODULE  NUMBER  E.8 
LAST  UPDATE:  14  Jan  1985 


(*  DESCRIPTION:  This  procedure  releases  a  linked  list  from  memory 
(*  allocation.  Specifically  used  on  Input  data  lists  destiructlon. 


(*  calling  ARGUMENTS:  List — pointers  to  front  of  list 


(*  INTEEINAL  VARIABLES:  P,Ptr — pointers  used  to  do  process 


(*  CONSTANTS:  None 


C*  FILES  USED 


(*  MODULES  CALLED:  None 


procedure  GetRidOf (var  List:  dataptr) 


var  PjPtr;  ^^ataptr; 

begin 

P:-  Ustj 
Ptr:*  P; 

while  P'^.neztl  O  nil  do  begin 
Ptr:-  P''.nextl; 
dlspose(P); 

P:"  Ptr; 

end;  (*whlie  loop*) 
dispose (P); 

List;"  nil; 

end;  (*procedure  GetRidOf*) 


^******************4r****4t***4t4r****it************************4r**4r***ilr*4r*4t**4r^ 

(*  *) 

(*  DlaposeOf  *) 

(*  *) 

(^iiiiii**********************************************************************') 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  E.7  *) 

(*  DESIGN  DATE:  Jan  1985  LAST  UPDATE;  14  Jan  1985  *) 

(* - *) 

(*  DESCRIPTION;  This  procedure  releases  a  linked  list  from  memory  *) 

(*  allocation.  Specifically  used  on  input  data  lists  destruction.  *) 

(* - *) 

(*  CALLING  ARGUMENTS;  List — pointers  to  front  of  list  *) 

(* - *) 

(*  INTERNAL  VARIABLES;  P,Ptr — pointers  used  to  do  process  *) 

(* - *) 

(*  CONSTANTS;  None  *) 

(* - - *) 

(*  PILES  USED:  None  *) 

(* - *) 

(*  MODULES  CALLED:  None  *) 

(* - : - *) 

procedure  DlaposeOf (var  List:  pointer); 

var  P,Ptr;  pointer; 

begin 

P;*  List; 

Ptr;-  P; 

irtille  P^'.nezt  O  nil  do  begin 
Ptr;*  P'^.next; 
dlsposeCP); 

P;*  Ptr; 

end;  (*while  loop*) 
dlsposeCP); 

List:*  nil; 

end;  (*procedure  DlaposeOf*) 


^  I  '.P 


'J  lPP*yp^"»"  i;wg»l!PTPTiit  >t 


whlleC  P  <>  nil)  and  (ABSCP*. thresh  -  thresh)  >  0.0005  )  do  begin 
Ptr:*  P; 

P:*  P^.Next; 
end;  (*whlle  loop*) 

If  P  ■  nil  then  begin  (*record  not  found*) 
new(P); 

Ptr''.Next:-  P; 

If  Number  >  1  then  begin 
P'‘.mlneven:"  0.0; 

P'^.aveeven:*  PctEven/Number; 

P'^.mlnodd:"  0.0; 

P'“.aveodd:"  PctOdd/Number; 

end 

else  begin 

P'^.mlneven:*  PctEven; 

P'*'  .aveeven:“PctEven; 
p'‘.mlnodd:“  PctOdd; 

P** .  aveodd ;  “PctOdd ; 
end;(*lf  Number  then  else*) 

P'^.maxeven:*  PctEven; 

P‘*.maxodd;“  PctOdd; 

P'‘. thresh:"  thresh; 

P'.Next  nil 
end  (*front  clause  of  If*) 
else  begin 
Ptr:“  P; 

If  PctEven  <  P'^.mlneven  then  P*.mlneven:“  PctEven; 

If  PctEven  >  P^'.maxeven  then  P''.maxeven:“  PctEven; 

P'^.aveeven  ((P'^.aveeven  *  (Number-1)  +  PctEven) /Number) ; 

If  PctOdd  <  P''.mlnodd  then  P''.mlnodd:“  PctOdd; 

If  PctOdd  >  P'‘.maxodd  then  P''.maxodd:“  PctOdd; 

P"*, aveodd  ((P''. aveodd  *  (Number-1)  +  PctOdd) /Number) ; 

end  (*lf  then  else*) 
end  (*lf  then  else*) 
end;  (*UpdateData*) 


(*  *) 

(*  UpdateDataList  *) 

(*  *) 

^**4r4t****4t**it********W1k*4r1lr*****************4c4t4t*4;*4t**4r4r*******4t*1k4r*4t4r4c**4t**^ 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  E.IO  *) 

(*  DESIGN  DATE:  Mar  1985  LAST  UPDATE:  4  Apr  1985  *) 


C’' 


0 


(*  DESCRIPTION:  This  procedure  takes  the  values  of  Threshold,  PctEven  *) 
(*  and  PctOdd  from  the  thresholding  process.  It  finds  the  record  that  *) 
has  stored  Minimum,  Maximum  and  Average  values  for  a  specific  thresh-*) 
old  value,  for  both  the  Odd  and  Even  functions,  and  updates  these  *) 
fields  with  the  appropriate  Info.  If  a  record  Is  not  found.  It  ere-  *) 
ates  the  record  and  places  appropriate  Intlal  values  In  the  fields.  *) 


(* 

(* 

(* 

(* 


(*  CALLING  ARGUMENTS:  List — pointer  to  front  of  list 


(* 

(* 

(* 

(* 


Even, Odd — real  values  of  correlation  functions 
Thresh — value  of  threshold 

Number — ^Number  of  thresholding  processes  already 
accomplished 


*) 

*) 

*) 

*) 

*) 


0 


(*  INTERNAL  VARIABLES; 

(*  Ptr,P — temp  pointers  used  to  do  the  process 


*) 

*) 


(*- 


0- 


(*  CONSTANTS:  None 
(* - 


*) 

-*) 

*) 


(*  PILES  USED;  None 


0 


(*  MODULES  CALLED: 
(* - 


None 


*) 

-*) 


procedure  UpdateDataLlstCvar  List  :  Pctptr;  Number: Integer;  PctEven, 

PctOdd,  Thresh:  real): 


var 

Ptr,P;  Pctptr; 

begin 

P;"  List; 

If  List  ■  nil  then  begin 
new(P); 

List  ;■  P; 

P'‘.mlneven;“  PctEven; 
P'^.maxeven;"  PctEven; 
P^.aveeven:"  PctEven; 
P'^.mlnodd;"  PctOdd; 
P'^.maxodd:"  PctOdd; 
P''.aveodd;"  PctOdd; 
P'*. thresh:"  thresh; 
P^'.Next  nil 
end  (*flrst  clause*) 
else  begin 


(*  *) 

(*  E  n  q  u  e  *) 

(*  *) 

(******«******************************************************************) 

(*  DESIGNER;  Richard  G.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  E.ll  *) 

(*  DESIGN  DATE:  Nov  1984  LAST  UPDATE:  14  Dec  1984  *) 


(*  DESCRIPTION;  This  procedure  enques  a  data  record  representing  the 
(*  values  of  the  Odd  and  Even  correlation  functions  at  a  discrete 
(*  value  in  an  ordered  linked  list. 


(*  CALLING  ARGUMENTS:  List, Back — pointers  to  front  and  rear  of  list 
(*  Even, Odd — values  of  correlation  functions 


(*  INTERNAL  VARIABLES: 

(*  p — temp  pointer  used  to  do  the  process 


(*  CONSTANTS:  None 


(*  FILES  USED 


(*  MODULES  CALLED 


procedure  EnqueCvar  List,  Back:  dataptr;  Odd,  Even:  real) 


P:  dataptr; 


begin 

if  List  "  nil  then  begin 
new(P); 

List  ;■  P; 

P'^.Even  ;■  Even; 
P^'.Odd  ;■  Odd; 

Back  :■  P; 

Back''. Next  1  nil 

end  else  begin 
new(P); 

P'*.Even  ;■  Even; 
P^'.Odd  :-  Odd; 
Back'‘.Nextl  :■  P; 

Back  ;■  P; 

Back^.Nextl  :■  nil 
end  (*lf  then  else*) 
end;  (*procedure  Enque*) 


Qi^lt*******ifrir1lr***il^*!l^**:fc*******!(lr*1lr*****i(lr*1lr******4r*^*!l^*1lr****4rir***4r4r**********^ 

(*  *) 

(*  GreateLlst  *) 

(*  *) 

(******************1k***************lk*********4t******;k1lc*ilr******************^ 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  E.2  *) 

(*  DESIGN  DATE;  Jan  1985  LAST  UPDATE;  4  Jan  1985  *) 


(*  DESCRIPTION:  This  procedure  will  read  in  a  file  of  char,  placing  each*) 
(*  character  into  a  linked  list.  It  terminates  at  the  EOF  of  the  file  *) 
(*  read.  *) 


(*  CALLING  ARGUMENTS:  Filename — string  of  char 

(*  List — pointer  to  the  beginning  of  linked  list 

(*  K — number  of  characters  read  into  list 


(*  INTERNAL  VARIABLES: 

(*  PiPtr — used  to  construct  list 

(*  bit — temp  location  for  character  read 


(*  CONSTANTS;  None 


(*  FILES  USED:  User  designated  Filename  for  code  storage 


(*  MODULES  CALLED: 


procedure  CreateLlst( Filename;  Name;  var  List:  pointer;  var  K:  Integer) 


P,  ptr:  pointer; 
bit:  char; 


begin 

asslgnCCodedta, Filename) ; 
reset (Coded ta) ; 
new(P) ; 

List  ;■  P; 
ptr  ;-  P; 

K  0; 

idille  not  eof(Codedta)  do  begin 
readCCodedta,  bit); 

K  ;-  K  +  1; 

P'‘.data  :■  bit; 
ptr''. Next  :■  P; 
ptr  ;■  P; 
new(P) 

end;  (*while  loop*) 
close ( Coded ta); 
dlspose(P); 
ptr''. Next;*  nil; 
end;  (*procedure  CreateLlst*) 


program  MassCorrelateC Input,  output); 
type 

Name  ■  strlng[14]; 

St3  ■  8tringt3]; 

StlO  ■  string[lO]; 
pointer  ■  datarecord; 
datarecord  ■ 
record 

data:  char; 

Next:  pointer 
end;  (*of  record*) 
dataptr  ■  Info; 
info  ■ 
record 

Even;  teal; 

Odd:  real; 

Nextl:  dataptr 
end;  (*record*) 

Pctptr  ■  Infol; 

Infol  " 
record 

thresh:  real; 
mlneven:  real; 
mazeven:  real; 
aveeven:  real; 
mlnodd:  real; 
mazodd:  real; 
aveodd:  real; 

Next:  Pctptr 
end;  (*record*) 


Number , Outside , Inside , S tartf lie , S topFlle , Phase , J , N :  Integer 
Llstl,  Llst2:  pointer; 

DstaLlst,  Back:  dataptr; 

PctLlst;  Pctptr; 

Result,  Fllenal,  Fllena2:  Name; 

Stl,St2:  St3; 

CodeName:  StlO; 

Resultfile,  Codedta:  text; 

PctEven,  PctOdd,  S, 

Sumleven,  Sum2even,  Sumlodd,  Sum2odd,  Threshold:  real; 
ch:  char; 

Stop: boolean; 


(*  *) 

(*  MassCorrelate  *) 

(*  *) 

^***4k*********************1lt**4t*************4r******************************^ 

(*  DESIGNER;  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR;  MAJ  Ken  Castor  MODULE  NUMBER  E.O  *) 

(*  DESIGN  DATE;  Jan  1985  LAST  UPDATE;  1  Apr  1985  *) 

(* - *) 


(*  DESCRIPTION;  This  program  will  perform  correlations  on  user  designs-  *) 
(*  ted  codesets  consisting  of  files  composed  of  ASCII  I's  and  O's.  The  *) 
i*  Codeset  Is  Identified  by  the  user,  and  bounded  by  the  extensions  l.e.*) 


(*  NAME. 001  NAME. 005.  If  these  are  the  user  Inputted  filenames,  the  *) 
(.*  program  will  do  the  cross  correlations  of  all  pairs  of  codes  in  this  *) 
(*  codeset.  It  determines  the  pet  out  of  threshold  function  for  each  *) 
(*  pair.  In  terms  of  the  even  and  odd  correlation  function  and  stores  *) 
(*  the  minimum,  maximum  and  average  value  for  each  threshold  value.  *) 

(*  When  the  process  Is  completed,  these  values  are  written  to  the  user  *) 
(*  Identified  output  file.  *) 

(* - *) 

(*  CALLING  ARGUMENTS;  Std  Input,  Output,  Designated  Text  File  *) 

(* - *) 

(*  INTERNAL  VARIABLES;  *) 

(*  Phase — the  phase  shift  of  desired  process  *) 

(*  PctEven,PctOdd — used  for  output  of  thresholding  values  *) 

(*  Sumleven,Sum2even,Suml0dd,Sum20dd--temp  storages  values  for  *) 

(*  thresholding  process  *) 

(*  Threshold — actual  threshold  passed  Into  threshold  process  *) 

(*  DataLlst, Back— pointers  used  to  create/manipulate  the  thresh  *) 

(*  old  functions  *) 

(*  Resultflle,Codedta — ^Flle  variable  for  file  I/O  *) 

(*  . Result, Fllenal,Fllena2;  string  for  filename  with  extension  *) 

(*  CodeName — User-Inputted  filename  for  storage  of  the  code  *) 

(*  ch — character  Input  to  the  program  from  user  *) 

(*  Llstl,Llst2  used  to  Identify  the  beginning  of  linked  lists  *) 

(*  representing  the  codes  to  be  correlated  *) 

(*  PctLlst — used  to  Identify  beginning  of  compiled  output  data  *) 

(*  Stop — boolean  used  to  stop  correlation  process  *) 

(*  Stl,St2 — string  expressions  for  filename  extensions  *) 

(*  Startflle,Stopflle — used  to  store  Integer  value  of  filename  *) 

(*  extensions  which  bound  the  entire  process  *) 

(*  Number — represents  the  number  of  cycles  the  program  has  ac-  *) 

(*  pllshed  and  Is  uaed  to  determine  the  average  values  *) 

(*  Inside, Outside — used  to  control  looping  structures  *) 

(* - ^ - *) 

(*  CONSTANTS;  None  *) 

(* - *) 

(*  FILES  USED;  User  designated  Filename  for  code  readln  and  data  output  *) 
(* - *) 


(*  MODULES  CALLED; 

CreateLlst 

,  UpdateDataLlst,  GetRldOf 

*) 

(* 

ShlftLlst, 

Correlate,  DetermlneName 

*) 

(* 

Wrltellst , 

CountUpper,  FlxLlst,  CountLower, 

DlsposeOf*) 

(* - 

- *) 

if  (-Ptrl''.Odd  <  Thresh)  end  (~Ptr2'‘.0dd  >■  Thresh)  then  begin 
R  (-Ptr2''.0dd  -  Thresh)  /  (-Ptr2'‘.0dd  +  Ptrl^.Odd); 

Count2  ;■  Count2  +  R 
end;  (*if*) 

if  (-Ptrl''.Odd  >■  Thresh)  and  (-Ptr2''.0dd  >■  Thresh)  then  begin 
Count2  Count2  +  1 
end;  (*if*) 

if  (-Ptrl''.Odd  >“  Thresh)  and  (-Ptr2'‘.0dd  <  Thresh)  then  begin 
R  :■  (-Ptrl''.Odd  -  Thresh)  /  (-Ptrl''.Odd  +  Ptr2''.Odd); 

Count2  ;■  Count2  +  R 
end; 

Ptrl  Ptr2; 

until  Ptrl''.Nextl  "  nil 


end;  (*procedure  CountLower*) 


(*  *) 

(*  Writelist  *) 

(*  *) 

(************************«*****«******************************************) 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  E.4  *) 

(*  DESIGN  DATE:  Nov  1984  LAST  UPDATE:  14  Mar  1985  *) 


(*  DESCRIPTION:  This  procedure  takes  a  linked  list  of  data  and  writes 
(*  the  data  fields,  mlneven,  mlnodd,  mazeven,  mazodd,  aveeven,  aveodd 
(*  and  Threshold,  to  the  designated  file. 


(*  CALLING  ARGUMENTS:  List — pointer  to  front  of  datallst 
(*  OutFlle-strlng  for  output  filename 


(*  INTERNAL  VARIABLES: 

(*  P — temporary  pointer  used  to  traverse  list 

(*  K — Counter  used  as  chip  position  In  list 


(*  CONSTANTS:  None 


(*  FILES  USED:  Writes  to  OutFlle 


(*  MODULES  CAT.T.Rn 


procedure  Writelist (List:  Pctptr;  var  OutFlle:  tezt); 


K:  Integer; 

P:  Pctptr; 

begin 

P  :-  List; 

GoToXY(20,l4); 

wrltelnC ’WRITING  THE  OUTPUT  FILE  NOW  ’); 

while  P  O  nil  do  begin 

wrlte( OutFlle,  F".thre8h:5:4, ' , ' ,P''.mineven:  8:  7, 

P''. mazeven:  8:  7,',',  P''. aveeven:  8:  7); 
wrltelnC  OutFlle,  ’ , '  ,P''.niinodd:  8:  7,  P''. mazodd:  8;  7 

P". aveodd:  8:  7); 

P  :-  P^'.Nezt 
end  (*whlle  loop*) 
end;  (*procedure  Writelist*) 


(*  *) 

(*  DetermineName  *) 

(*  *) 

(A************************************************************************) 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  E.9  *) 

(*  DESIGN  DATE:  Mar  1985  LAST  UPDATE:  14  Mar  1985  *) 

(* - *) 


(*  DESCRIPTION:  This  procedure  takes  two  strings  representing  filenames  *) 
(*  with  extensions,  and  separates  the  extensions  from  the  filenames .  It  *) 
(*  then  converst  the  extension  strings  to  Integer  values,  Startflle  and  *) 
(*  Stopflle.  It  should  be  noted  that  the  function  calls  In  this  module  *) 
(*  are  not  standard  and  the  Module  may  need  to  be  modified  for  different*) 


(*  systems.  *) 

(* - *) 

(*  CALLING  ARGUMENTS:  Fllenal,Fllena2 — strings  representing  filenames  *) 
(*  Startflle, Stopflle — Integer  value  of  filename  ex-  *) 

(*  tensions  *) 

(* - *) 

(*  INTERNAL  VARIABLES:  *) 

(*  Stl,St2 — strings  representing  filename  extensions  *) 

(*  K,I — temp  storage  for  character  positions  In  strings  *) 

(* - *) 

(*  CONSTANTS:  None  *) 

(* - *0 

<*  FILES  USED:  None  ♦) 

(* - *) 

(*  MODULES  CALLED:  None  *) 

(* - *) 


procedure  DetermlneName(var  Fllenal,Fllena2;Name;  var  Startflle, 

StopFlle: Integer;  var  CodeName:stlO) ; 

var  K,I:lnteger; 

Stl,St2:  St3; 


begin 

I:*  pos( ' . ' ,fllenal) ; 

Stl:“  copy (Fllenal,I+l, 1+4) ; 
St2:"  copy(Fllena2,I+l,I+4); 
CodeName:"  copy(Fllenal,l,I); 
val(Stl,Startflle,K) ; 
val(St2 , Stopflle ,K) ; 
end;  (*procedure  DetermineName*) 


(*  *) 

(*  FlxLlst  *) 

(*  *) 

(lllllcll*********************************************************************^ 

(*  DESIGNER:  Richard  C.  Gonder  MASTER  THESIS  WORK  *) 

(*  ADVISOR:  MAJ  Ken  Castor  MODULE  NUMBER  E.10.1  *) 

(*  DESIGN  DATE:  Mar  1985  LAST  UPDATE:  4  Apr  1985  *) 


(*  DESCRIPTION:  This  procedure  takes  a  linked  list  of  mlnlmums  mafimumg  *) 
(*  and  averages,  after  a  thresholding  process  has  been  completed  on  one*) 
(*  set  of  even  and  odd  fvmctlons,  and  places  appropriate  values  In  the  *) 
(*  records  that  were  not  updated  during  the  the  thresholding  process  .  *) 


(*  CALLING  ARGUMENTS:  List — pointer  to  front  of  list  *) 
(*  Thresh — value  of  threshold  *) 
(*  Number — ^Number  of  thresholding  processes  already  *) 
(*  accomplished  *) 


(*  INTERNAL  VARIABLES: 

(*  P — temp  pointers  used  to  do  the  process 


(*  CONSTANTS:  None 


(*  FILES  USED:  None 


(*  MODULES  CALLED:  None 


procedure  FlxLlst (Thresh: real;  Number: Integer;  Llst:Pctptr) 


var  P;  Pctptr; 


begin 

P  :■  list; 

while  (P  O  nil)  and  (P''. thresh  O  Thresh)  do 
P  :■  P''.next; 
while  P  O  nil  do  begin 

P'',aveodd  :■  P^.aveodd  *  (Ntmber  -  1)  /  Number; 
P'^.mlnodd  :■  0.0; 

P''.aveeven  :■  P'‘.aveeven  *  (Number  -  1)  /  Number; 
P'‘.mlneven  :■  0.0; 

P  :■  P''.next; 
end;  (*whlle  loop*) 
end;  (*procedure  FlxLlst*) 


(*************4t**********4r**4r**************ir*4r4t**********4t****4t************) 

(*  *) 

(*  MAINPROGRAM  *) 

(*  *) 

(********i[r***4t****4t**********1lr***4t4t**********4t*******4r**********4t**********^ 


M  « 


I  .V 

V,  \V, 


begin 

PctLi8t;*nil; 

wrlteln( 'This  program  allows  you  to  do  correlations  of  ASCII  code  files'); 
wrltelnC 'representing  a  set  of  pseudorandom  sequences.  You  Input  only'); 
wrltelnC'the  first  and  last  codeflle  In  the  set  Identified  by  Its 
extension'); 

wrltelnC 'and  the  program  does  the  crosscorelatlon  and  threshholdlng ' ) ; 
wrltelnC 'procedure  on  each  pair.  It  will  determine  the  minimum,'); 
wrltelnC 'maximum  and  average  values  of  the  Threshold  process  and  writes 
the'); 

wrltelnC 'results  to  the  desired  file'); 
wrlteln; 

wrltelnC 'Do  you  want  to  do  any  correlations  at  this  time? . .7/N' ); 
readlnCch); 

If  Cch  ■  'y')  or  Cch  ■  'Y')  then  begin 
repeat 

wrltelnC 'enter  the  first  file  you  wish  to  correlate'); 
wrltelnC 'with  extension  ..  l.e.  .001,  .002  etc'); 
readlnCFllenal ) ; 

wrltelnC 'now  enter  the  second  Filename  w/extentslon' ) ; 
readlnCFllena2) ; 

wrltelnC 'Enter  phaaeshlft  desired  for  run. ..0  to  length'); 
readlnCPhase) ; 

wrltelnC ' enter  the  file  to  write  results'); 
readlnCResult ) ; 

DetermlneNameCFllenal ,Fllena2,Startflle,StopFlle,CodeName) ; 

Number:"  1; 

ClrScr; 

for  Outside:"  CStartflle)  to  CStopflle  -1)  do  begin 
str COutslde : 3 , St 1 ) ; 

If  stl[l]  "  '  '  then  stl[ll:"  'O'; 

If  stl[2]  "  '  '  then  8tl[2l:"  *0'; 

Fllenal:"  CodeName  +  stl; 

CreateLlstCFllenal,  Llstl,  N); 
for  J  :"  1  to  Phase  do 
ShlftLlstCLlstl); 

for  Inside:"  COutslde  +  1)  to  Stopflle  do  begin 
str C Inside : 3 , st2) ; 

If  8t2[l]  "  '  '  then  st2[l]:"  'O'; 

If  8t2[2]  "  '  '  then  8t2[2l:"  'O'; 

Fllena2:"  CodeName  +  st2; 

CreateLlstCFllenal,  Llstl,  N); 
for  J  :"  1  to  Phase  do 
ShlftLlstCLlstl); 

DataLlst  :"  nil; 

Back  :"  nil; 

CorrelateCDataLlst ,  Back,  Llstl,  List!,  N); 


t'V  t  ^  ^  Jf  V  •  L  •  •  V 


Threshold  :■  0; 
stop  ;■  false; 

GoToXY(20,14); 

writelnC 'WORKING  ON  '.fllenal,'  and  ’.fllenaZ); 
Sumleven  ;■  0; 

Sum2even  :■  0; 

Sumlodd  :■  0; 

Sum2odd  :■  0; 

while  (Threshold  <■  N)  and  not  stop  do  begin 
CountUpper(DataList,  Threshold, Sumleven, Sumlodd) ; 
CotmtLower(DataLlst ,  Threshold, Sum2even, Sum2odd) ; 
PctEven  ;*  (Sumleven  +  Sum2even)  /  N; 

PctOdd  :■  (Sumlodd  +  Sum2odd)  /  N; 

S:-  Threshold/N; 

UpdateDataLlst(PctLlst,  Number,  PctEven,  PctOdd,  S) 
if  Threshold  /  N  >  0.10  then 

Threshold  Threshold  +  N  /  100 

else 

Threshold  :•  Threshold  +  N  /  500; 
if  (PctEven  ■  0.0)  and  (PctOdd  “  0.0)  then 
stop  :■  true; 
end(*while  loop*); 

Fizlii  St  (Threshold/N ,  Number ,  Pc tliist ) ; 

GetRidOf (DataList) ; 

DisposeOf (List2) ; 

Number:"  Number  +  1; 
end  (*foi;  Inside  loop*); 

DisposeOf (Listl) : 
end;(*for  Outside  Loop*) 
asslgn(Resultf lie, Result) ; 
rewrlte(Resultflle) ; 

Wrltellst(Pctllst,Resultflle) ; 
close(Resultflle) ; 

Mark(Pctlist); 

Release(PctLlst) ; 

ClrScr; 

wrltelnCDo  you  want  to  do  any  other  codesets?  7/N'); 
readln(ch); 

tmtll  (ch  ■  'N')  or  (ch  ■  'n'); 
end(*lf*) 

end.  {  PartlalCorrelate  } 
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